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ABSTRACT 


Design  guidelines  were  developed  to  provide  a  basis  for  conducting  design  trades 
for  a  homing  type  air  to  surface  missile  (ASM)  with  high  terminal  accuracy. 

Three  basic  homing  guidance  concepts:  proportional,  pursuit,  and  optimal  guid¬ 
ance  were  evaluated  on  the  basis  of  impact  error.  Two  nominal  trajectories  were 
investigated. 

An  optimal  guidance  law  was  developed  for  an  ASM  with  realistic  aerodynamic 
and  sensor  characteristics.  This  guidance  law  was  based  on  the  use  of  a  Kalman 
filter  to  obtain  best  estimates  of  the  ASM  state  variable  errors,  and  a  control 
criterion  that  minimizes  the  sum  of  the  mean  square  impact  error  and  the  integral 
of  a  quadratic  form  of  the  autopilot  control  variables. 

A  linearized  differential  equation  program  that  computed  the  mean  square  impact 
error  in  the  form  of  a  covariance  matrix  deviation  perpendicular  to  the  nominal 
trajectory  was  used  for  comparison  of  the  guidance  laws. 

A  normal  acceleration  autopilot  was  designed  to  meet  the  mission  requirements, 
and  advanced  bistable  controller  techniques  were  applied  to  obtain  a  quasiadaptive 
autopilot  that  required  no  gain  changes  throughout  the  ASM  midcourse  and  terminal 
phases. 

A  limited  state-of-the-art  survey  was  conducted  of  homing  and  inertial  sensors 
and  on-board  digital  computers  suitable  for  a  homing  ASM. 
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SECTION  I 


INTRODUCTION 

Recent  trends  in  requirements  for  high-accuracy  air-to-surface  missiles  (ASM) 
indicate  a  need  for  development  of  advanced  control  systems.  Even  with  perfect 
guidance  information,  the  control  task  is  formidable.  Three  major  control  areas 
and  questions  of  feasibility  with  regard  to  high-accuracy  ASM's  are  of  immediate 
concern  to  the  control  system  designer: 

•  Flight  Path  Control:  Are  conventional  techniques  such  as  pursuit  (veloc¬ 
ity  vector  aimed  toward  target)  and  proportional  (normal  acceleration 
proportional  to  the  rate  of  change  of  the  line  of  sight)  adequate  for  high- 
accuracy  systems,  or  are  new  techniques  required? 

•  Guidance  Sensors;  What  are  critical  sensor  characteristics;  what  are 
the  effects  of  bias  errors  and  output  noise?  What  are  the  best  methods 
for  guidance  during  the  final  phase  of  flight  where  the  guidance  sensor 
no  longer  provides  information  (i.e.,  the  sensor  blind  range)? 

•  Inner -Loop  Control:  Is  the  performance  of  such  conventional  control 
techniques  as  normal  acceleration  autopilots  adequate  in  the  presence 
of  rapidly  changing  flight  conditions  and  atmospheric  disturbances? 

The  primary  objective  of  this  program  was  to  provide  answers  to  some  of  these 
questions  in  the  form  of  control  system  design  guides  for  ASM's  with  high  termi¬ 
nal  accuracy.  It  was  intended  that  these  design  guides  would  delineate  those  sub¬ 
systems  and  components  that  limit  the  attainable  accuracy  of  the  ASM. 

A  secondary  objective  was  to  test  the  practicality  of  using  optimal  control  theory 
and  techniques  in  solving  ASM  control  problems  when  practical  missile  configur¬ 
ations  and  realistic  control  subsystems  characteristics  are  specified. 


I 


SECTION  n 


SUMMARY 

Design  guidelines  were  developed  to  provide  a  basis  for  conducting  the  following 
trades  for  a  homing  missile  control  system: 

(1)  Selection  of  the  most  suitable  guidance  concept  and  best  gain  for  that  concept; 

(2)  Selection  of  homing  and  inertial  sensors; 

(3)  Selection  of  autopilot. 

Three  basic  homing  guidance  concepts  — proportional,  pursuit,  and  optimal 
guidance  —  were  evaluated  on  the  basis  of  miss  distance.  Two  nominal  trajec¬ 
tories:  a  long-range  semiballistic  and  a  low-altitude  skip  trajectory,  were  in¬ 
vestigated  for  two  homing  sensor  acquisition  slant  ranges. 

In  addition  to  the  design  guides,  an  optimal  guidance  law  was  developed  based  on 
the  use  of  a  Kalman  filter  to  obtain  the  best  estimates  of  the  ASM  state  variable 
errors  from  imperfect  sensor  information.  This  guidance  law  minimizes  a 
performance  index  that  is  the  sum  of  the  mean  square  miss  distance,  plus  the 
integral  of  a  quadratic  form  of  the  autopilot  control  variables.  The  relative 
weighting  of  the  two  terms  is  adjusted  by  a  weighting  factor  on  the  control  integral 
term.  This  provides  a  means  for  minimizing  the  miss  distance  with  realistic 
restraints  on  control  action. 

A  normal  acceleration  autopilot  was  designed  to  meet  the  mission  requirements 
and  advanced  bistable  controller  techniques  were  applied  to  obtain  an  autopilot 
with  nearly  invariant  performance  without  gain  changes  throughout  the  ASM  mid¬ 
course  and  terminal  phases. 

Evaluation  studies  were  conducted  using  two  analytical  tools:  a  3-degree-of- 
freedom  linearized  differential  equation  digital  program  and  a  6-degree-of- 
freedcm  analog  simulation.  The  digital  program  provided  a  direct  computation 
of  the  mean  square  miBS  distance  for  optimal,  proportional,  and  pursuit  guidance 
laws.  The  effects  of  missile  nonlinearities  and  short-period  control  dynamics 
were  not  considered  in  the  digital  program,  but  tracker  dynamics,  tracker  errors 
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(noise  and  bias),  and  uncertainties  in  initial  position  and  velocity  due  to  boost 
and  midcourse  inertial  measurement  errors  were  included.  The  analog  simula¬ 
tion  was  used  primarily  to  provide  a  means  of  determining  the  effects  of  guidance 
gain,  missile  aerodynamic  nonlinearities,  atmospheric  turbulence,  and  autopilot 
characteristics;  however,  it  was  also  used  to  obtain  additional  data  for  compari¬ 
son  of  the  proportional  and  pursuit  guidance.  Quasistatistical  miss  distance  data 
were  obtained  from  the  analog  simulation  by  conducting  several  runs  and  calculat¬ 
ing  the  rms  for  each  data  point  when  stochastic  disturbances  were  Investigated. 
Additional  backup  studies  in  the  form  of  root-locus  stability  studies,  time  and 
frequency  response  analyses,  and  adaptive  autopilot  studies  provided  necessary 
autopilot  design  data. 

The  following  information  is  presented: 

(1)  Trades  for  selection  of  guidance  concept; 

(a)  Comparison  of  miss  distance  for  optimal,  proportional,  and  pursuit  guidance, 

(b)  Effect  of  guidance  gain  on  miss  distance, 

(c)  Effect  of  nominal  trajectory  on  miss  distance, 

(d)  Effect  of  atmospheric  turbulence  on  miss  distance, 

(2)  Trades  for  selection  of  sensors  (inertial  and  homing); 

(a)  Effect  of  dispersion  at  acquisition  (inertial  measurement  unit  error,  IMU 
error  vs  cost), 

(b)  Effect  of  homing  sensor  bias  errors  and  output  noise, 

(c)  Effect  of  homing  sensor  blind  range, 

(3)  Trades  for  selection  of  autopilot; 

(a)  Effect  of  short  period  control  dynamics  on  miss  distance, 

(b)  Achievement  of  frequency  and  time  response, 

(c)  Achievement  of  stability. 

The  results  of  the  guidance  law  comparisons  indicate  that  all  three  of  the  guidance 
concepts  are  capable  of  la  miss  distance  of  less  than  5  feet  with  nominal  TV 
tracker  bias  errors  of  1.74  mllllradlans  (0.1  degree)  and  output  noise  of  1.0 
mllliradlan  rms.  These  small  miss  distances  can  be  achieved  over  a  reasonably 
wide  range  of  guidance  gain,  which  indicates  some  latitude  in  the  implementation  of 
guldanee/autopllot  loops. 
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Miss  distance  was  found  to  be  essentially  invariant  with  both  homer  acquisition 
range  and  nominal  trajectory  type.  Note  that  miss  distance  was  measured  in  a 
plane  perpendicular  to  the  trajectory,  and  therefore,  miss  distances  for  the  dif¬ 
ferent  trajectories  will  be  different  when  resolved  into  the  ground  plane.  For 
example,  because  of  the  shallow  approach  angle  of  the  low-altitude  skip  trajectory, 
the  miss  distances  in  a  direction  along  the  trajectory  are  approximately  ten  times 
those  of  the  semiballistic  trajectories  that  have  more  nearly  vertical  approach 
angles. 

Initial  dispersion,  as  a  result  of  inertial  sensor  errors  generated  prior  to  target 
acquisition,  appears  to  have  little  effect  on  miss  distance.  However,  it  is  pos¬ 
sible  to  have  a  combination  of  large  initial  dispersion  and  short  slant  range  at  ac¬ 
quisition  that  would  affect  miss  distance  (the  shortest  slant  range  investigated 
was  15, 000  feet).  Initial  errors  were  essentially  nulled  and  the  missile  velocity 
vector  aligned  with  the  target  line  of  sight  with  2000  to  5000  feet  of  slant  range 
remaining  prior  to  impact.  This  indicates  that  the  most  important  aspect  of 
initial  dispersion  will  probably  be  in  considering  trades  between  inertial  guidance 
sensor  accuracy  and  the  acquisition  capability  of  the  homing  sensor.  The  effects 
of  atmospheric  turbulence  and  gusts  were  also  found  to  be  insignificant.  This 
was  not  unexpected,  because  the  high  missile  velocity  (1000  -  2500  fps)  combined 
with  relatively  low  gust  velocities  (<  12  fps  rms)  result  in  very  small  angle-of- 
attack  perturbations.  In  addition,  the  guidance  loop  tends  to  nullify  the  effect  of 
these  perturbations.  However,  the  effect  of  the  turbulence  does  become  significant 
for  large  tracker  blind  ranges  (greater  than  1500  feet). 

The  best  estimates  of  homing  sensor  (TV  tracker)  nominal  bias  and  noise  level 
(1,74  milliradians  and  1.0  mtlliradian  rms,  respectively)  had  a  negligible  effect 
on  miss  distance.  An  order  of  magnitude  increase  in  the  nominal  noise  level 
increases  the  miss  distance  considerably.  Pursuit  guidance  miss  distances  are 
proportional  to  the  bias  errors.  Tracker  blind  range  has  very  little  effect  up  to 
about  1500  feet.  Returning  the  acceleration  command  to  the  autopilot  to  zero 
during  blind  range,  rather  than  holding  the  last  guidance  sensor  command  value, 
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reduces  the  dispersion  for  large  tracker  blind  ranges.  For  proportional  guidance, 
the  allowable  minimum  slant  range  to  the  target  is  related  to  the  overall  autopilot/ 
guidance  stability  (and  thus  directly  to  the  missile  short-period  response). 

Controlled  missile  short-period  response  had  a  definite  effect  on  miss  distances. 
With  the  detailed  autopilot  configuration  investigated,  a  sharp  increase  in  miss 
distance  was  noted  when  the  equivalent  first-order  response  was  lowered  below 
6  radians  per  second.  The  increase  in  error  was  caused  by  an  inner-loop  sta¬ 
bility  problem  rather  than  a  'looseness"  of  control  problem.  Earlier  more  simpli¬ 
fied  studies  using  simplified  third  order  transfer  function  approximations  of  con¬ 
trolled  missile  response  indicated  that  2  radians  per  second  would  be  adequate . 

The  results  obtained  with  more  detailed  analog  simulation  emphasize  the  neces¬ 
sity  of  considering  guidance  gains,  autopilot  loops  and  gains,  and  body  bending 
together  when  the  homing  guidance  problem  is  investigated. 

Digital  computer  programs  were  developed  to  determine  die  coefficients  of  the 
optimal  feedback  control  law  and  the  Kalman  filter  through  the  solutions  of  Equa¬ 
tions  53  and  54  and  Equations  76  and  77 ,  respectively.  Programs  were  also  de¬ 
veloped  for  determining  the  mean  square  miss  distance  for  an  ASM  with  optimal, 
proportional,  or  pursuit  guidance.  These  programs  are  available  in  punched  card, 
program  listing,  or  magnetic  tape  form;  however,  no  user's  documentation  is 
available.  Requests  for  information  regarding  the  above  programs  should  be 
directed  to  L.  E.  Leistikow  of  The  Boeing  Company,  Missile  and  Information  Sys¬ 
tems  Division.  The  coefficient  matrices  (F(t)  and  G(t)  of  Equation  49)  required  as 
inputs  to  the  optimal  feedback  control  law  and  Kalman  filter  programs  were 
determined  with  an  existing  program  that  is  documented  in  Reference  14. 

Discussions  of  the  details  of  conditions  investigated  and  ihe  results  obtained  are 
presented  in  Ihe  following  sections. 
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SECTION  m 


GENERAL  DESCRIPTION  OF  STUDY 

The  program  was  aimed  at  effort  prior  to  the  preliminary  design  stage.  That  is, 
the  first  task  was  to  determine  whether  control  of  a  highly  accurate  ASM  is  even 
possible  considering  the  various  "real  world"  limitations  that  must  be  imposed. 
When  feasibility  was  established,  the  task  was  to  develop  design  guides  that  could 
be  applied  to  the  design  of  control  systems  for  this  general  class  of  vehicles.  The 
study  consisted  of  four  ipajor  efforts: 

(1)  A  limited  survey  was  conducted  of  the  current  state  of  the  art  of  ASM's  with 
regard  to  methods  of  control,  hardware  limitations,  terminal  accuracy,  and 
effects  of  atmospheric  disturbances.  The  results  of  the  survey  are  presented 
in  Appendix  I.  Applications  of  this  information  to  terminal  guidance  and  auto¬ 
pilot  studies  are  discussed  in  this  section. 

(2)  Optimal  control  theory  was  applied  to  the  development  of  an  optimal  control 
law  for  the  near-impact  phase  that  can  be  implemented  on-board  an  ASM.  The 
optimal  guidance  law  development  is  presented  in  Section  V.  Steepest  descent 
trajectory  optimization  techniques  were  applied  to  the  development  of  nominal 
trajectories.  The  trajectory  optimization  results  are  presented  in  the  classi¬ 
fied  supplement  to  this  report. 

(3)  The  effects  of  various  control  laws,  control  component  characteristics,  and 
natural  disturbances  on  terminal  accuracy  were  determined  and  the  results 
of  these  studies  were  used  to  develop  control  system  design  guides  for  the 
midcourse  and  near-impact  phases  of  the  ASM  mission.  These  guidelines  are 
presented  in  Section  IV . 

(4)  A  normal  acceleration  autopilot  was  designed  to  meet  the  mission  requirements 
and  advanced  bistable  controller  techniques  were  applied  to  obtain  a  quasi - 
adaptive  autopilot  which  required  no  gain  changes  throughout  the  midcourse 
and  terminal  ASM  flight.  Autopilot  design  study  results  and  discussions  of 
other  adaptive  approaches  are  presented  in  Section  VI. 

Following  are  descriptions  of  the  baseline  missile,  control  laws,  homing  sensors, 
and  atmospheric  disturbances  used  in  the  study. 

The  baseline  ASM  configuration  used  in  this  study  is  based  on  Boeing  air-to-surface 
missile  preliminary  design  studies.  As  a  result,  the  configuration  incorporates 
typical  vehicle  design  features  and  subsystem  integration  characteristics.  The 
assumptions  that  were  made  in  the  analytical  description  of  the  missile,  its  var¬ 
ious  subsystems,  and  environment  are  outlined  in  this  section.  Descriptions  are 
provided  for  the  following  items: 


(1)  Missile  configuration; 

(2)  Nominal  trajectories; 

(3)  Homing  guidance  laws; 

(4)  Inertial  guidance  sensor; 

(5)  Homing  sensor; 

(6)  Autopilot; 

(7)  Atmospheric  disturbances; 

(8)  Flight  simulation  and  miss  distance  comparison. 


Classified  details  of  the  missile  configuration  and  trajectory  optimization  are 
presented  in  the  supplement  to  this  report.  Separate  design  studies  were  made 
of  an  optimal  guidance  law  and  a  normal  acceleration  autopilot.  These  are  pre¬ 
sented  in  more  detail  in  later  sections. 

1.  Missile  Configuration.  A  missile  configuration  that  was  defined  in  a  Boeing 
ASM  program  provides  realistic  characteristics  for  use  in  this  study.  The  con¬ 
figuration  is  capable  of  a  variety  of  missions  through  the  use  of  a  two-pulse  solid- 
propellant  motor.  The  operational  flexibility  permits  the  use  of  either  optical  or 
radar -type  homing  sensors. 

The  geometric,  inertial,  aerodynamic  (including  roll-yaw  coupling  effects),  and 
structural  bending  characteristics  of  this  missile  are  contained  in  the  supplement 
to  this  report.  The  aerodynamic  characteristics  used  in  the  study  were  based  on 
wind  tunnel  data  for  this  configuration.  The  data  required  for  the  autopilot  studies 
were  reduced  to  transfer  function  form  and  are  presented  in  Section  VI. 

In  general,  the  aerodynamic  stability  derivatives  are  nonlinear  functions  of  Mach 
number  and  angle  'if  attack.  Significant  nonlinear  characteristics  were  included 
in  the  study.  However,  certain  simplifying  assumptions  were  made  in  cases  where 
they  provided  reasonable  representation  of  the  wind  tunnel  data.  The  body  normal 
force  curve  slope  is  strongly  dependent  on  angle  of  attack  and  Mach  number  ;  the 
body  normal  force  coefficient  is  closely  approximated  by: 

CNBody  -  CNaBody  **  +  CNa3Body  0  (1) 
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This  equation  was  used  with  CN^Body  anci  CN&3Body  Senerated  as  functions  of 
Mach  number.  The  wind  tunnel  data  indicated  that  these  two  contributions  to  total 
body  normal  force  may  be  assumed  to  act  at  two  distinct  centers  of  pressure,  with 
the  location  of  the  cubic  term  fixed  and  the  location  of  the  linear  term  a  function 
of  Mach  number.  This  assumption  was  used  in  the  analog  simulation. 

The  body  axial  force  coefficient  and  the  fin  normal  force  curve  slope  are  very  weak 
functions  of  angle  of  attack,  for  angles  of  attack  less  than  20  degrees.  They  were, 
therefore,  assumed  to  be  functions  of  Mach  number  only.  The  aerodynamic  center 
of  the  control  fins  varies  with  angle  of  attack  and  Mach  number,  but  since  the  vari¬ 
ation  represents  a  very  small  percentage  of  the  total  lever  arm,  a  constant  control 
arm  was  assumed. 

2.  Trajectories.  Three  reference  trajectories,  consisting  of  a  "low -altitude 
skip,"  a  "short-range  semiballistic,"  and  a  "long-range  semiballistic"  mission 
of  the  form  shown  in  Figure  1  were  initially  selected  for  study.  This  investigation 
was  primarily  concerned  with  the  homing  (or  terminal)  phase  of  the  trajectories. 
Tnis  phase  was  considered  for  acquisition  slant  ranges  of  15,000  to  30,000  feet 
to  the  target. 

The  low -altitude  skip  trajectory  utilizes  a  two-pulse  mode  solid -propellant  motor. 
The  short-range  semiballistic  trajectory  and  the  low -altitude  skip  trajectory  both 
achieved  satisfactory  results  without  violating  reasonable  angle  of  attack  or  mini¬ 
mum  dynamic  pressure  constraints.  These  trajectories  were  obtained  using  con¬ 
ventional  open -loop  performance  evaluation  trajectory  programs.  Study  of  the 
short-range  semiballistic  case  was  discontinued  after  preliminary  evaluation 
studies  indicated  insignificant  differences  in  miss  distances. 

For  the  long-range  trajectory  it  was  necessary  to  employ  trajectory  optimization 
techniques  to  provide  a  flight  path  with  sufficiently  large  dynamic  pressure  at  apogee 
so  that  aerodynamic  controls  could  function  to  keep  the  vehicle  from  tumbling.  The 
optimization  was  performed  with  a  fixed  range  requirement  using  constraints  on 
the  minimum  dynamic  pressure  at  apogee  of  10,  20,  35,  and  52  psf.  In  addition, 
appropriate  constraints  were  imposed  on  angle  of  attack  and  angle  of  attack  rate 
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RANGE 


Figure  1:  REFERENCE  TRAJECTORIES 


during  the  pull-up  maneuver  after  launch  from  the  launch  vehicle  and  the  velocity 
at  impact  was  maximized.  The  maximization  of  impact  velocity  is  an  important 
factor  for  missions  involving  penetration  to  defended  targets.  Data  and  results 
regarding  this  optimization  process  are  shown  in  Figures  1,2,  and  3  of  the 
supplement. 

3.  Homing  Guidance  Laws.  Conventional  proportional  and  pursuit  homing  guid¬ 
ance  laws  were  compared  with  an  "optimal"  homing  guidance  concept.  For  pro¬ 
portional  and  pursuit  guidance,  homing  guidance  operated  in  the  terminal  phase 
until  the  ASM  reached  the  homing  sensor  blind  range  from  the  target.  At  that 
point  options  were  provided  either  to  zero  the  guidance  commands  or  to  hold  the 
last  command  value.  Optimal  guidance  used  position  information  from  IMU 
measurements  to  guide  during  the  blind  zone. 

The  guidance  laws  are  based  upon  use  of  a  normal  acceleration  autopilot.  For 
proportional  guidance  the  normal  acceleration  command,  *z'c,  is  defined  as: 

zc  -  KX  V  7?  (2) 

where  Kj  is  the  guidance  gain,  V  is  the  missile  velocity,  and  r\  is  the  angular 
rate  of  the  target  line  of  sight.  Another  form  of  the  pro; ortional  guidance  law, 
which  is  useful  for  comparisons  with  pursuit  guidance,  can  be  derived  through 
small  angle  approximations.  This  is: 

zc  =  ^1  r  (ft*.  -  Of)  (3) 

^sl 

where  R8j  is  the  slant  range  to  the  target.  The  pursuit  guidance  law  is: 

zc  =  K2  <6t  -  a)  (4) 

The  basic  difference  between  the  two  laws  is  the  V^/Rsj  term ,  and  the  effective 
pursuit  guidance  gain  remains  constant  throughout  the  terminal  phase;  however, 
for  the  proportional  guidance  law,  gain  increases  as  the  missile  approaches  the 
target  (see  Figure  2  for  definitions). 

Proportional  guidance  can  readily  be  implemented  with  a  gyro -stabilized  gimbaled 
tracker  by  using  the  output  signals  from  the  rate  gyros  used  in  the  stabilization 
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Pursuit  guidance: 


2c  =  K2(  tj  +y) 


also 


2c  =  K2<  0T-a) 
Proportional  guidance: 


(5) 

(6) 


*c  =  KlV* 

=  K,V 


v(ij  +  -y) 


4  =K  v 

c  *1  —  ( V  +  y) 


{ H  +  Y)  =  small  angle  (7) 


(8) 


also 


Ksl 


(9) 


Figure  2:  GUIDANCE  LAW  AND  ANGLE  DEFINITIONS 
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Pursuit  guidance: 

Zc  =  k2(t)  +  y) 

(5) 

also 

W  *T“a) 

* 

(6) 

Proportional  guidance: 

r  i 

=  K.V  , 

*,l 

L  J 

(t)  +  y)  =  small  angle  (7) 

2  -K  v2 
c  i  T™  (»i  +  y) 

si 

(8) 

Figure  2:  GUIDANCE  LAW  AND  ANGLE  DEFINITIONS 


and  tracking  loops  of  the  tracker.  Pursuit  guidance  is  more  difficult  to  implement 
because  it  requires  a  measurement  of  the  target  line  of  sight  with  respect  to  the 
missile  velocity  vector.  This  is  not  a  directly  measurable  quantity  since  angle  of 
attack  must  be  subtracted  from  the  line -of -sight  angles  that  the  sensor  measures 
from  the  body  axes. 

The  optimal  guidance  law  investigated  minimizes  the  sum  of  the  mean  square  miss 
distance  plus  the  integral  of  a  quadratic  form  of  the  autopilot  control  variables.  A 
weighting  factor  on  the  integral  term  allows  the  relative  weighting  of  the  two  terms 
to  be  adjusted.  This  provides  a  means  of  minimizing  the  miss  distance  with  real¬ 
istic  restraints  on  control  action.  The  optimal  guidance  law  is  based  on  the  use 
of  a  Kalman  filter  to  obtain  the  best  estimates  of  the  ASM  state  variable  errors 
from  imperfect  sensor  information.  Sensor  information  in  the  form  of  accelera¬ 
tions  measured  with  an  on-board  inertial  platform  and  line-of -sight  information 
from  the  homing  guidance  sensor  are  used  in  the  computation.  Details  of  the  de¬ 
velopment  of  the  optimal  control  law  are  presented  in  Section  V. 

4.  Inertial  Guidance  Sensors.  A  conventional  gimbaled  type  of  inertial  measure¬ 
ment  unit  (IMU)  was  chosen  for  the  baseline  inertial  guidance  sensor.  This  choice 
was  made  in  part  because  this  was  the  only  type  of  instrument  for  which  good 
accuracy  and  cost  information  were  available.  Equipment  variations  for  this  type 
of  inertial  measurement  unit,  relating  cost  to  sensor  accuracy,  were  exercised 
with  an  upper  limit  on  cost  of  the  inertial  platform  and  associated  electronics  of 
$40,000  (excluding  the  airborne  computer).  Three  gradations  of  equipment  accu¬ 
racy  versus  cost  were  investigated. 

The  la  errors  in  three  inertial  position  and  velocity  components  were  evaluated 
as  measured  by  the  inertial  guidance  sensor  along  each  of  the  three  reference 
trajectories.  These  sets  of  accuracy  data  represent  the  errors  generated  during 
midcourse  flight  aDd  provided  the  initial  errors  for  the  terminal  homing  phase  of 
flight.  A  summary  of  this  data  in  the  form  of  la  root  sum  square  position  errors 
is  presented  in  Table  I.  The  data  Indicates  that  little  accuracy  is  gained  by  going 
beyond  the  $25,000  IMU.  Therefore,  this  IMU  was  selected  as  baseline.  Position 
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fixing  and  alignment  prior  to  launch  were  based  on  low-altitude  launch  from  the 
launch  vehicle.  Additional  details  of  the  inertial  guidance  sensor  survey  and 
error  analysis  are  presented  in  Appendix  I. 

This  error  data  was  used  in  the  ASM  miss  distance  studies  in  two  ways.  The 
IMU  la  errors  were  used  directly  in  the  initial  covariance  matrix  of  position  and 
velocity  errors  for  a  digital  linear  error  analysis  program.  In  a  nonlinear  ana¬ 
log  study,  the  IMU  Iff  positional  errors  for  each  component  were  used  as  finite 
x,  y,  and  z  component  errors  to  give  a  worst-case  deviation  from  the  nominal 
initial  conditions.  The  same  set  of  component  errors  were  used  for  all  the  ana¬ 
log  studies.  For  the  baseline  IMU,  inertially  referenced  angle  measurement 
errors  were  small  (less  than  1  degree).  However,  the  position  errors  result  in 
an  initial  velocity  vector  pointing  error  with  respect  to  the  target.  This  initial 
angular  error  was  approximately  3.0  to  4.0  degrees  for  the  cases  studied. 

Table  I:  INERTIAL  MEASUREMENT  UNJT  la  RSS  POSITION  ERRORS 
AT  TARGET  ACQUISITION 


Slant  Range 

RSS  Errors  (—feet) 

(feet) 

$13,000* 

$25,000* 

$40,000* 

Short-Range 

15 , 000 

1670 

1450 

1450 

Semiballistic 

Trajectory 

30,000 

1650 

1450 

1450 

60,000 

1620 

1450 

1450 

Low -Altitude 

15,000 

7910 

2450 

2330 

Skip  Trajectory 

30,000 

6680 

2250 

2130 

60,000 

4740 

2000 

1770 

Long-Range 

15 , 000 

5210 

2850 

2730 

Semiballistic 

Trajectory 

30,000 

5030 

2800 

2660 

60,000 

4690 

2650 

2540 

*  Inertial  Measurement  Unit  Cost. 

5.  Homing  Sensor.  Surveys  of  homing  guidance  sensors  were  conducted  to  ob¬ 
tain  operational  characteristics  of  various  TV,  infrared,  and  radar -type  sensors 
to  determine  which  of  these  sensors  would  be  suitable  for  the  study  missions.  As 
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a  result  of  these  surveys  (which  are  presented  in  Appendix  I),  a  TV -type  homing 
sensor  was  selected  for  the  study.  Because  of  a  lack  of  information  on  the  char¬ 
acteristics  of  wide  field  of  view  trackers,  and  because  of  the  conceptual  equivalence 
of  body -fixed  and  gimbaled  sensors  as  applied  to  the  optimal  guidance  task  (dis¬ 
cussed  in  Section  V) ,  the  studies  were  based  on  the  characteristics  of  narrow  field 
of  view,  gimbaled  TV  trackers  only.  A  simplified  single  axis  block  diagram  repre¬ 
sentative  of  this  type  of  tracker  is  shown  in  Appendix  I. 

Bias  errors  and  output  noise  were  considered  representative  of  the  types  of  errors 
to  be  encountered  in  TV  trackers.  The  errors  were  investigated  in  parametric 
form ,  because  detailed  performance  data  were  limited.  Bias  errors  typically 
would  result  from  boresight  alignment  errors  between  the  TV  vidicon  tube  and  its 
gimbal  structure  and  misalignment  between  the  tracker  base  and  the  missile  ref¬ 
erence  axes.  Major  sources  of  noise  are  probably  the  target  background,  vidicon 
tube  and  its  signal  processing  circuitry,  and  the  rate  sensors  used  to  stabilize  the 
tracker  gimbals.  The  noise  was  considered  to  be  "white  noise. "  In  the  analog 
studies  the  noise  was  obtained  from  a  low-frequency  noise  generator  with  a  30-cps 
cutoff  frequency.  The  noise  generator  was  scaled  to  a  nominal  tracker  noise  level 
of  1  milliradJan  (rms).  In  the  linear  analysis,  the  nominal  noise  was  considered 
to  be  "white  noiBe"  with  a  standard  deviation  of  1. 0  milliradian.  As  long  as  the 
response  of  the  ASM  system  is  significantly  below  the  30-cps  cutoff  frequency  of 
the  noise  generator,  the  two  techniques  are  comparable.  Nominal  values  of  1.74 
milllradians  (0.1  degree)  bias  error  and  1. 0  milliradian  (rms)  of  output  noise 
were  assumed.  (An  experimental  program  is  described  in  Appendix  I  in  which  a 
centroid  tracker  was  synthesized.  This  experiment  yielded  bias  error  values  of 
less  than  1  milliradian.) 

Nominal  tracker  gimbal  position  dynamics  were  considered  to  be  a  5-cps,  0.5- 
damped -quadratic  response  to  line -of -sight  errors  measured  by  the  optical  sensor. 
Because  the  program  was  concerned  primarily  with  the  terminal  homing  phase  of 
flight,  Investigation  of  the  dynamics  and  peculiarities  of  target  acquisition  were 
not  covered. 


6.  Autopilot.  A  normal  acceleration  autopilot  was  selected  for  control  of  the 
ASM.  This  choice  was  based,  primarily,  on  previous  Boeing  homing  dynamics 
studies  that  indicated  the  desirability  of  a  normal  acceleration  commanded  auto¬ 
pilot  because  of  the  simplification  (over  attitude  control)  in  guidance  laws  it  pro¬ 
vides.  In  addition,  the  normal  acceleration  autopilot  provides  a  parameter  for 
limiting  command  maneuvers  to  stay  within  vehicle  structural  limits. 

Autopilot  stability  and  time  response  studies  (discussed  in  Section  VI)  resulted  in 
a  normal  acceleration  plus  angular  rate  feer"  ack  system  with  forward  loop  com  - 
pensation  to  reduce  steady-state  errors.  Analog  terminal  phase  studies  were 
conducted  to  determine  the  effects  of  autopilot  characteristics  on  miss  distance. 
These  studies  used  a  detailed  representation  of  the  autopilot  as  shown  in  Figure  3. 
An  idealized  version  was  used  in  the  digital  evaluation  programs,  in  which  instan¬ 
taneous  response  to  commanded  accelerations  was  assumed. 

7.  Atmospheric  Disturbances.  Two  types  of  atmospheric  disturbances  were  con¬ 
sidered  in  the  study.  The  wind  shear  profile  as  shown  in  Figure  4  was  obtained 
by  adding  a  spike  to  Figure  9  of  Reference  1  and  is  considered  representative  of 
wind  shear  within  a  1%  probability.  Atmospheric  turbulence  with  a  power  spectral 
density  described  on  Page  63  of  Reference  1  was  added  to  the  analog  simulation. 

8.  Flight  Simulation  and  Miss  Distance  Comparison.  Two  approaches  were  used 
in  evaluating  the  effects  of  various  system  characteristics  on  miss  distance  error 
of  the  ASM.  The  effects  of  proportional  and  pursuit  guidance  gain,  variations  in 
sensor  and  autopilot  characteristics,  and  the  effects  of  atmospheric  disturbances 
were  determined  primarily  with  an  analog  simulation  of  the  homing  phase.  The 
optimal  guidance  law  is  not  amenable  to  simulation  on  the  analog  computer;  there¬ 
fore,  a  digital  computer  program  was  developed  to  permit  comparative  evaluation 
of  miss  distance  for  all  three  guidance  laws  (i.e. ,  pursuit,  proportional,  and 
optimal  guidance)  using  a  simplified  missile  representation  that  assumed  an  ideal¬ 
ized  autopilot.  The  use  of  the  two  analysis  techniques  and  the  detail  considered 
in  each  is  summarized  in  Table  II. 

The  analog  simulation  was  previously  developed  specifically  for  homing  guidance 
studies  (Reference  2).  Two  automatic  scale  changes  in  the  three  roissile-to-target 
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Figure  3:  GUIDANCE  AND  CONTROL  SYSTEM  -  ANALOG 


Table  II:  MISS  DISTANCE  EVALUATION  TECHNIQUES 


Purpose 

Equations 

6DOF 

3DOF 

Analog 

Digital 

Guidance  Law  Com¬ 
parison 

Linearized 

X 

X 

Optimal  Guidance 

Linearized 

X 

X 

Tracker  Errors 

Linearized 

Nonlinear 

X 

X 

X 

X 

Guidance  Gain,  Auto¬ 
pilot 

Nonlinear 

X 

X 

Atmospheric  Turbu¬ 
lence 

Nonlinear 

X 

X 

displacement  components  were  made  to  provide  satisfactory  definition  of  impact 
accuracy.  Six  degrees  of  freedom  were  simulated.  The  equations  of  motion  for 
pitch  and  yaw  rotations  and  three  translations  in  a  target  centered  cartesian  coor¬ 
dinate  system  were  solved.  The  missile  roll  axis  control  was  considered  to  be 
ideal  and  maintained  one  pair  of  the  cruciform  fins  in  the  vertical  plane .  Inertial 
cross -coupling  effects  were  included.  Nonlinear  aerodynamics,  proportional  or 
pursuit  guidance,  and  the  autopilot  as  shown  in  Figure  3  were  included  in  the  simu¬ 
lation.  The  terminal  guidance  TV  tracker  noise  and  bias,  and  tracker  dynamics 
as  described  in  5 .  Homing  Sensor  were  also  included.  Provisions  were  made 
for  atmospheric  disturbance  simulation  in  the  form  of  horizontal  wind  shear  pro¬ 
files  and  turbulence  that  acted  normal  to  the  missile  body  axis. 

The  homing  portion  of  the  missile  flight  was  Initiated  in  the  analog  simulation  at 
a  slant  range  of  15,000  feet  (from  approximately  1000  to  15,000  feet  correspond¬ 
ing  altitude)  for  the  three  nominal  trajectories  —  long-range  semiballistic ,  short- 
range  semiballistic ,  and  low -altitude  skip.  Simulation  during  the  region  of  TV 
tracker  blind  range  provided  the  normal  acceleration  commands  from  the  guidance 
system  to  be  either  set  to  zero  or  held  at  the  last  commanded  value  at  a  pre¬ 
selected  range  from  the  target.  The  simulated  autopilot  was  the  normal  accel¬ 
eration  plus  pitch  rate  feedback  system  with  forward  loop  compensation.  To 
obtain  adequate  and  uniform  response  and  stability,  the  two  autopilot  gains  were 


changed  continually  throughout  the  flight  proportional  to  the  Inverse  of  dynamic 
pressure  and  the  square  root  of  its  inverse.  Details  of  the  autopilot  character¬ 
istics  are  presented  in  Section  VI. 

Comparison  data  on  miss  distance  for  all  three  guidance  concepts  (proportional, 
pursuit,  and  optimal)  were  obtained  with  a  digital  program  for  computing  the  co- 
variance  matrix  of  position  and  velocity  errors  from  the  nominal  trajectories .  A 
simplified  description  of  the  ASM  was  used  in  this  program.  The  missile  was 
considered  to  be  a  point  mass  with  lift  and  drag  force  coefficients  as  nonlinear 
functions  of  angle  of  attack  and  Mach  number.  An  ideal  autopilot  was  assumed 
(i.e. ,  instantaneous  response  to  commanded  accelerations).  TV  tracker  dynamics 
were  included;  sensor  errors  were  assumed  to  consist  of  random  bias  plus  white 
noise. 

In  this  approach  the  equations  of  motion  of  the  ASM  were  linearized  about  a  nomi¬ 
nal  trajectory.  The  initial  conditions  for  the  covariance  matrix  were  the  position 
and  velocity  errors  at  target  acquisition  that  were  based  on  the  expected  disper¬ 
sions  developed  during  midcourse  from  inertial  instrument  errors.  The  covari¬ 
ance  matrix  of  the  state  of  the  linearized  missile  equations,  as  a  function  of  time, 
was  then  obtained  from  a  numerical  solution  of  the  linearized  differential  equations. 
The  integration  was  performed  on  a  Sperry  Rand  1107  computer.  The  miss  dis¬ 
tance  was  obtained  as  the  value  of  the  covariance  matrix  at  the  time  (he  nominal 
trajectory  hits  the  target.  This  evaluation  technique  is  described  in  detail  in 
Appendix  n. 

Miss  distances  are  defined  in  a  plane  passing  through  the  target  and  normal  to  the 
ASM  trajectory  as  shown  in  Figure  5.  These  errors  are  more  direct  than  errors 
measured  in  the  ground  plane,  because  they  are  independent  of  the  ASM  target 
approach  angle.  In  the  following  discussions  of  results,  errors  obtained  from 
the  linear  analysis  are  presented  as  Iff  (or  standard)  deviations  about  the  target 
as  shown  in  Figure  5 . 

Miss  distances  obtained  from  the  analog  studies  are  presented  in  two  forms:  an 
average  error  and  the  rms  deviation  from  the  average.  As  discussed  previously, 
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Figure  5:  MISS  DISTANCE  DEFINITION 


fee  initial  positional  errors  were  introduced  as  a  finite  set  of  errors  from  the 
nominal  trajectory.  Therefore,  the  analog  miss  distances  obtained  from  repeated 
simulation  runs  were  not  centered  about  the  target  like  those  obtained  from  the 
linear  analysis.  This  deviation  of  the  center  of  the  group  of  miss  distances  from 
the  target  is  presented  as  an  average  (or  arithmetic  mean)  error.  With  the  set 
of  initial  position  errors  used  in  the  analog  study,  this  average  error  was  always 
positive.  Obviously,  a  different  set  of  initial  conditions  could  have  been  selected 
that  would  have  resulted  in  negative  average  errors.  The  analog  rms  miss  dis¬ 
tances  that  are  presented  are  the  calculated  rms  deviations  from  the  average 
error  and  represent  the  effect  of  stochastic  disturbances  to  the  ASM.  From  15 
to  20  simulation  runs  were  made  for  each  analog  rms  data  point. 


SECTION  IV 


CONTROL  SYSTEM  DESIGN  GUIDES 

The  control  system  design  guides  are  presented  in  three  categories: 

(1)  Selection  of  Guidance  Concepts 

These  data  present  a  comparison  of  achievable  miss  distances  for  three 
guidance  concepts  as  a  function  of  nominal  trajectories  and  atmospheric 
disturbances.  The  effect  of  gains  for  each  guidance  concept  is  also  dis¬ 
cussed. 

(2)  Selection  of  Inertial  and  Homing  Sensors 

Effects  of  dispersion  at  target  acquisition,  homing  sensor  bias  errors  and 
output  noise,  and  homing  sensor  blind  range  are  discussed. 

(3)  Selection  of  Autopilot 

Effects  of  autopilot  short-period  control  dynamics  on  miss  distance  are 
presented. 

Discussion  of  the  details  of  conditions  investigated  and  the  results  obtained  arr 
presented  in  the  following  paragraphs. 

1.  Selection  of  Guidance  Concept.  Three  homing  guidance  concept!  (propor¬ 
tional,  pursuit,  and  optimal  guidance)  were  evaluated  for  impact  accuracy. 
Guidance  gains  were  optimized  for  proportional  and  pursuit  guidance;  the  effect 
of  the  weighting  factor  on  the  integral  term  of  the  optimal  guidance  performance 
index  was  explored.  (This  performance  index  minimizes  the  sum  of  mean  square 
miss  distance  plus  the  integral  of  the  quadratic  form  of  the  control  variables. ) 

Evaluations  were  made  for  two  nominal  trajectories.  Two  acquisition  slant 
ranges  were  evaluated  for  one  of  the  trajectories.  A  range  of  atmospheric  tur¬ 
bulence  from  0  to  12  fps  rms  was  simulated  in  the  analog  computer  evalus.tr 
of  proportional  guidance. 

a.  Comparison  of  Guidance  Concepts.  Table  III  compares  the  obtainable  miss 
distances  for  proportional,  pursuit,  and  optimal  guidance  with  the  nominal  ASM 
system  errors.  Guidance  gains  were  set  at  or  near  the  best  values  for  each 
guidance  technique.  Error  data  is  shown  only  from  the  linear  analysis  program 
and,  therefore,  represents  the  Id  (or  standard)  deviations  about  the  target  (as 
discussed  in  Section  III).  Average  error  data  (as  shown  in  Figures  6  and  7)  from 
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Figure  6:  EFFECT  OF  GUIDANCE  GAIN 

FOR  PROPORTIONAL  GUIDANCE 
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Figure  7: 


•  Nominal  Dispersion  at  Target  Acquisition 

•  Slant  Range  at  Acquisition  =  15,000  Ft 

•  Tracker  Noise  =  1  .0  Milliradian  (rms) 

•  Tracker  Bias  Error  =  1 .74  Milliradians 
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the  6-degree-of- freedom  analog  simulation  were  reasonably  close  to  the  Iff  errors 
shown  in  Table  in.  This  indicates  that  for  nominal  conditions,  the  vehicle  re¬ 
sponse  is  not  degrading  the  ASM  performance  and  the  simplified  linear  analysis 
is  adequate. 

Each  of  the  guidance  concepts  is  capable  of  achieving  small  miss  distances;  the 
spread  between  minimum  and  maximum  Iff  errors  is  less  than  2  feet.  The  results 
show  optimal  guidance  to  be  better  than  either  proportional  or  pursuit,  and  pro¬ 
portional  to  be  slightly  better  than  pursuit. 


Table  III:  GUIDANCE  LAW  COMPARISON  — 

MISS  DISTANCE  MEASURED  NORMAL  TO  TRAJECTORY 


Iff  Miss  Distance  (feet) 

Trajectory 

Proportional 

Guidance 

Pursuit 

Guidance 

Optimal 

Guidance 

Long-Range  Semiballistic 

Acquisition  at: 

30, 000- foot  Slant  Range 

2.7 

2.9 

1.6 

15,000-foot  Slant  Range 

3.1 

3.3 

2.2 

Low- Altitude  Skip 

Acquisition  at: 

15, 000- foot  Slant  Range 

— 

1.74 

1.44 

Conditions: 

Baseline  IMU  Iff  dispersion  at  acquisition 
Tracker  noise  —  I  milliradian 
Tracker  bias  —  I.  74  milliradians 
Tracker  blind  range  —  500  feet 

It  can  be  seen  from  Table  HI  that  the  differences  in  miss  distance,  as  measured 
perpendicular  to  the  trajectory,  between  trajectories,  and  for  different  slant 
ranges  at  target  acquisition  are  very  small.  A  third  nominal  trajectory  (short- 
range  semlballistlc)  was  dropped  from  evaluation  when  initial  studies  indicated 


no  significant  differences  between  trajectories.  When  the  measured  miss  dis¬ 
tances  are  resolved  into  the  ground  plane  through  the  trajectory  approach  angle, 
however,  the  deviation  oi  impact  on  the  ground  for  the  low-altitude  skip  trajectory 
becomes  approximately  ten  times  that  of  the  long-range  semiballistic  trajectory, 
due  to  its  very  shallow  target  approach  angle.  For  the  range  of  trajectories 
investigated,  the  deviation  normal  to  the  trajectory  can  be  considered  independent 
of  the  trajectory.  This  is  significant  because  it  indicates  that  the  ASM  perform¬ 
ance  is  not  seriously  affected  by  the  terminal  phase  flight  conditions.  On  the  two 
trajectories  investigated,  terminal  velocities  are  different  by  approximately  a 
factor  of  two. 

b.  Effect  of  Guidance  Gain.  Figures  6  and  7  present  the  effects  of  guidance 
gain  on  proportional  and  pursuit  guidance  accuracy  with  the  nominal  ASM  system 
errors.  As  discussed  in  the  presentation  of  tracker  error  effects,  tracker  bias 
errors  were  not  included  in  the  proportional  guidance  simulations.  These  data 
were  obtained  from  the  6-degree-of-freedom  analog  simulation  described  in 
Section  in  that  included  complete  autopilot  and  guidance  systems.  Only  the 
average  error  (as  described  in  Section  IE)  is  shown  for  clarity.  The  rms  devi¬ 
ation  from  the  average  is  less  than  0.  5  foot  for  the  nominal  conditions  and  best 
gains.  Note  that  the  errors  for  the  best  gains  are  comparable  to  those  shown 
in  Table  III  for  the  ljflpar  analysis. 

The  increase  in  miss  distance  at  low  guidance  gains  for  both  concepts  is  the 
result  of  inadequate  gain  to  remove  the  effect  of  initial  position  errors  (and  hence 
velocity  vector  pointing  errors).  At  low  gains  the  ASM  is  essentially  operating 
in  an  open-looo  or  unguided  condition.  As  discussed  in  Section  III,  the  velocity 
vector  pointing  error,  for  the  nominal  initial  position  error  of  approximately 
1500  feet,  is  about  3  to  4  degrees. 

The  effect  of  proportional  guidance  gain,  K^,  on  miss  distance.is  shown  in 
Figure  6  and  indicates  that  best  accuracy  is  obtained  with  a  guidance  gain  between 
3  and  4.  To  explain  the  degradation  of  performance  as  the  guidance  gain  is  in¬ 
creased  beyond  4,  remember  that  the  effective  proportional  guidance  gain  is 
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V  /Rgi  (shown  in  Section  HI).  Thus,  as  the  ASM  approaches  the  target  the 
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effective  gain  increases,  and  at  some  range,  Rgj,  stability  problems  will  occur. 
In  Figure  6,  the  range  is  shown  at  which  the  combined  autopilot/guidance  gain 
margin  is  reduced  to1 6  db  for  the  ASM  configuration  studied.  In  obtaining  the 
data  for  Figure  6,  guidance  was  terminated  (z"c  =  0)  at  the  indicated  ranges.  The 
long  blind  ranges  allow  insufficient  time  to  remove  the  effect  of  initial  position 
errors  and  large  miss  distances  result.  A  similar  increase  in  miss  distance 
would  be  noted  at  guidance  gains  beyond  4  or  5  because  of  decreased  stability  if 
guidance  was  not  terminated. 

Stability  was  not  a  problem  for  proportional  guidance  gains  less  than  4  and  slant 
ranges  down  to  500  feet,  and  all  remaining  data  runs  were  conducted  without  the 
6-db-gain  margin  restriction.  As  will  be  shown  later  in  the  discussion  of  homing 
sensor  blind  range  effects,  long  blind  ranges  are  not  desirable  when  atmospheric 
turbulence  is  present. 

The  effect  of  pursuit  guidance  gain  on  missile  accuracy  is  shown  in  Figure  7.  It 
may  be  seen  that  the  low-altitude  skip  trajectory  requires  less  guidance  gain  than 
does  the  long-range  trajectory.  The  reason  for  this  difference  is  that  the  veloci¬ 
ties  for  the  two  trajectories  differ.  The  effective  gain  for  proportional  guidance 
is  KjV2/R8j.  The  best  Kj  was  found  to  be  independent  of  the  missile  velocity 
for  the  trajectories  investigated.  Because  the  effective  guidance  gain  is  a  con¬ 
stant  for  pursuit  guidance,  it  follows  that  for  these  trajectories  the  gain  for  best 
accuracy  will  be  proportional  to  the  square  of  the  velocity.  The  ratio  of  the 
squares  of  the  initial  velocities  of  the  two  trajectories  is  6  to  1,  which  co¬ 
incides  very  closely  to  the  ratio  of  the  minimal  gains  tor  acceptable  accuracy 
as  shown  in  Figure  7. 

Because  pursuit  guidance  gain  is  independent  of  range,  missile  stability  does 
not  decrease  as  the  missile  approaches  the  target.  Therefore,  pursuit  guidance 
cutoff  range  can  be  selected  independent  of  guidance  gain.  For  the  pursuit  guid¬ 
ance  runs  represented  by  Figure  7,  the  guidance  cutoff  slant  range  (blind  range) 
was  set  at  500  feet. 
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The  termination  of  the  low-altitude  skip  trajectory  curve  at  K  =  2500,  as  com- 
pared  to  12, 500  for  the  long-range  trajectory,  is  the  result  of  an  instability 
caused  by  excessive  acceleration  commands  at  acquisition.  The  problem  can 
probably  be  overcome  with  acceleration  command  and/or  acceleration  error 
limiters.  Howeve:',  this  extension  was  not  attempted  because  the  terminal  gains 
in  each  case  represent  a  gain  approximately  twice  the  minimum  acceptable  value. 
Acceleration  command  limiting  was  not  required  at  other  times  in  the  trajectory. 

The  effect  of  the  weighting  factor  on  the  control  integral  term  of  the  optimal 
control  performance  index  was  not  investigated  in  detail.  (This  term,  which  is 
discussed  in  Section  V,  essentially  limits  the  ASM  maneuver  capability. )  How¬ 
ever,  it  is  believed  that  considerable  latitude  is  possible  in  the  selection  of  this 
factor.  Emphasis  in  this  investigation  was  on  obtaining  a  weighting  factor  that 
would  provide  satisfactory  miss  distances.  As  was  shown  in  Table  in,  this 
objective  was  accomplished.  The  weighting  factor  that  yielded  these  miss  dis¬ 
tances  was  two  orders  of  magnitude  below  the  value  of  the  first  factor  tried.  The 
larger  factor  had  given  completely  unsatisfactory  errors.  Detailed  simulation 
studies  of  the  optimal  control  concept  would  be  required  to  compare  the  optimal 
guidance  maneuver  requirements  with  those  of  a  conventional  guidance  concept. 
Studies  of  this  type  were  beyond  the  scope  of  this  investigation. 

c .  Effect  of  Atmospheric  Disturbances.  Figure  8  shows  the  effect  of  turbulence 
on  impact  error  for  the  long-range  trajectory,  using  proportional  guidance  with 
a  blind  range  of  500  feet.  (The  6-db-gain  margin  restriction  on  overall  gain  was 
ignored. )  The  gust  characteristics  obtained  from  Reference  1,  and  discussed  in 
Section  III,  were  applied  in  the  vertical  plane.  It  may  be  seen  that  miss  distance 
insensitive  to  rms  gust  velocity  for  the  blind  range  considered.  Although  the 
probability  of  occurrence  is  very  slight,  digital  runs  were  made  with  approxi¬ 
mately  20-fps  rms  turbulence  with  no  effect  on  miss  distance.  This  is  due  pri¬ 
marily  to  four  factors:  First,  the  high  missile  velocity  combined  with  relatively 
low  gust  velocities  results  in  small  angle- of- attack  perturbations;  Second,  the 
random  nature  of  gusts  tends  to  make  the  average  effect  on  the  trajectory  small; 
Third,  with  a  500- foot  blind  range,  any  transient  due  to  gusts  existing  at  the 
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time  of  guidance  cutoff  has  little  time  (approximately  0.  25  second)  to  affect  the 
impact  point;  Fourth,  the  guidance  system  tends  to  compensate  for  gusts. 

However,  turbulence  does  contribute  to  miss  distance  when  blind  rmge  is 
increased,  as  is  shown  in  the  discussion  of  blind  range  (Paragraph  2-d  of  this 
section). 

The  wind  shear  profile  shown  in  Figure  4  was  applied  in  both  a  headwind  and 
tail-wind  condition  in  a  simplified  homing  simulation.  The  contribution  of  the 
wind  profile  to  miss  distance  for  either  condition  was  insignificant. 

2.  Selection  of  Inertial  and  Homing  Sensors.  Variations  in  several  sensor 
anomalies  were  investigated  for  effect  on  miss  distance.  These  were:  inertial 
instrument  errors  (i.  e. ,  effect  on  dispersion  at  target  acquisition),  TV  tracker 
bias  errors  and  output  noise  level,  and  tracker  blind  range.  Investigations  of 
the  various  effects  on  proportional  and  pursuit  guidance  were  conducted  primarily 
with  the  analog  simulation;  the  digital  program  was  used  to  provide  error  infor¬ 
mation  for  optimal  guidance  and  checks  on  the  analog  results. 

a.  Effects  of  Dispersion  at  Target  Acquisition.  There  are  two  important  con¬ 
siderations  affecting  allowable  initial  dispersions  at  the  beginning  of  the  terminal 
phase  after  target  acquisition;  the  time  required  to  correct  the  ASM  heading 
errors  with  respect  to  the  target,  and  the  homing  sensor  target  acquisition  capa¬ 
bility.  Initial  dispersion  of  up  to  twice  the  nominal  values  (nominals  were  ap¬ 
proximately  1500  feet)  had  no  effect  on  miss  distance  for  optimal  guidance  for  the 
minimal  acquisition  range  of  15,000  feet.  Examination  of  the  analog  simulation 
data  indicated  that  the  remaining  range  to  the  target,  when  the  effects  of  initial 
errors  were  corrected  by  homing  guidance,  was  2000  feet.  It  is  expected  from 
these  results  that  initial  dispersion  of  the  magnitude  considered  will  have  a  minor 
effect  unless  the  acquisition  range  is  decreased  significantly  from  15,000  feet. 

Initial  dispersion  may  have  a  very  significant  effect  on  the  acquisition  of  the 
target  with  the  homing  sensor.  The  effect  of  initial  position  and  velocity  vector 
errors  must  be  considered  within  the  limitations  in  homing  sensor  field  of  view, 
acquisition  range,  and  allowable  time  for  acquisition.  The  study  of  the  target 
acquisition  was  not  conducted  in  this  program. 
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b.  Effect  of  Tracker  Bias  Errors.  Only  pursuit  and  optimal  guidance  were 
investigated  for  bias  effects.  It  was  assumed  that  gimbaled  tracker  bias  errors 
would  be  removed  during  acquisition  with  proportional  guidance,  because  pro¬ 
portional  guidance  would  be  implemented  using  a  line-of-sight  rate  signal.  (The 
method  of  implementation  was  discussed  in  Section  III. ) 

As  indicated  in  Figure  9,  other  ASM  errors  were  nominal.  (The  contribution 
of  the  nominal  tracker  noise  of  l.C  milliradian  to  miss  distance  is  less  than 
0.  5  foot,  so  bias  error  effects  are  essentially  isolated. ) 

The  effect  of  angular  bias  errors  in  tracker  look  angle  is  shown  in  Figure  9. 

The  bias  error  indicated  was  applied  simultaneously  to  both  the  pitch  angle  and 
the  yaw  angle;  results  are  shown  in  terms  of  analog  average  miss  distance. 

Cross- range  and  down-range  errors  are  comparable.  Pursuit  guidance  errors 
are  approximately  proportional  to  bias  error.  The  miss  distances  for  pursuit 
guidance  are  larger  than  those  that  would  be  expected  from  a  simple  propagation 
of  the  bias  error  over  the  blind  range.  This  larger  error  is  caused  by  the  in¬ 
herent  dynamic  characteristic  of  pursuit  guidance  that  causes  it  to  lag  a  moving 
target,  because  a  pointing  error  is  required  to  generate  a  guidance  command. 

The  tracker  bias  error  causes  the  target  to  have  an  apparent  velocity.  Results 
for  pursuit  guidance  errors  obtained  from  the  digital  linear  analysis  were 
comparable. 

Bias  errors  of  up  to  3  degrees  had  no  effect  on  optimal  guidance  miss  distance. 

c.  Effect  of  Tracker  Noise.  Expected  values  for  TV  tracker  noise  (representing 
line-of-sight  angle  errors)  are  from  1  to  3  milliradians  (rms).  A  range  of  up  to 
50  milliradians  was  examined.  The  effect  of  tracker  noise  on  miss  distance  is 
shown  in  Figure  10  and  Table  IV.  The  analog  data  in  Figure  10  is  presented  as 
average  data  plus  the  rms  deviation  from  the  average.  The  rms  deviation  is 
represented  by  the  shaded  areas.  As  discussed  in  Section  HI,  the  error  data 

"for  each  point  was  obtained  from  the  reduction  of  miss  distance  from  15  to  20 
analog  runs.  The  linear  analysis  data  presented  in  Table  IV  is  the  lv  deviation 
from  the  target. 


32 


•  Long-Range  Semiballistic  Trajectory 

•  Nominal  Dispersion  at  Target  Acquisition 

•  Slant  Range  at  Target  Acquisition  -  15,000  Ft. 

•  Tracker  Bias  -  0  Degree 

•  Blind  Range  (Guidance  Termination)  -  500  Ft. 

•  Sti  1 1  Atmosphere 
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Cross-Hatched  Are«>  Indicates  -  RMS 
Deviation  from  Average  Value 


TRACKER  NOISE  MILLIRADIANS  (RMS) 


Figure  10:  EFFECT  OF  TRACKER  N0I SE  -  LONG-RANGE  TRAJECTORY 


Tracker  noise  was  represented  in  the  analog  simulation  by  a  noise  generator 
that  supplies  an  approximate  white  noise  output  in  the  frequency  range  from  0 
to  35  cps,  while  the  statistical  description  of  tracker  noise  used  in  the  digital 
analysis  included  all  frequencies.  Other  conditions  were  nominal  as  shown  in 
Figure  10  and  Table  IV,  and  were  introduced  in  the  manner  presented  in  Section 
HI. 

The  effect  of  noise  levels  beyond  the  nominal  value  of  1  milliradian  (rms)  are  not 
shown  for  proportional  guidance  in  Figure  10  because  higher  levels  of  noise 
completely  saturated  the  analog  simulation.  The  pursuit  guidance  simulation 
provided  slightly  better  inherent  system  filtering  and  levels  to  5  milliradians 
were  tried.  This  situation  is  comparable  to  what  could  be  encountered  in  a  real 
system.  If  sensor  outputs  with  high  noise  levels  are  not  filtered  so  they  contain 
only  the  frequency  spectrum  required  for  control ,  the  high  frequency  components 
can  saturate  the  autopilot.  For  the  values  investigated  on  the  analog,  the  noise 
does  not  have  a  large  effect  on  miss  distance. 

Effects  of  tracker  noise  beyond  10  milliradians  were  evaluated  with  the  digital 
linear  analysis  for  proportional,  pursuit,  and  optimal  guidance.  The  digital 
results  are  shown  in  Table  IV.  An  interesting  result  is  that  optimal  guidance 
does  not  produce  significantly  better  results.  Apparently  the  inherent  filtering 
of  the  ASM  syt-.em  dynamics  is  almost  as  good  as  the  Kalman  optimal  filter 
with  respect  to  resultant  system  error. 

d.  Effect  of  Tracker  Blind  Range.  The  effect  of  tracker  blind  range  is  shown 
in  Figures  11  through  14.  Data  are  shown  for  two  autopilot  command  techniques 
following  guidance  termination:  holding  the  last  commanded  normal  accelerations 
during  blind  range  traverse,  and  commanding  zero  normal  acceleration  (z^,  =  0) 
during  blind  range  traverse. 

As  in  the  other  analog  studies,  several  runs  were  made  for  a  given  trajectory 
(i.  e. ,  long-range  or  skip)  with  the  nominal  aet  of  initial  position  deviations. 
Because  of  the  initial  deviations,  a  maneuver  is  required  to  Impact  the  target 
even  in  the  absence  of  disturbances.  Atmospheric  turbulence  of  10  fps  rms  was 
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TABLE  IV 

EFFECT  OF  TV  TRACKER  NOISE  -  \a  MISS  DISTANCE 


RMS  NOISE  LEVEL 
MILLIRADIAN 

*  MISS  DISTANCE  (FEET) 

PROPORTIONAL 

GUIDANCE 

PURSUIT 

GUIDANCE 

OPTIMAL 

GUIDANCE 

1 

3.1 

3.3 

2.2 

10 

15 

13.3 

14.3 

50 

72 

64.8 

59.4 

CONDITIONS: 

NOMINAL  DISPERSION  AT  ACQUISITION 
SLANT  RANGE  AT  ACQUISITION  —  15,000  FEET 
TRACKER  BIAS  —  1 .74  MILLIRADIANS 
TRACKER  BLIND  RANGE  —  500  FEET 


*  ERROR  MEASURED  PERPENDICULAR  TO  TRAJECTORY 
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Figure  13:  EFFECT  OF  BLIND  RANGE  FOR  PROPORTIONAL  GUIDANCE 
LONG-RANGE  TRAJECTORY 
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included  as  a  disturbance.  Two  effects  are  noted  when  blind  range  is  increased. 
The  average  (or  arithmetic  mean)  miss  distance  is  increased.  This  effect 
appears  to  be  primarily  associated  with  the  initial  deviation.  Some  minimal 
time  is  required  to  reduce  initial  dispersion  and  to  get  the  velocity  vector 
directed  at  the  target.  This  time  varies  between  trajectories  and  guidance  laws, 
but  generally  for  the  15,000  feet  of  slant  range  at  acquisition,  is  accomplished 
with  a  minimum  of  2000  feet  of  slant  range  remaining.  Figures  11  through  14 
show  that  appreciable  increases  in  average  miss  distance  do  not  occur  until  the 
blind  range  is  extended  beyond  1500  feet.  The  second  effect  is  the  rms  deviation 
about  the  average  miss  distance.  (The  rms  deviation  is  indicated  by  the  shaded 
areas.)  This  appears  to  be  the  result  of  dispersion  in  the  acceleration  command 
at  the  beginning  of  the  blind  range  caused  by  turbulence  and  tracker  noise.  If 
the  acceleration  command  is  zeroed,  rather  than  held  at  the  value  occurring  at 
the  beginning  of  the  blind  range,  the  rms  deviation  is  reduced.  The  only  excep¬ 
tions  to  this  rule  occur  at  blind  ranges  of  2000  feet  or  less.  At  these  blind 
ranges,  the  difference  in  rms  deviation  between  holding  the  last  command  and 
setting  the  command  to  zero  is  very  small  (less  than  1. 5  feet) ,  and  the  choice 
of  one  mode  of  command  over  the  other  is  not  critical. 

3.  Selection  of  Autopilot.  This  study  was  limited  to  the  investigation  of  one 
type  of  autopilot.  A  normal  acceleration  autopilot  was  selected  because  it  pro¬ 
vides  more  direct  flight  path  control  than  other  types  (such  as  an  attitude  auto¬ 
pilot).  By  limiting  commands  to  the  normal  acceleration  autopilot,  a  means  of 
meeting  structural  load  limit  requirements  is  also  obtained.  The  normal  accel¬ 
eration  autopilot  is  compatible  with  a  gimbaled  homing  sensor.  A  body  fixed 
homing  sensor  may  require  a  different  type  of  autopilot  if  it  has  a  small  field 
of  view. 

The  effect  on  miss  distance  of  variations  in  the  autopilot  response  frequency 
was  investigated  on  a  simplified  3  -degree -of -freedom  simulation  and  with  the 
6 -degree -of -freedom  analog  simulation.  The  simplified  studies  used  a  third 
order  autopilot  with  well  damped  quadratic  poles  at  20  radians/sec.  With  slant 
ranges  and  initial  position  offsets  comparable  to  those  used  in  the  6 -degree -of  - 
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freedom  error  studies,  the  first  order  pole  could  be  reduced  to  2  radians/sec 
with  no  effect  on  miss  distance.  A  response  of  one  radian/sec  resulted  in  sig¬ 
nificantly  larger  miss  distance  so  the  autopilot  response  requirement  was  initially 
set  at  2  radians/sec. 

The  effect  of  autopilot  response,  as  obtained  from  the  6 -degree -of -freedom 
analog  simulation,  is  shown  in  Figure  15.  The  abrupt  increase  in  miss  distance 
below  6  radlans/sec  is  caused  by  autopilot  instability  rather  than  "looseness" 
of  control.  In  the  final  detailed  autopilot  design,  the  ASM  rigid  mode  poles  were 
significantly  lower  than  the  20  radians/sec  anticipated  in  the  preliminary  studtes. 
Thus,  with  guidance  loops  closed,  it  was  not  possible  to  lower  the  response 
below  6  radians/sec  without  the  ASM  becoming  unstable.  Actually,  in  the 
nominal  Type  1  normal  acceleration  autopilot  simulated,  the  time  response 
during  the  terminal  phase  of  the  nominal  trajectories  ranges  from  0. 1  to  0. 17 
second.  These  values  were  adequate  for  the  range  of  initial  dispersions  and 
slant  ranges  investigated  in  this  study.  However,  as  mentioned  in  the  discussion 
of  blind  range,  ASM  response  (time  to  remove  initial  errors),  slant  range  at 
acquisition,  and  blind  range  are  very  much  interrelated.  Shorter  slant  ranges, 
larger  initial  errors,  and  longer  blind  ranges  may  require  faster  autopilot 
response. 

The  effect  of  autopilot  steady  state  error  on  accuracy  wan  not  investigated. 

Steady  state  autopilot  errors  tend  to  look  like  changes  in  guidance  gain.  Early 
in  the  study  it  was  believed  that  it  would  be  very  difficult  to  obtain  small  impact 
errors,  and  that  appreciable  steady  state  autopilot  errors  could  not  be  tolerated 
because  of  anticipated  tight  tolerances  on  autopilot  gain.  From  the  evaluation  of 
guidance  gain  in  paragraph  l.b.  of  this  section  it  can  be  seen  that  not  much 
tolerance  is  available  for  proportional  guidance,  while  pursuit  guidance  is  not 
quite  as  critical.  A  "Type  O"  autopilot  (with  constant  gain)  that  was  studied  had 
steady  state  errors  to  step  commands  which  varied  from  2  to  16%  during  the 
homing  phase.  The  missile  configuration  that  was  m,od  in  the  study  had  near 
neutral  aerodynamic  stability.  Therefore,  It  was  relatively  easy  to  add  Integral 
forward  loop  compensation  to  remove  this  steady  state  error,  and  at  the  same 
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MISS  DISTANCE  (FT) 
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•  Nominal  Dispersion  at  Target  Acquisition 

•  Slant  Range  at  Acquisition  =  15,000  Ft. 

•  Tracker  Noise  =  0.125  Milliraaians  (rms) 

•  Blind  Range  (Guidance  Termination)  =  500  Ft 

•  Long-Range  Traiectory 

•  Atmospheric  Turbulence  =  10  Ft/Sec  (rms) 
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EQUIVALENT  1ST  ORDER  RESPONSE  FREQUENCY  -  RAD /SEC 


Figure  15:  EFFECT  OF  FIRST  ORDER  RESPONSE  FREQUENCY  - 

PROPORTIONAL  GUIDANCE;  LONG-RANGE  TRAJECTORY 


time  satisfy  the  requirements  for  adequate  gain  and  phase  margins.  The  resulting 
"Type  1"  system  had  adequate  response  and  damping  as  noted  by  the  near  0. 1 
second  response  during  the  terminal  phase.  This  relatively  easy  solution  is  not 
possible  in  all  instances,  and  the  effect  of  steady  state  autopilot  error  may  be 
critical  for  other  configurations. 


SECTION  V 


ASM  OPTIMAL  GUIDANCE 

An  optimal  guidance  law  for  a  given  ASM  is  one  that  yields  the  minimum  mean 
square  miss  distance.  It  gives  the  best  possible  performance  attainable  with 
the  ASM  and  can  serve  as  a  standard  of  comparison  for  other  guidance  laws. 

A  well-developed  theory  exists  for  optimal  control  of  linear  systems  with  ran¬ 
dom  errors.  However,  real  missile  systems  are  nonlinear,  and  an  adequate 
theory  of  optimal  control  is  not  available  for  nonlinear  systems  with  random 
errors.  In  this  study,  the  realistic  nonlinear  ASM  equations  of  motion  are 
linearized  about  a  nominal  trajectory.  Linear  random  control  theory  is  used 
to  compute  the  optimal  control  law  for  the  linearized  ASM  equations  of  motion. 
Because  these  equations  are  an  accurate  approximation  of  actual  nonlinear 
equations,  this  control  law  is  a  good  approximation  of  the  optimal  nonlinear 
ASM  control  law.  The  ASM  system  and  system  errors  considered  are 
described  in  Section  EH. 

The  optimization  theory  on  which  the  ASM  optimal  guidance  was  based  states  that 
the  optimal  guidance  law  consists  of  two  parts.  One  .is  a  Kalman  optimal  filter 
that  supplies  the  best  possible  estimate  of  position  and  velocity  coordinates  in  the 
presence  of  sensor  and  system  noise.  The  other  is  an  optima  ,  linear  feedback 
control  law  that  converts  these  estimates  into  steering  commands.  This  type  of 
guidance  system  is  represented  in  Figure  16.  The  Kalman  optimal  filter  and  the 
optimal  feedback  control  law  contain  sets  of  time- varying  coefficients  that  are 
dependent  on  the  nominal  trajectory.  These  coefficients  are  precomputed  and 
stored  in  the  onboard  ASM  computer.  The  performance  criterion  for  the  control 
law  was  the  minimization  of  the  sum  of  the  mean  square  miss  distance  and  the 
integral  of  a  quadratic  function  of  the  control  vector.  The  form  of  the  guidance 
law  depends  on  the  type  of  sensors  used  in  the  ASM.  Optimal  guidance  laws  were 
considered  for  an  ASM  with  both  a  TV  tracker  and  an  inertial  platform  (IMU)  and 
with  only  a  TV  tracker. 
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Figure  16:  OPTIMAL  TERMINAL  GUIDANCE 


The  succeeding  sections  discuss  in  detail  development  of  the  optimal  guidance 
law  and  sizing  an  airborne  computer  required  to  implement  it,  and  draw  some 
conclusions  regarding  the  performance  of  the  optimal  guidance  law. 

1.  ASM  Description.  The  nominal  ASM  characteristics  are  described  in  Section 
HI.  The  ASM  equations  of  motion  and  sensor  equations  as  used  in  the  develop¬ 
ment  of  the  optimal  guidance  law  are  presented  in  this  section. 


The  ASM  was  represented  as  a  point  mass.  The  autopilot  was  considered 
ideal;  i.e.,  normal  acceleration  commands  are  transcribed  directly  into 
acceleration  of  the  ASM.  Both  body  fixed  and  gimbaled  TV  trackers  were 
considered  in  providing  line-of-sight  error  information  to  the  optimal 
guidance  laws.  ASM  acceleration  information  was  obtained  with  an  inertial 
platform. 

a.  Equations  of  Motion.  Because  the  terminal  portion  of  the  ASM  trajectory  is 
only  a  few  miles,  the  equations  of  motion  of  the  ASM  are  those  of  a  vehicle  with 
lift  and  drag  moving  in  a  constant  gravitational  field.  The  ASM  coordinate  vari¬ 
ables  are  indicated  in  Figure  17.  Expressed  in  flight  path  coordinates,  the  point 
mass  equations  of  the  ASM  are: 


x  =  ♦  V  cos  y  cos  X 


y  =  +  V  cos  y  sin  X 
z  =  +  V  sinT 

.  D  (a  +  Aoew,  v  +  AVW,  z) 

V  =  - - -  g  sin  y 


X  = 

y  = 


L  («  +  Ac*w,  V  +  AVW,  z) 
mV  cos  y 

L  (or  +  Aaw,  V  +  Ayw,  Z) 
mV 


sin  (P  +  A0W) 
cos  (P  +  A  fa#)  - 


& 

V 


cos  y . 


(10) 

(11) 

(12) 

(13) 

(14) 

(15) 


Note  that  the  angle  of  attack,  a,  and  bank  angle,  P,  defined  in  Figure  17,  are  dif- 
erent  variables  than  those  usually  used  in  aerodynamic  studies.  In  this  analysis, 
a  is  the  missile  body  total  angle  of  attack  and  p  the  bank  angle  relative  to  the 
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alrstream.  Tbe  aerodynamic  drag  (D)  and  lift  (L)  forces  are  expressed  ss 
functions  of  velocity,  angle  of  attack,  and  Mach  number  with  provisions  for  the 
inclusion  of  wind  effects  incorporated  by  the  perturbed  expressions  V  +  AVW, 

®+  Attw,  and  0  + 

The  functional  form  of  the  expressions  L  (a,  V,  z)  and  D  (a,  V,  z)  is: 

L  (a,  v.  Z)  - 1  p  rn  v2  cL  fa,  i  (16) 

D  («.  V,  z)  -  §  P  (2)  V2  CD  L,  ^ 

In  these  expressions  the  air  density,  p  ,  is  expressed  as  a  function  of  altitude  by 
standard  ARDC  tables;  Cl  and  Cq  are  tabulated  functions  of  angle  of  attack  and 
Mach  number  obtained  from  wind  tunnel  data  and  the  velocity  of  sound,  V8,  is 
expressed  as  a  tabulated  function  of  altitude. 

b.  TV  Tracker.  TV  trackers  measure  line -of- sight  angles  (LOS)  to  the  target. 

Either  body  fixed  or  gimbaled  TV  trackers  can  be  considered  for  instrumentation 
with  the  optimal  guidance  system. 

] 

t 

A  body  fixed  TV  tracker  measures  the  azimuth  and  elevation  of  the  LOS  for  a 

coordinate  system  fixed  in  the  missile,  aligned  with  the  missile  axis.  This  can  j 

be  combined  with  missile  attitude  information  to  provide  measurements  of  LOS  I 

azimuth  and  elevation  for  an  inertial  coordinate  system.  There  will  be  two  typer- 

of  error  in  this  measurement:  a  bias  error  due  to  a  misalignment  of  the  sens;: 

and  a  noise,  which  will  be  assumed  to  be  a  white  Gaussian  noise,  due  to  the  ei 

fects  of  target  background  and  the  measurement  process.  If  i>  and  are  tbe 

measured  LOS  azimuth  and  elevation,  the  equations  for  and  $  may  be  written: 

♦  =  tan-1  ^  +  bx  +  Nj,  (18) 

+  b2  +  N2  (19) 
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in  which  bj  and  bg  are  the  random  biaaes  and  and  N2  are  white  noises. 

In  a  basic  implementation  of  a  gtm baled  tracker,  the  gimbals  are  inertially  rate 
stabilized.  This  rate  stabilized  system  is  commanded  to  point  at  the  target  with  an 
optical  sensor  loop  that  measures  the  tracker  head  deviation  from  the  target  line  of 
sight.  As  described  in  Section  III,  this  tracker  was  assumed  to  have  a  quad¬ 
ratic  response.  If  the  azimuth  and  elevation  of  the  target -are  designated  tjj  and 
t)2.  the  following  equations  represent  the  gimbaled  tracker  system: 

*/'l  +  Cx  \  =  C2(tan_1  *  +Nj_  -  t^)  ,  (20) 


n2  +C3n2=  C4 


'(V7T72)  * N2 '  ”2} 


The  arctangent  functions  represent,  in  target  coordinates,  the  optical  sensor 
measurements  of  tracker  head  deviation  from  the  target,  and  and  C2  repre¬ 
sent  the  traditional  quadratic  damping  and  natural  frequency  terms. 

The  quantities  yx(t)  =  hi  (t)  and  y2(t)  =  f)2(t)  are  measured,  it  is  possible  to 
express  the  LOS  angles  4>  and  4>  in  terms  of  these  measurements  in  the  following 
way.  From  the  definition  of  y-^ft)  and  y^(t) 

^(t)  =  \<to)  +  f  yi(s)  ds  (22) 


V>  =  ^(tQ)  +  /  y2(s)  ds. 


It  follows  that 


-1  y  i  /•* 

tan  *  +  Nx  =  —  C^i  +  C^yp  +  J  y1(s)  ds  +  l?^) 


N2  ‘^*2+  C3V  */  y2(s»  *  +  Vo> 

\f  x  +  y  /  4  t 
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The  Initial  alignment,  tj^  <t^) ,  ij^ft^,  of  the  seeker  axis  is  known  except  for  the 
initial  alignment  bias  errors  bj  and  b2;  that  is,  the  quantities  +  bj, 

rj^fto)  +  l>2  can  be  measured.  This  Implies  that  and  ♦  may  be  expressed  in 


terms  of  the  measurement  quantities: 

-t  r* 

yl*  y2*  yi*  y2’  J  yl(8)  d8'  J  y2(S)  d8’  +  bi*  aaA%^  +  b2  P6) 

in  t  - 


Thus,  the  optimal  filter  for  a  gimbaled  TV  seeker  may  also  be  designed,  assuming 
that  <t>  and  4>are  measured. 


Note  that  in  the  gimbaled  tracker,  the  tracker  input  is  assumed  to  be  the  differ¬ 
ence  between  the  line-of-sight  angle  plus  white  noise  and  the  tracker  axis  angle. 
The  only  biases  considered  are  those  introduced  by  error  in  measurement  of  the 
initial  alignment  77^ <tQ)  and  (tQ)  of  the  tracker  axis.  If  there  are  biases  in  the 
input  to  the  tracker,  this  can  be  handled,  using  the  same  set  of  equations,  by 
reinterpreting  ^(t)  and  ?^><t)  as  the  tracker  axis  angles  plus  these  biases.  Be¬ 
cause  the  derivatives  of  fjjft)  and  Tjgft)  are  measured,  Equations  27  and  28  hold 
without  change  in  this  situation. 

c.  Inertial  Platform.  A  conventional  inertial  platform  attempts  to  measure 
missile  accelerations  in  a  coordinate  system  fixed  in  space.  If  a^ft),  a^), 
a  (t)  are  the  measured  accelerations,  missile  position  and  velocity  are  obtained 


by  integrating  the  equations 


(29) 


with  initial  conditions  given  by  initial  estimates  of  position  and  velocity.  The 
"g"  in  the  equation  represents  the  bias  needed  to  represent  the  effect  of  the 
gravitational  acceleration  on  the  vehicle. 

The  two  sources  of  error  in  the  computation  of  position  and  velocity  are  in  the 
initial  estimates  of  position  and  velocity ,  and  in  the  acceleration  measurements.  The 
major  source  of  error  in  the  computation  of  missile  positions  and  velocities  from 
integrating  Equation  29  is  in  the  knowledge  of  the  missile's  initial  position  and 
velocity.  The  errors  in  measurement  caused  by  the  rotation  of  the  platform, 
scale  factor,  and  bias  are  of  a  lower  order  of  magnitude.  This  conclusion  is 
substantiated  by  the  inertial  platform  error  analysis  studies  discussed  in  Appen¬ 
dix  I.  For  this  reason,  the  acceleration  measurements  a^),  a^ft),  a^(t)  will  be 
treated  as  exact.  This  assumption  considerably  reduces  the  complexity  of  the 
optimal  filter  implementation. 


If  measurement  bias,  scale  factor,  and  platform  tilt  errors  must  be  taken  into 
account,  the  true  missile  accelerations  and  measured  missile  accelerations  are 
related  by  the  formula: 


-  bx 
a2(t)  -  b2 
,a3(t) 


(30) 


In  this  formula,  bp  b2,  and  bg  are  measurement  bias  errors  and  C  is  a  matrix 
that  expresses  the  scale  factor  errors  in  measurements  and  the  amount  the  iner¬ 
tial  platform  has  rotated  out  of  alignment  with  its  desired  inertial  coordinate 
system.  To  take  into  account  the  random  effects  of  scale  factor  and  bias  errors 
in  the  optimal  filter,  the  elements  of  the  matrix  C  and  b  ,  b0,  b  would  have  to 
be  introduced  as  new  state  variables.  This  greatly  increases  the  order  of  the 
optimal  filter  equations. 
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2.  Optimal  Control  of  Linear  Systems  with  Gaussian  Noises.  The  techniques 
used  to  compute  the  optimal  guidance  law  will  be  based  on  the  theory  given  In 
Reference  3.  These  results  are  summarized  in  this  section. 


Let  a  linear  system  be  defined  by: 

x  =  F (t)x  +  G(t)u  +  C(t)  ,  (31) 

where  x(t)  is  an  n-vector  of  state  variables,  and  u(t)  an  m-vector  of  control 
variables;  F(t)  and  G(t)  are  known  n  x  n  and  n  x  m  matrices;  £(t)  is  an  n- 
vector  of  Gaussian  white  noises  with  zero  means.  The  covariance  matrix  of 
C(t)  is  aesumed  known  and  given  by: 

E  [c  (t)  C  fr)T}  =  Q(t)  (32) 

The  symbol  E  |  J  denotes  the  expected  value  of  the  quantity  in  the  brackets. 

3  un  quantity  6(i  -  T)  is  the  Dirac  6-function.  Suppose  the  p-vector 

y(t)  =  H(t)  x(t)  +  rj(t)  (33) 


is  observed.  H(t)  is  a  known  p  x  n  matrix  and  r)(t)  a  p-vector  of  Gaussian  white 
noises  with  zero  means  and  known  covariance  matrix 

E(rXt)  rj(nT)=  R(t)  6<t  -  r>.  (34) 

In  addition,  the  cross  correlations  of  C(t)  and  T)(t)  are  assumed  to  be  given  by: 

Eject)  r)(r)T|  =  S(t)  6(t  -  r).  (35) 

Let  tQ  be  the  initial  and  t^  the  final  times.  Let  the  performance  index  be: 


n  rl  m 

Z  X)  xi  +  J  Z  bi,  j  ui  &  uj  <s>  68 


(36) 


The  first  term  in  the  performance  index  is  a  quadratic  function  of  the  state  vector 
at  the  final  time,  while  the  second  term  in  the  time  integral  of  a  quadratic  func¬ 
tion  of  the  control  vector.  This  performance  index  penalizes  final  error  and 
control  effort  used  to  achieve  the  final  state. 
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Consider  the  problem  of  choosing  u(t)  as  a  function  of  the  past  observations: 

u(t)  =  u[y(n),u(s),t0  s  s  <  t  s  txJ  (37) 

so  that  the  performance  index  is  minimized. 

THEOREM:  Let  the  matrices  R  and  M  =  (b^)  be  positive  definite.  Then  the 
optimal  control  law  is  given  by: 

u(t)  =  -M_1  GT  Ux  (38) 

in  which  x,  an  estimate  of  x,  is  the  n-vector  output  of  the  filter  whose  input 
is  the  observation  vector,  y(t),  and  whose  equations  are  given  by: 

x  =  (F  -  SR_1  H  -  EHT  R_1  H)  x  -  GM’1  GT  Ux  +  (SR-1  +  £HT  R_1)y 

(39) 


The  quantity  U  is  an  n  x  n  matrix  which  satisfies  the  differential  equations 

U  =  UGM_1  GT  U  -  UF  -  FT  U  (40) 

with  the  terminal  condition 


Uty  =  (a..)  (41) 

The  quantity  £  is  an  n  x  n  matrix  which  satisfies  the  differential  equation 
£  =  -  £HT  R-1  H£  +  (F  -  SR_1  H)  £  +  £  <FT  -  HT  R_1  ST)  +  Q  -  SR_1  ST  (42) 
The  matrix  £  must  satie  fv  the  initial  condition 

£(to)  =  E|x(to)x(to)TJ  (43) 

Remarks 

The  optimal  feedback  control  law  for  the  deterministic  system  with  equation 

x  =  F(t)  x  +  G(t)  u  (44) 

and  performance  criteria  given  by 

n  *1  m 

Z  ai,  j  *i  %)  x^i)  +  J  2  bi  i  ut  (s)  u.(s)  da  (45) 

U=1  J  tQ  i, j=l  1 
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Is  given  by 


u(t)  =  -  M_1  GT  Ux  (46) 

in  which  U  is  a  solu+ion  of  Equation  40.  The  optimal  least  squares  estimate  of 
the  state  variables,  that  is  the  Kalman-Wiener  estimate,  for  the  system 

x  =  F(t)  x  -  G(t)  M_1  G(t)T  U(t)  x  +  C(t)  (47) 

with  measurements 

y(t)  =  H(t)  x(t)  +  r)(t)  (48) 

is  given  by  x(t)  which  is  the  solution  of  Equations  39  and  42. 

These  two  remarks  show  that  the  results  of  the  theorem  can  be  stated  intuitively 
by  saying:  The  optimal  filter  may  be  computed  by  ignoring  the  noises  and  com¬ 
puting  the  optimal  feedback  control  law,  Equation  46,  for  the  deterministic  sys¬ 
tem.  Substitute  this  control  law  for  the  control  law  of  the  system  and  compute  the 
Kalman-Wiener  optimal  estimates  of  the  state  variables,  Equation  39.  The  op¬ 
timal  control  law  is  then  obtained  by  substituting  these  estimates  for  the  state 
variables  in  the  deterministic  optimal  feedback  control  law. 

The  matrix  £  (t),  which  is  the  solution  of  Equation  42,  and  is  deed  in  defining  the 
filter  Equation  39,  has  another  important  property.  The  matrix  E  (t)  is  the  co- 
variance  matrix  of  the  difference  between  the  state  vector,  x(t),  and  the  state 
vector  estimate,  x(t). 

3.  Application  of  Optimal  Control  Theory  to  ASM  Guidance — In  this  subsection 
the  theory  of  the  previous  subsection  is  applied  to  determine  the  equations  for 
the  ASM  optimal  guidance  law.  The  approximations  of  the  system  model  that 
are  necessary  to  apply  the  theory  to  this  problem  are  carried  out.  The  theory 
of  Section  V,  2.  implies  the  optimal  guidance  law  consists  of  two  parts:  a  Kal¬ 
man  optimal  filter,  and  an  optimal  feedback  control  law.  The  equations  that 
determine  each  of  these  are  specified.  The  equations  for  the  optimal  filter 
are  obtained  for  an  ASM  that  uses  both  a  TV  tracker  and  an  Inertial  platform,  and 
an  ASM  that  uses  only  a  TV  tracker.  Some  of  the  relative  advantages  of  these 
types  of  filters  are  discussed.  In  one  of  these  filters  the  equations  call  for  the 
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derivative  of  a  measurement  whose  derivative  contains  white  noise.  Because  it 
is  not  practical  to  instrument  a  system  to  compute  such  a  derivative,  a  proce¬ 
dure  is  given  that  avoids  the  problem. 

a.  Linearization  of  the  Equations  of  the  ASM.  The  theory  of  the  preceding  sec¬ 
tion  applies  to  linear  systems.  In  order  to  apply  this  theory  to  the  ASM,  nom¬ 
inal  trajectories  were  selected  and  the  equations  of  motion  linearized  about 
them.  The  two  different  nominal  trajectories  that  were  used  are  discussed 
in  Section  in. 

Equations  10  through  15  were  linearized  about  the  nominal  trajectories  by  approxi¬ 
mating  the  ASM  trajectory  by  the  nominal  trajectory  plus  variations  from  the 
nominal  trajectory.  The  variations  are  solutions  of  the  equations: 


in  which  the  matrices  F(t),  G(t)  and  D(t)  are  respectively  matrices  of  partial 
derivatives  of  the  right  sides  of  Equations  10  through  15  with  respect  to  the  state 
variables,  the  control  variables,  and  the  variations  &VW,  due  to 

wind.  In  the  matrices,  the  partial  derivatives  are  expressed  as  functions  of 
time  by  substituting  values  of  the  state  and  control  variables  of  the  nominal 
trajectory  at  the  corresponding  time  in  these  expressions.  These  three  matrices 
of  partial  derivatives  are  given  in  Figures  18,  19,  and  20.  The  approximation 
of  the  trajectory'  by  a  nominal  plus  variations  from  the  nominal  is  a  standard 
technique  for  approximating  a  nonlinear  system  by  a  linear  system.  A  good 
approximation  is  obtained  if  the  actual  trajectory  is  close  to  the  nominal  trajectory. 

To  illustrate  the  variational  equations  more  concretely,  the  equation  for  6x  could 
be  obtained  as  follows:  Let  x(l),  y(t),  z{t),  v(t),  \(t),  y(t),  ar(t),  /j(t)  denote  the 
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Figure  18:  MATRIX  F  (t)  OF  INFLUENCE  COEFFICIENTS  OF  THE 
STATE  VARIABLES  IN  THE  LINEARIZED  EQUATIONS 
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state  and  control  variables  of  the  nominal  trajectory.  The  first  rows  of  F(t), 

G(t),  and  D(t)  are  the  vectors 

[o,  0,  0,  cos  y  (t)  cos  x  (t),  -  V(t)  cos  V  (t)  sin  x(t)»  -  V(t)  sin  y(t)  cos  x  (t)]  ,  (50) 

[o,  o]  and  [o,  0,  o] 

Hence  the  equation  for  6x  is  given  by: 

6x  =  cos  y(t)  cos  x  (1)  6V(t)  -  V(t)  cos  y(t)  sin  X  (t)  §X  (t)  (51) 

-  V(t)  sin  y(t)  cos  X  (t)  6v(t) 

b.  Performance  Index.  To  use  the  theory  of  Section  V.2,  the  performance 
index  must  be  a  quadratic  function  of  the  state  variables  at  the  final  time  plus 
the  integral  of  a  positive  definite  quadratic  form  in  the  control  variables.  The 
true  miss  distance  of  the  ASM  is  the  distance  along  the  ground  from  the  impact 
point  to  the  target.  The  final  time  is  considered  to  be  the  time  when  the  nominal 
trajectory  hits  the  target.  The  magnitude  of  the  position  deviation  of  the  actual 
trajectory  from  the  nominal  trajectory  at  this  time  does  not  give  the  true  miss 
distance.  Suppose  the  position  deviation  at  this  time  is  resolved  into  a  vector 
parallel  to  the  nominal  trajectory  and  a  vector  perpendicular  to  the  nominal  tra¬ 
jectory.  The  true  miss  distance  may  be  approximated  by  the  magnitude  of  the 
deviation  perpendicular  to  the  nominal  trajectory  times  the  cosine  of  the  angle 
the  nominal  trajectory  makes  with  the  ground. 

A  trajectory  with  a  nonzero  component  of  deviation  parallel  to  the  nominal  trajec¬ 
tory  will  have  approximately  the  same  miss  distance  as  one  with  a  zero  com¬ 
ponent.  Because  the  ASM  velocity  vector  will  be  nearly  aligned  with  the  nominal 
trajectory  near  the  target,  this  component  of  position  deviation  will  not  contribute 
significantly  to  the  miss  distance. 

Based  on  this  reasoning,  the  performance  criterion  will  be  taken  to  be  the 
expected  value  of  the  square  of  the  component  of  deviation  perpendicular  to  nom¬ 
inal  trajectory  at  the  final  time  plus  the  integral  of  a  constant  multiple  of  the 
sum  of  squares  of  angle  of  attack  and  bank  angle.  (The  reciprocal  of  the  cosine 
of  the  angle  the  nominal  trajectory  makes  with  the  ground  has  been  incorporated 


into  the  constant  multiple.)  It  is  necessary  to  have  the  integral  term  in  the 
control  variables  in  the  performance  criterion  in  order  to  ensure  that  a  solution 
exists.  The  constant  multiple  of  this  factor  was  chosen  to  make  this  factor 
small  compared  to  the  perpendicular  component  of  deviation.  This  implies  that 
the  dominant  term  in  the  performance  criterion  is  the  square  of  the  perpendic¬ 
ular  component  of  deviation.  This  performance  criterion  is  given  by  the  expected 
value  of  the  expression: 


sin  y(t^)  - 


6z(tj)  cos 


y^i))2  +  Sy^)2 


o 


(52) 


c.  Computation  of  the  Optimal  Feedback  Control  Law.  The  theorem  of  Refer¬ 
ence  3  that  is  stated  in  Section  V.2  asserts  that  the  optimal  feedback  control 
law  for  system  (Ell)  with  performance  criteria.  Equation  52,  is  given  by: 

^6x  \ 

6y 


-1  T 
=  -  M  G  U 


in  which  U  is  the  solution  of  the  matrix  differential  equation 


U  =  UGM*1  GTU  -  UF  -  FTU 


(53) 


(54) 


with  terminal  condition 


U(tx)  - 


sin2  y(t1) 


0 

1 


-  sin  yft^  cos  y(t  )  0 


0 

0 
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0 
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-  sin  y(tx)  cos  ytf^ 


cos2  y^) 
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The  matrix,  M,  of  Equations  53  and  54  is  a  2  x  2  diagonal  matrix  with  the  con¬ 
stant  "a"  of  the  performance  criteria  in  each  diagonal  position.  The  only  unknown 
quantity  in  Equation  53  is  the  matrix  U(t),  which  must  be  computed  by  numerically 
integrating  Equation  54. 

In  Figures  21,  22,  and  23,  typical  nonzero  coefficients  of  the  optimum  control 
law  are  given.  The  variation  in  angle  of  attack  is  a  linear  sum  of  these  coef¬ 
ficients  times  the  variations  in  the  state  variables. 

The  coefficients  illustrated  are  the  elements  of  the  first  row  of  the  matrix 
-M  given  in  Equation  53.  Because  there  is  a  wide  variation  in  the  values 

of  these  coefficients,  they  have  been  plotted  with  a  scale  change. 

The  coefficients  are  relatively  small  until  near  the  time  the  nominal  trajectory 
reaches  the  target.  Part  of  this  behavior  may  be  accounted  for  by  noting  that  if 
there  is  still  a  position  deviation  just  prior  to  the  terminal  time,  a  large  control 
force  is  needed  to  correct  it.  As  the  terminal  time  is  approached,  the  coefficients 
approach  zero.  This  may  be  explained  by  noting  that  very  close  to  the  terminal 
time  there  is  not  sufficient  time  left  to  make  an  appreciable  position  correction. 
Hence,  the  controller  merely  minimizes  the  integral  term  in  performance  cri¬ 
terion,  Equation  36,  by  letting  the  control  approach  zero. 

d.  Optimal  Guidance  Filter.  The  optimal  filter  depends  on  the  type  of  sensors 
used  on  the  missile.  In  the  following  discussion,  the  application  of  the  theory  of 
Reference  3  to  two  combinations  of  sensors  will  be  discussed.  An  optimal  filter 
based  on  measurements  from  both  a  TV  tracker  and  an  Inertial  platform  will 
give  better  estimates  of  position  than  one  based  only  on  a  TV  tracker.  Because 
an  optimal  filter  based  on  a  TV  tracker  alone  would  be  less  expensive  to  imple¬ 
ment,  and  could  be  designed  in  a  simple  way  from  the  theory  in  Section  V.  2, 
both  types  of  filters  are  discussed  in  the  following  sections. 

(1)  Optimal  Filter  Based  on  a  TV  Tracker.  In  Section  V.  1.  b. ,  it  was  shown 
that  either  a  body  fixed  wide  angle  field  of  view  or  a  gimbaled  narrow  field  of 
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WEIGHTING  COEFFICIENT  OF  X-POSITION  ON  ANGLE  OF  ATTACK 


WEIGHTING  COEFFICIENT  OF  FLIGHT  PATH  ANGLE  ON  ANGLE  OF  ATTACK 


view  TV  tracker  may  be  used  to  produce  measurements  of  line-of-sight  azimuth 
and  elevation  angles  with  additive  biases  and  Gaussian  white  noises.  The  equa>- 
tions  of  these  measurements : 


-1 


(p  =  tan  ^  +  by  +  Ny 


4> 


tan 


-1 


V  <>?  +  y2) 


+  bo  +  No 


(56) 


in  which  by  and  b2  are  ranaom  biases  with; 

by  -  b2  =  0 


(57) 


and  Ny  and  N2  are  Gaussian  white  noises  with  mean  zero  and  covariance  matrix; 


R  (t) 


e|ni  (t)2j 

j  N2  (t)  Ny  (t)j 


(58) 


These  line-of-sight  measurement  equations  were  linearized  about  the  nominal 
vrajectory  to  obtain  measurements  that  are  linear  in  the  variations  from  the 
nominal  trajectory.  When  this  is  done  the  equations, 


6i>  =  Ay  (t)  6X  +  A2  (t)  6y  +  by  +  Ny 

6  $  =  By  (t)  6x  +  B2  (t)  6y  +  B3  (t)  6z  +  b2  +  N2 


(59) 


are  obtained.  (t)  and  Bj  (t)  are  the  partial  derivatives  of  the  appro¬ 

priate  arctangent  expressions  evaluated  on  the  nominal  trajectory.  In  order  to 

apply  the  theory  of  Section  V.2,  the  biases  by  and  b2  are  considered  as  extra 

•  • 

state  variables  and  the  equations  by  =  0  and  b2  =  0  are  added  to  the  equations  of 
motion.  Let  H  (t)  denote  the  matrix 

Ay  (t)  A2  (t)  0  0  0  0  1  0 
By  (t)  B2  (t)  B3  (t)  0  0  0  0  1 

Let  Avw,  Aaw,  Af^be  the  vector  of  white  noises  which  represer  ts  the  three 
components  of  force  on  the  missile  caused  by  winds.  Let  the  3x3  covariance 
matrix  of  these  white  noises  be  denoted  by  Q  (t). 
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From  Theorem  1,  Section  V.2,  Equations  39  and  42,  it  is  seen  that  the  optimal 


filter  satisfies  the  equations : 


(61) 


(64  -  Ai  (t)  §x  -  A2  (t)  £>y 
6$  -  Bj  (t)  6x  -  B2  (t)  §y 


A  V 

b‘  ) 

A  A  I 

B3  (t)  6z  -  b2 ) 


The  matrix  right  of  the  equal  sign  is  indicated  in  terms  of  submatricies.  This 
notation  is  used  frequently  in  succeeding  sections. 

The  hat  symbols  are  used  to  denote  estimated  values.  The  expression  K  (t)  is 
given  by: 

K  (t)  =  £  (t)  H  (t)T  R  (t)"1  (62) 


The  matrix  £  is  the  solution  of  the  equation: 

£  =  -  £Ht  R"1  HE  +  F£  +  £Ft  +  DQDt 
Note  that  if  the  vehicle  was  at  the  position 


(63) 


(64) 
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using  the  control  law: 


and  the  wind  variations  Aofw,  and  AjSw  were  zero,  the  linearized  equa¬ 

tions  of  motion.  Equation  49,  imply  that  the  derivatives  of  the  variation.!,  6x, 
6y,  etc.,  are  given  by: 


Hence,  Equation  66,  with  the  variation  fix,  fly,  etc.  replaced  by  these  esti¬ 
mated  values  6x,  oy,  etc.,  is  the  estimated  value  of  this  derivative  based  on 
past  observations. 

The  expression 

(60  -  Ax  <t$x  -  A2  (t)$y  -  bj  \ 

A  A  A  a  <67> 

6$  -  Bi  (t)6x  -  B2  (t)fly  -  B3  (t)6z  -  b2  / 

is  a  weighting  by  the  matrix  K  (t)  of  the  difference  between  the  actual  measured 
variations  in  llne-of-sight  angles  and  estimates  for  these  quantities  based  on  past 
obrervations.  These  two  quantities  are  added  together  in  Equation  61  to  produce 
the  derivative  of  the  estimates  of  the  deviations. 
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(2)  Optimal  Filter  Based  TV  Tracker  and  Inertial  Platform.  In  Section  V.l.c, 
a  discussion  of  inertial  platlorm  measurements  is  given.  The  conclusion  of  this 
discussion  is  that  the  platform  accelerations  should  be  treated  as  exactly  meas¬ 
ured  quantities.  The  equations  for  the  missile  accelerations  maybe  written: 

x  =  ax  (t)  y  =  a2  (t)  z  =  a3  (t)  (68) 

in  which  a^t),  a2(t),  and  ag(t)  are  the  three  measured  missile  accelerations 
in  an  inertial  coordinate  system.  Because  these  measurements  are  exact,  no 
additive  white  noises  are  involved  in  them.  Thus,  they  do  not  exactly  fit  the 
theory  given  in  Section  V .  2.  To  handle  this  situation,  the  Kalman  optimal  esti¬ 
mate  of  the  positions  and  velocities  will  be  computed  in  a  different  maimer  using 
techniques  from  Reference  4.  A  similar  situation  in  a  simpler  case  is  discussed 
in  Reference  4,  Pages  78  to  82. 

Because  exact  measurements  of  the  missile  accelerations  are  being  made,  the 
expressions  for  the  accelerations  given  by  the  equations  of  motion  are  redundant. 
Thus,  measured  accelerations  of  the  platform  will  be  used  to  replace  the  equa¬ 
tions  of  motion.  The  equations  for  the  missile  accelerations  and  line-of-sight 
measurements  are  given  by: 


x  =  ax  <t) 

(69) 

y  =  *2  w 

(70) 

z  =  a3  <t)  -  g 

(71) 

=  tan'1  *•  +  bx  +  Nj 

(7->) 

_T  Z 

V  (x'z  +  ya)  2 

(73) 

The  Kalman  optimal  filter  for  this  system  will  be  computed  considering  Equa¬ 
tions  69,  70,  and  71  as  equations  of  motion  and  Equations  72  and  73  as 
measurements. 

To  apply  the  theory  of  Reference  4,  Pages  54  -  55,  formulas  11c  and  111c,  the 
equations  must  be  linear,  first  order,  homogeneous  equations.  Let  x  (t),  y  (t), 


69 


1 


l 


z  (t),  u  (t),  v  (t),  and  w  <t)  denote  the  positions  and  velocities  in  a  rectangular  ■ 

i 

coordinate  system  of  an  object  whose  Initial  position  and  velocity  ure  the  same 
as  the  initial  estimates  of  the  missile's  position  and  velocity  and  is  acted  on  by 
accelerations  a!  (t),  a2  (t),  a3  (t)  -  g.  Let  x  <t),  y  (t),  z  (t),  u  (t),  v  (t),  and 

w  (t)  denote  the  actual  positions  and  velocities  of  the  missile.  Denote  x  (t)  - 

^ 

x  (t),  y  (t)  -  y  (t),  etc.,  by  fix,  6y,  etc.  Then  the  equations 

•  •  j 

fix  =  6u  6u  =  0  bj  =  0  : 

*  • 

fiy  =  6v  6v  =  0  b2  =  o  (74) 

•  •  | 

6z  =  fiw  6w  =  0 

are  satisfied.  Equations  59  are  approximately  valid  with  fix,  fiy,  fiz,  replaced 

by  fix,  fiy,  fiz.  ' 

From  Equations  11c  and  111c  of  Reference  4,  Page  54,  the  optimal  estimates  of 

fix,  Zy,  etc. ,  satisfy  the  equations  \ 


in  which. 


F  = 


0  10  0  0 
0  0  10  0 
0  0  0  1  0 
0  0  0  0  0 
0  0  0  0  0 
0  0  0  0  0 
0  0  0  0  0 
0  0  0  0  0 


(78) 


A  _  A  A _  A 

Letting  x  (t)  =  x  (t)  +  ox  (t),  y  (t)  =  x  (t)  +  ox  (t),  the  equation  of  the  filter  may 

be  expressed  in  terms  of  x,  y,  ...  by: 


(3)  Discussion  of  Optimal  Filters.  In  the  filter  based  on  the  tracker  and  IMU, 
exact  measurements  are  used  for  the  derivatives  of  the  estimated  velocities.  In 
the  filter  based  on  the  tracker,  these  accelerations  were  estimated  using  past 
estimates  of  the  missile's  position  and  velocity  and  the  equations  of  motion.  The 
filter  for  the  tracker  and  EMU  should  give  better  estimates  in  that  there  will  not 
be  an  accumulation  of  error  from  the  estimation  of  accelerations  as  in  the  filter 
for  the  tracker  only.  In  the  filter  for  the  tracker  and  EMU,  the  effects  of  winds 
are  taken  directly  into  account  by  the  measured  missile  accelerations. 

Note  that  in  the  filter  for  the  tracker  and  IMU,  Equation  77  for  the  matrix  £ 
does  not  involve  the  covariance  matrix,  Q,  of  the  winds,  while  in  Equation  63, 
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concerning  the  filter  for  the  ASM  with  tracker  only,  the  matrix  DQDT  appears  as 
an  additive  term.  The  matrix  £  is  the  covariance  matrix  of  the  difference  be¬ 
tween  the  true  and  estimated  position  deviations.  An  extra  nonnegative  term, 
DQDt,  in  the  derivative  of  this  quantity  will  make  the  covariances  larger. 

Notice,  though,  that  a  precise  comparison  of  the  two  covariance  matrices,  L, 
cannot  be  made  on  this  basis  because  the  F  matrices  in  Equations  77  and  63  are 
different. 

In  each  of  these  filters  discussed,  the  biases  bj  and  b£  in  the  line-of-sight  angle 
measurements  are  estimated.  These  quantities  are  called  bo  resight  errors. 

In  a  conventional  tracking  seeker  the  angular  rate  of  turning  of  the  seeker  axis 
is  used  as  an  estimate  of  the  angular  rate  of  change  of  the  line  of  sight.  This 
quantity  is  used  primarily  to  eliminate  the  boresight  error  problem. 

Both  the  estimate  of  line-of-sight  angle  rate  from  tracker  axis  rate  and  the  esti¬ 
mate  of  bias  error  in  the  optimal  filter  will  have  a  transient  response.  The 
transient  response  of  the  axis  rate  estimate  will  be  due  to  the  tracker  dynamics. 
The  transient  responses  in  the  estimates  of  the  biases  are  governed  by  the  filter 
equations.  Estimating  the  biases  with  the  filter  equations  has  the  advantage  that 
the  errors  due  to  the  response  of  the  tracker  to  missile  motions  are  eliminated. 

In  instrumenting  the  filter  when  a  gimbaled  TV  tracker  is  used,  a  difficulty  is 
encountered  in  that  to  compute  4>  and  $,  the  measurements  yj  and  y2  of  Section 
V.l.b  must  be  differentiated.  The  derivatives  of  y^  and  y2  contain  white  noise. 

It  is  not  practical  to  instrument  a  device  to  obtain  such  derivatives.  This  dif¬ 
ficulty  may  be  avoided  by  using  a  technique  suggested  in  Reference  5.  This  tech¬ 
nique,  as  applied  to  the  present  situation,  follows.  The  equations  of  the  optimal 
filter  are  given  by  Equation  79. 
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in  which  and  y2  are  the  measured  tracker  angle  rates  that  are  defined  in 
Section  V.  l.b  and  K  (t)  is  the  matrix  given  by  Equation  76. 

Then,  using  Equations  79  and  80  and  eliminating  4>  -  yj/c2  and  $  -  Vc4 
by  use  of  Equations  27  and  28: 
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The  hatted  quantities  may  be  eliminated  from  the  right  side  of  Equation  81  by 
using  Equation  80.  Thus,  a  filter  for  the  starred  variables  is  obtained  whose 

inputs  are  ylt  y2,  J  y^  (s)  ds,  J  y2  (s)  ds,  r»i  (tQ)  +  bx,  and  +  b2. 

to  t0 

Thus,  derivatives  of  and  y2  are  not  involved  as  inputs,  to  this  filter.  Equation 

80  implies  that  the  desired  estimates  are  obtained  by  adding  K  (t)  ^  y  ) t° 

the  output  of  this  filter.  1 

A  set  of  typical  weighting  factors  for  the  optimal  filter  are  plotted  in 
Figures  24  and  25.  These  weighting  factors  are  the  elements  of  the  matrix  K  (t) 
defined  by  Equation  76,  and  used  in  the  optimal  filter,  Equation  79.  Examination 
of  the  right  side  of  Equation  79  shows  that  the  matrix  K  (t)  gives  the  relative 
weighting  of  the  two  bracketed  vector  quantities  on  the  right  side  of  this  equation 

in  the  derivative  of  the  position  and  velocity  estimates.  The  nonzero  elements  of 

l 

the  first  vector  in  Jthe  right-hand  side  of  Equation  79  are  the  estimated  ASM  velocity 
components  and  th^  measured  ASM  acceleration  components.  The  second  vector 
is  the  difference  between  the  estimated  and  measured  llne-of-sight  angles.  The 
matrix  K  <t)  is  the  "assignment  of  relative  importance  of  these  quantities"  in  the 
derivative  of  estimated  position  and  velocity. 

In  Figure  26  the  weighting  coefficient  of  llne-of-slght  azimuth  on  y  velocity  Is 
plotted  for  each  of  the  three  nominal  trajectories  used  to  illustrate  the  dependence 
of  a  typical  weighting  coefficient  on  the  nominal  trajectory. 

4.  Optimal  Navigation  Implementation.  To  discuss  the  implementation  of  the 
optimal  navigation  system,  consider  the  case  in  which  both  an  inertial  platform 
and  a  gimbaled  TV  tracker  are  used  as  sensors  by  the  ASM.  In  the  cases  in 
which  the  gimbaled  tracker  is  replaced  by  a  body  fixed  tracker  or  the  Inertiai 
platform  omitted,  the  equations  of  the  optimal  navigation  system  are  much  simpler. 

Equations  80  and  81  give  the  optimal  filter  and  Equation  53  the  optimal  control 
law.  These  make  up  the  optimal  navigation  system.  In  implementing  this  system, 
the  inertial  platform  acceleration  measurements  a^t),  a^it) ,  a^(t) ,  the  TV  tracker 
line-of -sight  angle  rates  y^t),  y2(t),  and  the  measurement  of  initial  tracker  angles 
T)j(t0)  +  bj,  iJ2(tQ)  +  b2,  are  inputs  to  the  differential  Equations  81.  Initial 
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(B)  WEIGHTING  COEFFICIENT  — 
Y- POSITION  ESTIMATE 


Figure  24:  OPTIMAL  FILTER  COEFFICIENTS  FOR  POSITION  ESTIMATE 
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LINE-OF-SIGHT  AZIMUTH  COEFFICIENT  UNE-OF-SIGHT  ELEVATION  COEFFICIENT 


•  Long-Range  Semiballistic  Trajectory 

•  Slant  Range  at  Acquisition  =  15,000  Ft. 

•  Tracker  Noise  =  10.0  Millirodion  (rms) 
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(A)  WEIGHTING  COEFFICIENT  — 
X-  VELOCITY  ESTIMATE 


0  2  4  6  8 

TIME- (SEC) 
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<C)  WEIGHTING  COEFFICIENT — 
Z  -  VEL0C ITY  ESTIMATE 


Figure  25:  OPTIMAL  FILTER  COEFFICIENTS  FOR  VELOCITY  ESTIMATES 
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•  Tracker  Noise  “  1  Milliradian  (rms) 


0  2  4  6  8  10 

TIME- (SEC) 

Figure  26:  EFFECT  OF  NOMINAL  TRAJECTORY 
ON  OPTIMAL  FILTER  COEFFICIENTS 


conditions  for  this  differential  equation  are  the  estimates  of  position  and  velocity 
of  the  ASM  at  target  acquisition  supplied  by  the  inertial  platform.  The  matrices 
K  (t)  and  k  (t),  which  are  coefficients  of  Equations  81  and  80,  may  be  precomputed 
from  Equations  76  and  77  and  stored  in  the  onboard  ASM  computer.  Typical 
examples  of  these  filter  coefficients  are  shown  in  Figures  24,  25,  and  26.  Equa¬ 
tions  80  and  81  would  be  solved  by  the  computer  to  obtain  estimates  of  the  mis¬ 
sile's  position  and  velocity  in  rectangular  coordinates.  These  position  and 
velocity  estimates  would  be  subtracted  from  values  of  the  same  variables  of  the 
nominal  trajectory  to  obtain  estimated  deviations  from  the  nominal  trajectory. 

The  nominal  trajectory  in  rectangular  coordinates  would  be  stored  in  the  ASM 
on  board  computer. 

The  estimated  deviations  from  the  nominal  trajectory,  in  rectangular  coordinates, 
would  be  multiplied  by  the  linearized  transformation  matrix  between  rectangular 
and  flight  path  coordinates  to  express  these  estimated  deviations  in  flight  path  co¬ 
ordinates.  The  transformation  is  given  by  the  matrix  3  that  is  used  in  Appendix 
II,  Equation  107. 

-1  T 

These  variations  are  then  multiplied  by  the  matrix  -M  G  U  of  Equation  53 
to  obtain  the  optimal  control  variations  in  angle  of  attack  and  bank  angle.  These 
variations  are  then  added  to  the  nominal  values  of  these  variables  stored  on  the 
ASM. 

In  this  computation  the  variations  in  position  and  velocity  are  multiplied  by  the 
two  by  six  product  matrix  -M  G  U  E  and  the  result  added  to  the  nominal 
bank  angle  and  angle  of  attack.  Typical  control  law  coefficients  for  this  matrix 
are  shown  in  Figures  21  through  23.  The  matrix  U  is  precomputed  from  Equa- 
tion  54  with  boundary  condition  from  Equation  55.  The  matrices  M  ,  G  ,  and 
E  are  known.  Hence,  this  product  matrix  may  be  precomputed  and  stored  on¬ 
board  the  ASM.  The  angle  of  attack  and  bank  angle  are  then  fed  as  command 
variables  to  the  autopilot. 

The  missile-borne  digital  computer  requirements  for  implementation  of  the  ASM  opti¬ 
mal  navigation  were  obtained  by  coding  Equations  53,  80,  and  81  in  a  whole  num- 


ber  digital  computer  language.  The  results  of  the  analysis  performed  on  the 
navigation  loop  were  combined  with  estimates  of  the  remaining  computational 
load  to  size  the  missile-borne  computer. 

A  timing  summary  and  an  airborne  computer  storage  requirement  for  the  various 
quantities  of  ‘he  computation  are  presented  in  Tables  V  and  VI. 

—1  T 

In  obtaining  this  table,  the  functions  K(t),  K(t),  and  -M  G  U  H were  re¬ 
presented  by  tabular  functions  of  time  with  fifty  points  for  each  variable.  The 
nominal  trajectory  was  expressed  as  a  tabular  function  of  time  with  twenty  five 
data  points.  The  estimated  2650  word  storage  capability  is  about  1500  words 
greater  than  similar  estimates  for  conventional  proportional,  and  pursuit 
guidance  data.  Any  improvement  is  not  expected  to  be  so  great  as  to  change  the 
class  of  the  airborne  digital  computer  required. 

Due  to  the  large  quantities  of  numerical  data,  a  whole  number  general  purpose 
(or  hard  wired  special  purpose)  computer  is  dictated  for  the  optimal  control 
application. 


Table  V:  STORAGE  SUMMARY 


Navigation  Words 

K  (t)  table  400 

K  (t)  table  400 

-M"1  GTU  Stable  250 

Nominal  State  and  Control  Variables  150 

Constants  38 

Variables  30 

Incremental  Inputs  5 

Instructions  581 

SUB  TOTAL  1854 

Nonnavigation  Estimates 

Resolver  Inputs  5 

Gyro  Torquing  100 

Flight  Control  Functions  50 

Check  and  Calibration  150 

Initial  Alignment  40 

Discrete  Functions  30 

Monitor  25 

Status  50 

SUB  TOTAL  450 

TOTAL  2304 

15%  uncertainty  346 

RESULTANT  TOTAL  2650 


Table  VI:  TIMING  SUMMARY 
Navigation  Instructions 

per  major  cycle  3397 

Other  instructions 

per  major  cycle  250 

TOTAL  3647 

For  a  0.2  second  major  cycle,  which  is  quite  fast,  the  average  instruction 
execution  time  must  be  less  than  54. 8  microseconds. 
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In  the  weight  category  of  10  to  30  pounds,  the  following  is  a  brief  survey  of 
some  available  computers  that  satisfy  the  requirements  of  the  missile -borne 
digital  computer. 


Manufacturer 

A.  C.  Sparkplug 

Anna 

Autonetics 

Control  Data  Corp. 

Honeywell 

Hughes 

Lear  Siegler 

Litton 

Nortronics 


Model 

Magic  series 
M 169 

D-26  series 
5300  series 
ALERT 

HCM  205  and  206 

DIVIE 

C-221 

NDC-1051 


The  following  numerical  approximations  were  made  in  coding  the  optimal  navi¬ 
gation  problem . 


1.  Tto-'V-V-  (f) 


2.  Given  y  as  an  estimate  of  y  =  Vz"|  the  second  order  recursion  formula 
°  „  .  2 


yi  +  1  =  yi  ‘ 


\  (yt  -z) 

3y£2  +  z 


is  assumed  to  be  adequate  after  two  iterations.  Yq  is  given  for  the  first 
At  and  thereafter  the  starting  yQ  is  obtained  by  using  the  result  of  the 
previous  guidance  cycle. 


3.  Second  order  interpolation  is  required  to  extract  values  of  the  tabular 
functions  from  their  tables . 
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Second  order  interpolation  was  used  for  the  table  look-up  functions.  Improve¬ 
ments  in  computer  requirements  could  result  if  data  was  supplied  as  segments 
of  polynomials  pieced  together  to  provide  an  adequate  fit  to  the  numerical  data. 

a.  Conclusions  from  the  Optimization  Computation.  Tables  III  and  IV  of  Section 
IV  compare  the  standard  deviation  of  miss  distance  for  the  optimal  system  to  the 
standard  deviation  of  miss  distance  for  pursuit  and  proportional  guidance  for 
several  nominal  trajectories  and  sensor  noise  levels.  Appendix  n  contains  the 
method  used  to  make  the  calculations  for  these  tables. 

In  Appendix  II,  it  is  shown  that  the  miss  distance  can  be  resolved  into  a  term 
that  is  due  to  the  error  in  estimating  the  ASM  position  deviation  and  a  term  that 
depends  on  the  controller.  These  are  called  filter  error  and  controller  error, 
respectively.  Examples  of  the  time  history  of  these  errors  are  given  in  Figure 
27.  The  square  root  of  the  sum  of  the  squares  of  these  quantities  at  the  final 
time  gives  the  miss  distance.  In  all  of  the  calculations  carried  out,  the  contri¬ 
bution  of  the  controller  error  to  miss  distance  was  negligible.  This  shows  that 
almost  all  the  miss  distance  is  due  to  error  in  estimating  the  ASM  position 
deviation. 

The  comparison  with  proportional  and  pursuit  guidance,  given  in  Table  IV,  shows 
that  these  guidance  laws  have  almost  the  same  miss  distances  as  the  optimal 
system  for  a  variety  of  sensor  noise  levels.  This  shows  that  these  simple  sys¬ 
tems  are  using  the  Bensor  information  and  controlling  in  a  nearly  optimal  manner. 

In  Appendix  n,  it  is  shown  that  the  filter  error  is  completely  independent  of  the 
controller,  and  that,  if  a  linear  feedback  controller  of  the  form 
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1 0  MISS  DISTANCE 


Long-Range  Semiballistic  Trajectory 
Slant  Range  at  Acquisition  =  15,000  Ft. 
Tracker  Noise  =  50  Milliradians 


0  2  4  6  8 

TIME  -  (SEC) 


Figure  27 :  COMPARISON  OF  CONTROL  AND  FILTER  ERROR 
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was  used  with  the  optimal  filter,  a  modification  of  Equation  110  governing  con¬ 
troller  error  still  holds.  It  can  be  shown  that,  if  the  controller  A(t)  is  chosen, 
so  that  the  system 


drives  the  component  of  position  deviation  normal  to  the  nominal  trajectory  to 
zero,  that  controller  error  will  be  driven  to  zero.  There  are  many  choices  of 
A(t)  that  will  fulfill  this  requirement.  Therefore,  there  are  many  linear  feed¬ 
back  controllers  that  will  behave  in  a  nearly  optimal  fashion. 

In  the  computation  of  the  optimum  filter  coefficients,  the  matrix  differential 
Equation  63  for  the  matrix  I  must  be  solved.  This  matrix  is  the  covariance 
matrix  of  the  difference  of  the  actual  state  vector  of  the  ASM  and  the  value 
estimated  for  this  vector  by  the  filter.  Examination  of  the  elements  of  this 
matrix  at  various  times  will  show  how  rapidly  the  optimal  filter  is  estimating 
corresponding  components  of  the  state  vector.  This  is  illustrated  in  Figures 
28  and  29. 

In  Figure  28  the  standard  deviations  of  the  differences  between  the  true  values 
and  the  estimated  values  of  the  variations  in  x  and  z  are  plotted.  Note  that 
these  standard  deviations  do  not  substantially  decrease  with  time.  In  Figure  29 
these  same  values  are  stressed  In  a  different  coordinate  system .  In  this  figure 
the  standard  deviations  of  the  same  differences  resolved  into  components  par¬ 
allel  to  and  perpendicular  to  the  nominal  trajectory  are  plotted.  Note  that  the 
component  perpendicular  to  the  Dominal  trajectory  decreases  very  rapidly  to 
zero,  while  the  component  parallel  to  the  nominal  trajectory  decreases  very 
slowly. 
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This  may  be  explained  by  noting  that  only  line-of-sight  angles  and  ASM  acceler¬ 
ations  are  being  measured.  Therefore,  very  little  Information  about  the  missile's 
position  along  the  nominal  trajectory  is  being  supplied.  These  measurements  do 
supply  information  about  the  missile's  position  normal  to  the  nominal  trajectory. 
Fortunately,  as  was  pointed  out  in  Section  V.3.b,  the  deviation  in  position  normal 
to  the  nominal  trajectory  is  the  imporlant  component  in  computing  miss  distance. 
This  question  is  explored  in  further  detail  in  Appendix  m. 


•  •  ing-Range  Semi  ballistic  Trajectory 

•  Slant  Range  at  Acquisition  =  15,000  Ft. 
•Tracker  Noise  =  1.0  Millirodian 


O 

s 

< 

9  400 


X  Deviation 


§  600 


Z  Deviation 


0  2  4  6  8 

TIME  ~  (SEC) 

Figure  28:  STANDARD  DEVIATIONS  OF  X  AND  Z 
POSITION  ESTIMATES 


I 


STANDARD  DEVIATION~(FT.) 


#  Long-Range  Semi  ballistic  Trajector/ 

#  Slant  Range  at  Acquisition  =  15000  Ft. 

#  Tracker  Noise  =  1.0  Milliradianj(rms) 


Figure  29:  COMPARISON  OF  THE  ACCURACY  OF  ESTIMATION  OF  POSITION 
DEVIATION,  PARALLEL  TO  AND  PERPENDICULAR  TO  THE 
NOMINAL  TRAJECTORY 
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SECTION  VI 

NORMAL-ACCELERATION  AUTOPILOT  STUDIES 

A  normal-acceleration  commanded  autopilot  was  selected  for  the  ASM  because 
normal  acceleration  is  a  direct  command  parameter  for  flight -path  corrections. 
Because  of  this  feature,  it  leads  to  simpler  guidance  laws  for  the  near-impact 
phase  than  attitude  control.  Thus,  the  practicality  of  instituting  optimal  control 
techniques  in  the  near-impact  phase  is  increased  through  simplification  of  guid¬ 
ance  law  expressions  using  this  control  parameter.  Normal  acceleration  control 
also  provides  a  parameter  for  limiting  commanded  maneuvers  to  stay  within  the 
structural  limitations  of  the  missile. 

The  normal-acceleration  autopilot  was  investigated  for  launch,  midcourse,  and 
near- impact  (flight  from  target  acquisition  to  impact)  phases  of  three  nominal 
trajectories — long-range  aemiballistic,  short-range  semiballistic,  and  low- 
altitude  skip.  Dynamic  pressures  from  10  to  5500  lb/ft2  were  considered.  The 
Mach  number  range  was  0. 59  to  4.2.  Seven  representative  flight  conditions 
(Figure  30)  were  investigated  in  fixed  point  studies.  Because  the  major  interest 
of  the  program  was  in  the  homing  phase  of  the  ASM,  emphasis  was  placed  upon 
selecting  autopilot  gains  and  compensation  to  obtain  satisfactory  control  system 
response  during  the  terminal  phases  of  nominal  trajectories  (approximated  by 
Flight  Conditions  1,  2,  and  3). 

Stability  and  time  response  criteria  for  the  autopilot  were: 

(1)  The  rigid  mode  must  exhibit,  as  minimal  stability  characteristics,  a  6-db  gain 
margin  and  a  30-degree  phase  margin; 

(2)  All  body-bending  modes  must  be  at  least  6-db  gain  stable  (to  account  for 
uncertainty  in  the  phase  of  the  servo  response  at  the  body-bending  frequencies); 

(3)  The  controlled  missile  must  have  a  time  response  as  fast  as  a  first-order  sys¬ 
tem  with  a  time  constant  of  0. 5  second. 

The  time  response  criteria  were  based  on  the  results  of  the  preliminary  near-impact 
phase  digital  simulations  that  indicated  that  a  0.5  second  constant  made  a  negligible 
contribution  to  impact  error.  This  preliminary  simulation  included  position  offset 


89 


errors  and  wind  shear  effects.  As  was  indicated  in  the  discussion  of  autopilot 
response  in  Section  IV,  this  criterion  was  not  adequate  when  the  more  complete 
simulation  was  used.  No  criterion  was  set  for  acceptable  steady-state  errors 
in  response  to  commands.  The  nominal  autopilot  design  was  a  Type  1  controller, 
and  therefore  exhibited  zero  steady-state  error  for  step  inputs. 

The  autopilot  study  was  conducted  in  two  phases.  The  first  phase  consisted  of 
the  development,  through  analytical  methods  (root  locus  and  frequency  response 
analyses) ,  of  a  programmed  gain  nominal  autopilot  that  could  be  used  in  the  guid¬ 
ance  evaluation  studies.  This  phase  resulted  in  the  Type  1  autopilot  with  com¬ 
pensation  selected  to  give  adequate  response  during  the  near-impact  phase  of 
flight.  Gains  were  developed  to  meet  the  stability  criteria  ior  all  flight  conditions 
(except  Condition  7,  a  condition  of  very  low  dynamic  pressure,  which  was  unstable). 
The  synthesis  of  a  normal  acceleration  autopilot  presented  no  unusual  difficulties 
for  the  configuration  studied.  However,  care  was  required  in  selecting  autopilot 
gains  that  provided  acceptable  time  response  and  steady  error  performance  and 
did  not  violate  the  stability  requirements  when  body-bending  and  tail-*.ags-dog 
effects  were  considered. 

A  somewhat  unconventional  root  locus  technique  was  used  and  aided  in  arriving 
at  a  satisfactory  compromise.  This  approach  is  presented  in  the  discussion  that 
follows. 

The  second  phase  consisted  of  analog  computer  verification  of  the  analytical  studies 
of  the  nominal  Type  1  autopilot  design,  and  the  investigation  of  three  nonlinear 
effects  on  the  controlled  missile  performance.  The  nonlinear  quantities  studied 
were:  (1)  control  servomechanism  rate  limiting;  (2)  pure  transport  lag  in  the  servo¬ 
mechanism;  and  (3)  variations  in  control  fin  effectiveness  due  to  "masking"  effects 
at  angle  of  attack. 

Additional  analytical  and  computer  Investigations  were  made  of  an  advanced  con¬ 
troller  (a  quasladaptive  concept)  that  eliminates  the  need  for  gain  changing  or 
programming  during  the  flight.  This  controller  concept  was  based  on  inhouse 
extensions  of  optimal  bistable  controller  investigations  initiated  by  Gieseking 
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(Reference  7)  and  others.  This  study  resulted  in  an  autopilot  which  combined 
a  bistable  controller  and  a  Type  0  (no  integral  compensation)  autopilot.  The 
need  for  gain  and  bistable  controller  coefficient  changes  was  eliminated,  and 
satisfactory  response  and  steady-state  error  performance  was  obtained  for  all 
six  of  the  major  flight  conditions.  A  discussion  is  also  presented  on  other  advanced 
control  concepts  that  appear  applicable  to  the  ASM  autopilot. 

1.  Nominal  Autopilot  Design.  Tor  all  autopilot  studies,  the  ASM  configuration 
representation  included  two  flexible  body-bending  modes,  a  first  order  fin  servo, 
the  tail-wags-dog  effect  of  fin  inertia,  and  a  linearized  representation  of  the  rigid 
body  dynamics  for  Flight  Conditions  1  through  6,  which  were  shown  in 
Figure  30. 

The  resulting  autopilot  design  is  shown  in  Figure  31.  Feedback  signals  for  the 
functional  autopilot  were  taken  to  be  linear  body  normal  acceleration  and  body 
angular  rate.  The  location  of  the  normal  accelerometer,  dictated  by  the  constraint 
of  available  space  in  the  missile's  internal  arrangement,  was  about  5.3  feet  ahead 
of  the  vehicle's  center  of  gravity.  This  instrument  location  yielded  a  feedback 
signal  having  more  than  a  desired  content  of  vehicle  angular  acceleration  and  body¬ 
bending  signal  i—apiltude.  However,  accelerometer  placement  was  near  the  sec¬ 
ond  body-bendbig  node,  so  it  did  not  sense  second  body  bending  motion.  The  at¬ 
titude  rate  sensor  was  located  in  the  autopilot  electronics  section,  near  the  anti¬ 
node  of  the  first  body  bending  mode,  and  hence  sensed  the  flexible  body  modal 
motions  of  only  the  second  and  higher  modes.  A  hydraulic  fin  actuator  was  as¬ 
sumed,  having  a  first  order  response  with  a  60  rad/sec  characteristic  frequency 
and  (for  the  final  analog  computer  studies)  a  deflection  rate  limit  of  150  degrees/ 
second. 

The  compensation  used  for  body  bending  stabilization  consists  of  two  sets  of  second- 
order  filters  employing  complex  zeros  in  conjunction  with  real-axis  poles.  This 
form  of  body  bending  compensation  has  a  less  adverse  effect  on  rigid  mode  re¬ 
sponse  than  simple  lag  compensation,  and  can  be  implemented  using  RC  networks 
and  a  single  operational  amplifier  for  each  of  the  quadratic  filters.  The  zeros 
of  this  quadratic  compensation  are  placed  at  the  middle  of  the  range  of  variation 
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0*01 


DYNAMIC  PRESSURE  (LB/FT2) 

Figure  30:  ASM  FLIGHT  CONDITIONS 


of  first  body  bending  frequency  from  launch  to  burnout,  so  that  the  compensation 
will  ensure  first  body  bending  stability  throughout  the  flight.  The  compensation 
network  parameters  are  fixed  at  values  that  permit  satisfactory  controlled  sys¬ 
tem  response  for  all  flight  conditions  in  the  terminal  homing  phases  of  the  three 
reference  trajectories.  The  required  values  for  Kj/Kg,  Kg,  and  Kg  are  shown 
in  Table  VII  for  Flight  Conditions  1  through  6. 


Table  VII:  REQUIRED  AUTOPILOT  GAINS 


Flight  Condition 

Kj/Kj 

JSL 

1 

3 

-0.0274 

-0.042 

2 

3 

-0.0371 

-0. 050 

3 

3 

-0.0525 

-*  '56 

4 

3 

-0.0525 

.  \J 

5 

3 

-0.  0909 

-O.tl 

6 

3 

-0. u909 

-0.08 

A  somewhat  unconventional  procedure  was  used  for  selecting  satisfactory  auto¬ 
pilot  gains  (Kj,  Kg,  and  K^).  The  procedure  used  was  to  select  values  of  Kj  and 
Kg  based  on  the  roots  of  the  numerator  of  the  6(/6  transfer  function.  It  can  be 
shown  that  the  roots  of  the  numerator  become  zeros  in  the  final  loop  closure  for 
z/zc  using  the  Kg  gain.  Placing  these  zeros  in  a  well  damped  position  assists  in 
meeting  the  overall  system  stability  and  response  requirements  with  the  final  loop 
closure. 

The  intermediate  missile  dynamics  transfer  functions,  0/5  and  2/6,  r  ’  the  appro¬ 
priate  pole  and  zero  values  as  a  function  of  flight  condition  are  sho\  •<  in  V  ible  VO. 
Note  that,  because  of  sensor  location,  Zg  ^  and  P^  g  cancel  in  0/6  and,  similarly 
in  z/6,  Zlg  19,  and  P50  cancel.  The  Inclusion  of  tail-wags-dog  (fin  inertia)  ef¬ 
fects  in  the  vehicle  dynamics  transfer  functions  results  in  higher  order  numerators 
than  denominators.  When  the  remainder  of  the  servo  and  sensor  dynamics  and 
loop  closures,  as  shown  in  Figure  31,  are  completed,  final  transfer  functions  have 
the  expected  lower  order  numerators.  Root  loci  with  Kg  as  the  gain  parameter 
were  first  obtained  for  the  numerator  of  the  0^/6  transfer  function  for  a  wide  range 
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Table  VHI:  MISSILE  DYNAMICS  TRANSFER  FUNCTIONS 
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of  the  ratio  Kj/K|  at  each  flight  condition.  (The  composition  of  this  transfer  func¬ 
tion  can  be  determined  from  Figure  31.)  Figure  32  shows  typical  loci  of  system 
zeros  (or  roots  of  the  (5f/6  numerator)  for  Flight  Condition  3  with  values  of 
Kj/Kg  =  2,3,  and  4  respectively.  Because  the  roots  on  thesS  loci  become  zeros 
of  the  overall  system,  it  is  desirable  to  place  them  in  a  well  damped  position. 

The  real  open- loop  pole,  located  at  -2,  -3,  and  -4  respectively  in  these  figures, 
is  given  directly  by  the  ratio  Kj  /Kg*,  and  indicates  that  system  response  becomes 
faster  for  higher  values  of  Kj/Kg  .  However,  as  shown  in  these  figures,  when 
Kj/Kg  is  greater  than  3,  the  complex  branch  of  the  locus  fails  to  progress  in  an 
arc  down  to  the  real  axis,  and  these  roots  are  left  in  a  position  representing  less 
damping.  This  results  in  a  lighter  damped  rigid  mode,  and  therefore  is  undesir¬ 
able.  These  considerations  lead  to  the  choice  Kj/Kg  =  3. 

Figures  33  through  37  show  the  roots  of  the  0?/6  numerator  for  the  remainder  of 
flight  conditions,  except  Condition  7.  These  root  loci  are  based  on  the  ratio 
Kj/  Kg  =  3,  and  indicate  satisfactory  values  for  Kg.  From  these  plots,  values 
for  Kg  were  selected  that  place  the  quadratic  pair  in  a  well-damped  position. 

The  total  system  root  loci,  z/zc,  for  determining  K g  are  shown  in  Figures  38 
through  43.  The  zeros  appearing  in  these  plots  are  those  determined  by  the  loci 
of  Figures  32b  and  33  through  37,  plus  the  zeros  introduced  by  compensation  in 
the  inner  forward  loop.  At  the  left  of  each  of  these  plots  is  an  enlargement  of 
the  region  near  the  first  body  bending  pole,  and  on  the  right  is  an  enlargement 
of  the  region  about  the  origin  showing  the  rigid  mode  details.  The  gain-stable 
requirement  for  the  first  body  bending  mode  implies  that  the  operating  point  on 
the  locus  emanating  from  the  first  body  bending  mode  pole  be  at  least  6-db  below 
the  gain  at  which  the  locus  crosses  the  circle.  Uncertainty  in  phase  of  the  servo 
response  at  this  frequency  implies  an  uncertainty  in  the  direction  of  the  locus 
emanating  from  the  first  body  bending  pole.  Examination  of  these  figures  shows 
that  the  critical  design  requirements  arise  from  the  high  Mach  number  condition 
(Figure  41)  where  system  response  requirements  require  high  gains.  Limitations 
on  high  gains  arise  from  the  locus  that  crosses  the  imaginary  axis  in  the  vicinity 
of  140  rad/sec,  and  from  the  gain-stable  requirement  for  the  first  body  bending 
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Figure  34:  LOCI  OF  SYSTEM  ZEROS  —RIGHT  CONDITION  2 
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Figure  41:  ROOT  LOCI  FOR  z/z— FLIGHT  CONDITION  4 
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mode.  For  the  high  Mach  number  condition,  the  requirement  for  high  gain  arises 
in  providing  adequate  damping  of  the  rigid  mode  response,  rather  than  from  the 
requirement  for  response  time.  The  high  dynamic  pressure  condition  (Figure  38) 
does  not  impose  any  critical  design  requirements;  adequate  system  response  can 
be  achieved  with  gains  that  are  not  critical  from  a  stability  standpoint. 

Amplitude  versus  phase  plots  of  the  system  are  shown  in  Figures  44  through  49. 
The  stability  criteria  for  gain  and  phase  margins  are  shown  in  a  combined  form, 
or  "vector  margin"  requirement,  represented  by  the  ellipse.  The  frequency  locus 
in  the  vicinity  of  the  first  body  bending  mode  frequency  must  lie  6-db  or  more 
below  the  0-db  point  for  the  operating  gain  because  its  location  laterally  on  the 
plot  is  uncertain  due  to  phase  response  uncertainty  of  the  fin  servos  at  this  fre¬ 
quency.  These  plots  are  useful  in  determining  system  gains  that  satisfy  the  sta¬ 
bility  criteria  because  a  change  in  operating  gain  in  effect  shifts  the  0-db  re¬ 
ference  point  in  these  plots  and  does  not  alter  the  locus  shape.  In  general,  the 
chosen  operating  gain  is  different  than  the  gain  for  which  the  plot  was  originally 
made.  The  location  of  the  0-db  line,  corresponding  to  the  chosen  operating  gain, 
is  indicated  by  the  position  of  the  horizontal  axis  of  the  vector  margin  ellipse. 

Figure  50  shows  step  response  characteristics  for  the  chosen  operating  gains  as 
obtained  using  a  digital  computer  time  response  program.  A  first-order  0.5- 
second  time  constant  step  response  is  shown  for  reference.  The  use  of  integral 
plus  proportional  compensation  in  the  outer  forward  loop  yields  zero  steady 
state  error.  The  initially  chosen  time  response  requirement  of  0.5  second  is 
achieved  or  exceeded  for  all  cases.  The  high  frequency  oscillatory  response 
occurring  early  in  the  response  time  is  associated  with  the  residues  of  the  roots 
on  the  locus  that  cross  the  Imaginary  axis  in  the  vicinity  of  140  rad/sec.  These 
responses  were  obtained  without  rate  saturation  in  the  fin  servo.  (When  a  rate 
limit  of  150  degrees  per  second  was  added  in  the  analog  studies,  this  oscillatory 
response  was  essentially  eliminated.)  Stability  limitations  on  upper  gain  margin 
for  the  high  Maoh  number  condition  restrict  the  damping  achievable  for  the  rigid 
mode  response;  consequently,  the  overshoot  characteristic  is  high  compared  to 
the  response  achievable  for  flight  conditions  where  the  missile  is  statically  stable. 
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Figure  46  :  FREQUENCY  RESPONSE  -  FLIGHT  CONDITION  3 


Figure  47  :  FREQUENCY  RESPONSE  -  FLIGHT  CONDITION  4 
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Figure  48:  FREQUENCY  RESPONSE  -  FLIGHT  CONDITION  5 


Figure  49:  FREQUENCY  RESPONSE  -  FLIGHT  CONDITION  6 
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time  response  for  nominal  autopilot 


The  six  major  flight  conditions  were  simulated  on  art  analog  computer  to  check 
against  the  digital  gain  margin  data  and  time-response  data  and  to  study  several 
nonlinear  hardware  characteristics.  Gain  margins  and  time  responses  were 
verified. 

Fin  control  servomechanism  rate  limiting  of  160  degrees /second  had  no  degrad¬ 
ing  effect  on  autopilot  performance  and  served  to  nearly  eliminate  the  high- 
frequency  oscillatory  response  that  was  shown  in  the  digital  computer  time  re¬ 
sponses.  The  difference  in  response  Is  shown  for  Flight  Condition  5  in  Figure 
5L  The  effect  of  a  pure  time  delay  (such  as  could  result  from  a  digital  auto¬ 
pilot  and/or  the  fin  control  servo)  on  autopilot  stability  was  investigated.  It  was 
found  that  a  0. 01-second  delay  could  be  tolerated  for  Flight  Conditions  1  through 
5,  and  0. 015-second  could  be  tolerated  for  Flight  Condition  6. 

The  effects  of  aerodynamic  fin  masking  and  high  fin  angles  of  attack  on  autopilot 
stability  were  found  to  be  small  for  this  configuration. 

The  low  dynamic  pressure  (10  psf)  apogee  condition  (Flight  Condition  7)  pre¬ 
sented  special  design  problems.  Inspection  of  the  missile  acceleration  transfer 
function  (shown  in  Table  vm)  for  this  condition  shows  that  the  "accelerometer 
zeros, "  —  those  rigid  mode  zeros  whose  location  is  determined  by  the  acceler¬ 
ometer  position  along  the  missile  longitudinal  axis  —  lie  on  the  real  axis  in  the 
left  and  right  half  planes .  This  situation  prevents  the  design  of  a  stable  Type  1 
normal  acceleration  system  of  the  type  considered.  These  zero  locations  are, 
in  conditions  of  extremely  low  dynamic  pressure,  dependent  not  only  on  the  ac¬ 
celerometer  location,  but  also  on  the  tail-wags-dog  characteristic.  If  the  effec¬ 
tive  tail-wagB-dog  frequency  can  be  made  sufficiently  high,  the  accelerometer 
zeros  can  be  made  to  lie  on  the  imaginary  axis  as  they  do  for  the  other  flight 
conditions  considered  in  the  terminal  phase  of  flight.  The  tail-wags-dog  fre¬ 
quency  can  be  influenced  by  mass  balancing  of  the  control  fins,  and  this  influences 
the  minimal  dynamic  pressure  for  which  a  satisfactory  normal  acceleration  auto¬ 
pilot  can  be  designed.  Study  of  Condition  7  was  dropped  during  the  nominal  auto¬ 
pilot  design  because  the  primary  interest  was  in  the  terminal  homing  phase  of 
flight. 
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It  is  significant  to  note  the  importance  that  booty  bending  plays  In  the  autopilot 
design.  The  missile  is  relatively  stiff  structurally  with  a  ratio  of  rigid  mode  to 
first  body  bending  mode  frequencies  of  approximately  1  to  50.  Yet  when  practi¬ 
cal  consideration  is  made  of  the  phase  uncertainty  of  high  frequency  fin  servo 
response,  along  with  the  choice  of  integral  plus  proportional  control  to  produce 
zero  steady  state  error,  the  consideration  of  body  bending  directly  influences 
the  achievable  system  response  and  involves  the  designer  in  trades  on  the  choice 
of  body  bending  compensation.  If  integral  plus  proportional  control  is  not  used, 
relatively  fast  rigid  mode  response  can  be  achieved.  However,  the  design  of 
such  a  system  to  keep  steady  state  error  bounded  to  less  than  45%  requires  the 
use  of  high  gains;  these  high  gains  conflict  with  requirements  for  body  bending 
stability.  These  considerations  lead  to  the  conclusion  that  control  analyses  that 
neglect  the  effects  of  body  bending  will  very  likely  reach  invalid  conclusions  re¬ 
lative  to  the  achievable  time  response  and  steady  state  error  characteristics. 

2.  Bistable  Controller.  The  nominal  Type  1  autopilot  required  gain  changes 
over  the  flight  regime  to  provide  satisfactory  performance.  As  discussed  in 
Section  III,  this  was  accomplished  for  the  terminal  (or  homing)  phase  in  the 
analog  simulation  by  programming  Kg  as  a  function  of  q.  It  would  be  desirable 
from  a  simplification  standpoint  to  avoid  gain  programming.  A  substantial  in- 
house  research  program  on  the  application  of  advanced  control  theory  to  auto¬ 
pilots  for  defense  missiles  was  in  progress  at  the  time  of  the  ASM  investigation. 
Extensions  of  the  bistable  controller  concept  described  by  Gieseking  (Reference 
7),  and  employing  Lyapunov's  second  method  were  found  to  be  quite  successful 
in  providing  constant  gain  autopilots  with  nearly  invariant  response  for  defense 
missile  configurations.  (This  work  is  documented  in  Referer.ce  8.)  The  bistable 
concepts  provide  a  state  variable  dependent  bias  command  signal  to  a  bistable 
control  element  in  the  autopilot  forward  loop.  Because  of  the  success  of  the  de¬ 
fense  missile  application,  it  appeared  possible  that  these  same  concepts  might 
have  the  potential  for  simplifying  the  ASM  autopilot  even  though  the  ASM  is  a  dif¬ 
ferent  type  of  missile.  A  block  diagram  of  a  bistable  augmented  autopilot  resulting 
from  the  application  of  these  concepts  is  shown  in  Figure  52.  The  following  is  a 
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discussion  of  the  synthesis,  simulation,  and  performance  of  this  bistable  controller 
concept  for  the  ASM  configuration  used  in  this  program. 

A  bistable  controller  was  first  applied  to  the  nominal  Type  1  autopilot.  Results 

with  the  Type  1  autopilot  were  not  entirely  successful.  The  bistable  controller 

did  speed  up  time  responses  to  step  commands  and  improved  the  stability  -in 

Flight  Condition  4.  However,  the  system  suffered  from  relaxation-type  limit- 

♦ 

cycle  problems  and  was  quite  sensitive  to  the  6  gain  entering  the  bistable  com¬ 
mand  channel.  Further  investigations  were  conducted  with  a  Type  0  system 
(no  forward  loop  integrator).  With  this  system,  a  bistable  controller  concept 
was  synthesized  which  exhibited  satisfactory  stability  and  time  response,  and  had 
less  than  15%  steady  state  error.  The  final  controller  concept  deviated  con¬ 
siderably  from  Gieseking's  concept;  however,  his  work  did  provide  considerable 
insight  into  what  was  required  to  synthesize  the  final  bistable  controller  concept. 

The  Type  0  autopilot  was  first  studied  without  the  bistable  controller  for  perform¬ 
ance  comparison  purposes.  The  Type  0  autopilot  is  identical  to  the  Type  1  auto¬ 
pilot  except  the  forward  loop  integration  and  compensation  networks  shown  in 
Figure  31  are  removed.  Gain  limits  for  this  system  were  approximately  the 
same  as  those  found  for  the  Type  1  system.  The  Type  0  system  met  the  response 
requirement  but  required  about  the  same  gain  scheduling  as  the  Type  1  system. 

It  also  exhibited  about  10%  steady  state  acceleration  error  with  the  best  gains 
fo-  '-‘ach  condition.  With  the  autopilot  gains  set  for  Flight  Condition  1  (a  terminal 
o  ming  phase  condition)  and  held  constant,  the  system  response  was  very  poor 
in  some  flight  conditions  and  the  steady  state  acceleration  error  was  as  large 
as  45%. 

A  bistable  controller  was  then  added  to  the  Type  0  system  using  the  approach 
of  Gieseking.  As  shown  in  Figure  52,  the  bistable  contnjller  adds  a  positive  or 
negative  acceleration  command  (U)  in  parallel  with  zf  .  The  sign  of  (U)  is  chosen 
by  tho  sum  of  the  state  variables,  which  are:  acceleration  error,  zf;  attitude 
rate,  0;  and  fin  angle,  6.  The  weighting  gains  of  these  summed  state  variables 
were  determined  using  the  second  method  of  Lyapunov.  The  bistable  controller 
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shown  in  Figure  52  worked  well  and  met  response  time  requirements.  It  reduced 
steady  state  acceleration  error  to  approximately  oi  e-third  of  that  associated  with 
a  fixed  gain  Type  0  autopilot. 

Figure  53  illustrates  the  time  response  of  the  bistable  augmented  Type  0  system 
and  compares  this  response  to  that  of  an  unmodified  Type  0  autopilot  without 
gain  adjustment  to  a  Type  1  autopilot  that  has  ideally  adjusted  gains. 

The  final  bistable  controller  represents  a  deviation  from  that  described  by 
Gieseking.  The  low  pass  filter  was  added  to  the  bistable  output  to  avoid  exciting 
body  bending  modes.  Without  the  filter,  the  bistable  gain  can  not  be  made  high 
enough  to  improve  the  ASM  response  at  high  q  conditions.  The  addition  of  the 
filter  induced  a  limit  cycle.  Two  more  modifications  were  necessary  to  elimi¬ 
nate  this  limit  cycle.  The  bistable  output  magnitude  was  made  proportional  to 
ze  while  the  sign  was  controlled  by  the  sum  of  state  variable  parameters  z(  and 
9  each  weighed  by  a  constant  multiplier  (see  Figure  52).  The  gain  of  this  0 
state  variable  signal  was  about  half  the  magnitude  and  the  negative  of  what 
Lyapunov's  second  method  showed  it  should  be.  The  use  of  the  (5  state  variable 
in  this  case  was  to  stabilize  a  limit  cycle  and  not  to  improve  system  response 
time  as  in  the  basic  bistable  controller  design. 

A  6-db  tolerance  was  demonst rated  on  parameters  Bgg,  Tgg,  and  c$  with  this 
controller.  A  much  faster  response  time  could  have  been  obtained  by  reducing 
Tbs  but,  in  this  case,  the  allowable  tolerance  on  was  very  small.  Work  with 
this  control  system  indicated  that  several  bistable  controllers  performing  differ¬ 
ing  tasks  would  further  improve  performance.  One  such  controller  would  be  used 
to  improve  response  time  and  another  to  stabilize  limit  cycle  oscillations 
with  perhaps  a  third  controller  to  reduce  steady  state  errors.  Further  study 
may  indicate  which  combination  of  state  variables  is  best  for  switching  and  which 
is  beat  for  absolute  values  of  (U). 

The  haidware  necessary  to  implement  the  bistable  augmented  autopilot  is  con¬ 
ventional  and  inflight  gain  adjustment  is  not  required.  The  absolute  value  functions 
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and  simple  summing  of  signals  to  control  the  bistable  output  sign  can  be  imple¬ 
mented  using  relays  or  solid  state  switches  and  amplifiers,  or  digital  components. 

Even  though  the  final  bistable  controller  is  not  of  the  same  form  as  that  studied 
by  Gieseking,  it  was  developed  only  through  the  insight  gained  through  the  in¬ 
vestigation  of  his  technique.  The  synthesis  of  the  controller  required  a  consider¬ 
able  amount  of  cut  and  try  effort;  however,  analytical  techniques  are  being  de¬ 
veloped  which  will  make  the  synthesis  more  efficient  (Reference  8). 

3.  Self-Adaptive  Autopilots.  Several  adaptive  control  concepts  were  examined 
briefly  for  potential  application  to  the  ASM  considered  in  the  study.  The  exami¬ 
nation  was  conducted  in  light  of  known  characteristics  of  the  ASM  configuration 
(e.  g. ,  body  bending,  aerodynamic  stability,  tail-wgs-dog  effects,  and  required 
autopilot  gains)  to  outline  possible  advantages  and  potential  problem  areas  asso¬ 
ciated  with  implementing  these  systems.  The  adaptive  concepts  which  were 
considered  as  likely  candidates  were: 

•  Honeywell  high  gain  system  (Reference  9) 

•  G.  E.  adaptive  system  (Reference  10) 

•  Pitch,  yaw,  or  roll  axis  dither  to  determine  control  effectiveness  (Reference  11) 

•  Calculation  of  Ma  during  flight  (References  12  and  13) 

The  Honeywell  system  tends  to  hold  the  system  gain  to  a  value  near  the  upper 
gain  stability  limit.  TheG.  E.  system  can  have  a  preselected  gain  margin. 

Both  concepts  utilize  a  pair  of  system  poles  which  become  unstable  at  gains 
lower  than  critical  vehicle  modes  such  as  body  bending.  These  poles  may 
be  inherent  in  the  basic  flight  control  system  design  or  may  be  introduced 
especially  to  provide  a  dynamic  characteristic  for  the  adaptation  loop 
design. 

In  the  Honeywell  system,  a  limit  cycle  motion  results  when  a  pair  of  poles  are 
driven  across  the  j  w  axis;  autopilot  gain  is  adjusted  to  maintain  a  preselected 
limit  cycle  amplitude.  The  G.  E.  system  does  not  permit  a  limit  cycle  of  the 
compensation  pole,  and  consequently,  for  the  same  compensation  poles,  would 
tend  to  operate  at  a  lower  gain  than  the  Honeywell  system. 


The  difficulty  with  both  of  these  concepts  lies  in  the  selection  of  satisfactory 
poles,  invariant  with  flight  conditions  both  in  frequency  and  the  gain  at  which 
they  become  unstable,  for  monitoring  the  stability  of  the  control  system.  For 
some  vehicles,  a  critioal  vehicle  mode  may  place  a  lower  stability  limit  on  the 
autopilot.  This  problem  was  encountered  in  X-20  studies  and  may  place  very 
close  bounds  on  the  acceptable  autopilot  gains.  For  the  study  vehicle,  which  has 
relatively  high  body  bending  frequencies,  the  complex  rate  gyro  poles  could  be 
considered  suitable  for  monitoring  the  stability  of  the  autopilot  for  either  of  the 
approaches.  (The  gyro  poles  would  have  to  be  sufficiently  invariant  with  temper¬ 
ature  to  allow  their  use.)  For  Flight  Conditions  1  through  6,  this  pair  of  poles 
(at  140  rad/sec)  become  unstable  prior  to  the  first  body  bending  mode  (at  370 
rad/sec).  However,  at  Condition  7,  as  was  previously  discussed,  the  body 
bending  poles  become  unstable  first.  To  make  either  concept  work  on  the  study 
vehicle,  It  would  be  necessary  to  select  different  compensation  poles,  or  make 
modifications  (such  as  changes  in  the  tail-wags-dog  frequencies  as  was  previously 
discussed)  so  that  the  gyro  poles  would  become  unstable  first.  In  either  case, 
the  designer  becomes  involved  in  the  entire  regime  of  making  trades  between 
system  stability,  time  response,  sensor  location,  and  required  gains  and 
compensation. 

For  the  Honeywell  concept,  the  effects  of  the  limit  cycle  must  be  considered. 

First,  a  determination  must  be  made  of  the  tolerable  limit-cycle  amplitudes  from 
the  standpoint  of  structural  integrity.  Secondly,  the  effect  of  limit  cycle  ampli¬ 
tude  on  control  servo  power  requirements  must  be  considered.  Installed  power 
capabilities  on  ASM  configurations  tend  to  be  marginal  because  space  and  weight 
a.e  at  a  premium.  The  power  required  to  maintain  a  limit  cycle  such  as  is  re¬ 
quired  in  the  Honeywell  concept  can  become  a  predominant  factor  in  determining 
servo-power  requirements . 

The  third  adaptive  approach  is  to  program  autopilot  gain  as  a  function  of  control 
effectiveness.  A  dither  signal  is  applied  to  the  autopilot  vehicle  axis.  By  moni¬ 
toring  the  effect  of  the  signal  on  vehicle  amplitude  about  that  axis  the  control 
effectiveness  can  be  determined.  Then  if  the  required  autopilot  gains  to  maintain 
the  stability  margins  are  known  as  a  function  of  control  effectiveness  (Ng),  the  gain 


can  be  adjusted.  More  vehicle  parameter  variation  effects  can  be  removed  if  the  dither 
is  applied  to  the  axis  in  which  the  gain  adjustment  is  required.  However,  if  this  is  unde¬ 
sirable  with  some  configurations ,  it  may  La  possible  to  apply  the  dither  signal  to  one 
axis  and  adjust  the  gain  in  another  axis  (e.  g, ,  dither  roll  axis  and  adjust  pitch  axis) . 

An  acceptable  pitch  rate  gain  program  (Kg)  as  a  function  of  control  effectiveness 
(actually  shown  as  Ng/I),  for  the  baseline  Type  1  autopilot  discussed  previ¬ 
ously  is  shown  in  Figure  54.  This  program  provides  the  desired  gains  which 
were  determined  in  the  Type  1  autopilot  design  (Table  VII)  for  the  terminal  phase 
conditions  (1,  2,  and  3).  Gains  for  4,  5,  and  6  are  acceptable.  The  upper  and 
lower  gain  stability  limits  for  Flight  Conditions  1  through  6  are  also  shown.  It 
can  be  seen  that  none  of  these  points  fall  within  6  db  of  the  programmed  gain.  It 
should  be  noted  that  only  Condition  4  has  a  lower-gain  limit.  The  missile  is 
statically  unstable  at  this  condition  and  a  minimum  gain  is  therefore  required  for 
stability.  This  can  be  seen  from  the  amplitude-phase  plot  in  Figure  47.  An  ad¬ 
vantage  of  this  approach  to  an  adaptive  autopilot  is  the  possibility  of  conducting 
all  the  testing  and  monitoring  of  response  in  the  roll  axis  which  is  usually  the 
least  critical  of  the  vehicle  axes.  A  complicating  factor  is  that  relationships 
between  roll-axis  response  and  control  effectiveness  in  pitch  and  yaw  are  not  always 
straightforward. 

The  final  concept  Involves  the  calculation  of  the  autopilot  gains  from  a  single 
parameter,  M^/l.  As  shown  in  Figure  55,  there  is  an  explicit  autopilot  gain 
for  given  M^/l.  (This  brings  together  all  the  effects  of  Mach  number,  dynamic 
pressure,  angle  of  attack,  and  c.g.  shift).  This  concept  Involves  the  calculation 
of  Ma  from  measurements  of  0,  if,  z,  and  6.  References  12  and  13  establish 
the  feasibility  of  applying  the  technique  to  an  adaptive  autopilot.  It  has  been  shown 
to  be  possible  to  calculate  the  values  of  the  stability  derivatives  based  on  a  com¬ 
parison  of  the  assumed  vehicle  equations  of  motion  to  actually  measured  motions 
of  the  vehicle.  Two  equations  are  usually  needed  for  each  axis  beoause  of  the 
requirement  for  measurable  state  variables.  This  computation  can  be  carried 
out  continuously,  with  values  of  stability  derivatives  being  continuously  updated. 
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The  application  of  any  of  the  above  concepts  to  an  ASM  of  the  type  being  considered 
appears  feasible  but  would  require  considerable  analysis  and  simulation  to  con¬ 
firm  feasibility  and  to  evaluate  potential  performance  improvements.  Perform¬ 
ance  improvements  would  have  to  be  weighed  against  implementation  complexity 
to  determine  the  best  approach  for  a  given  application. 


-M«  /  I~1/SEC 

Figure  55:  TYPE  1  AUTOPILOT  BEST  GAIN  vs  M_/ 1 
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CONCLUSIONS  AND  RECOMMENDATIONS 
In  summary,  the  results  of  the  study  have  shown  that: 

(1)  Assuming  that  satisfactory  homing  sensors  are  available,  high  accuracy 
terminal  guidance  of  ASM  is  feasible  and  considerable  control  system  design 
flexibility  is  possible; 

(2)  An  optimal  guidance  system  can  be  implemented  to  meet  practical  missile 
requirements  and  characteristics,  and  can  be  implemented  with  state-of-the- 
art  onboard  digital  computers; 

(3)  Nominal  inertial  and  homing  sensor  characteristics,  which  were  selected 
for  the  study,  are  compatible  with  ASM  miss  distances  of  approximately 
5  feet; 

(4)  For  the  selected  ASM  configuration,  a  normal  acceleration  autopilot  provided 
satisfactory  inner-loop  performance  for  the  homing  phase. 

Based  upon  the  effects  of  nominal  equipment  anomalies  and  disturbances  con¬ 
sidered,  the  choice  of  homing  guidance  concept  should  be  made  from  the  stand¬ 
point  of  simplicity,  cost,  and  development  time.  In  terms  of  miss  distance,  no 
really  significant  performance  advantages  or  disadvantages  were  found  for  any 
of  the  concepts.  The  one  exception  is  that  pursuit  guidance  should  not  be  used  if 
large  bias  errors  are  present.  If  an  inertial  rate  stabilized  gimbaled  tracker  is 
used,  proportional  guidance  can  be  implemented  easily  using  the  output  of  the 
rate  gyro  for  sensor  gimbal  stabilization.  Optimal  guidance  can  be  implemented 
with  either  body-fixed  sensors  or  gimbaled  sensors. 

implementation  of  the  optimal  guidance  concept  was  shown  feasible  with  realistic 
system  characteristics.  The  preliminary  digital  computer  sizing  (Section  V) 
shows  that  several  state-of-the-art  computers  have  sufficient  capacity  for  imple¬ 
mentation  of  the  optimal  guidance  law.  Based  on  the  assumptions  made  for  sensor 
performance  and  estimates  of  the  increased  computer  requirements,  use  of  the 
optimal  guidance  laws  does  not  appear  justified  at  this  time.  The  2600-word 
memory  computer  needed  for  optimal  guidance  is  about  1500  words  larger  than 
estimates  for  implementing  the  more  conventional  proportional  or  pursuit 
guidance. 
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The  optimal  guidance  concept  showed  no  significant  advantages  over  proportional 
and  pursuit  guidance  for  the  range  of  sensor  errors  and  disturbances  that  were 
investigated.  Examination  of  optimal  guidance  error  data  shows  that  the  main 
contribution  to  the  miss  distance  is  the  error  the  optimal  (Kalman)  filter  makes 
in  estimating  the  position  deviations  normal  to  the  trajectory.  Because  the  miss 
distances  for  proportional  and  pursuit  guidance  were  almost  as  small,  it  appears 
that  the  filtering  action  of  the  vehicle  and  autopilot  dynamics  are  almost  as  good 
as  the  optimal  filtering  action  of  the  optimal  filter. 

As  more  homing  sensors  are  developed  and  performance  data  are  obtained,  it 
may  be  found  that  the  performance  of  some  of  the  sensors  may  not  be  as  good  as 
assumed  in  this  study.  However,  these  sensors  may  have  tactical  advantages 
that  make  consideration  of  them  necessary. 

V 

The  Kalman  filter  that  was  developed  for  the  study  was  based  on  white  noise; 
however,  other  filters  can  be  developed  for  nonwhite  noise  if  the  noise  can  be 
described  statistically.  Thus  the  optimal  techniques  have  the  potential  of  coping 
with  poorer  sensor  characteristics  than  the  conventional  techniques.  Trades 
could  then  be  conducted  between  sensor  performance  and  cost,  and  optimal- 
guidance  performance  and  implementation  cost.  The  results  that  show  the  op¬ 
timal  filter  to  be  the  predominate  source  of  error  suggest  the  possibility  that 
more  simple  linear  guidance  laws  could  be  used  with  the  filter  with  little  increase 
in  error. 

There  are  tactical  situations  in  which  optimal  guidance  concepts  would  have  ad¬ 
vantages.  If  there  were  long  blind  zones  such  as  might  be  eucountered  if  the 
TV  tracker  lost  the  target;  the  ability  of  the  optimal  system  to  guide  on  position 
information  from  inertial  platform  measurements  updated  while  the  target  seeker 
was  still  operating  would  allow  accurate  guidance.  If  proportional  or  pursuit 
guidance  did  not  have  enough  time  to  get  the  missile  accurately  aimed  at  the  tar¬ 
get  before  the  blind  zone,  or  if  the  blind  zone  was  long  enough  for  the  effects  of 
wind  to  be  appreciable,  proportional  or  pursuit  guidance  would  be  substantially 
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Although  Investigations  of  equipment  anomalies  indicated  no  serious  restrictions 
on  impact  accuracy,  the  results  of  the  state-of-the-art  surveys  indicate  that  the 
homing  sensor  may  be  a  major  weakness  in  the  design  of  an  ASM.  This  weakness 
is  not  necessarily  due  to  poor  performance  on  the  part  of  the  sensor  (in  fact  the 
results  of  the  TV  tracker  experiment  were  surprisingly  good)  but  is  primarily 
due  to  the  lack  of  adequate  performance  data  for  actual  trackers.  The  tracker 
error  model  (bias  and  white  noise)  used  in  this  study  was  based  on  the  best  avail¬ 
able  information  but  is  nevertheless  idealistic.  Sensor  errors  are  seldom  this 
simple —  rate  gyros  that  may  be  used  in  a  tracker  gimbal  stabilization  loop  are 
known  to  have  output  noises  that  may  contain  dominant  frequencies  as  low  as  1 
to  2  cps.  This  type  of  noise,  because  it  approaches  the  missile  control  frequency 
regime,  makes  a  high  accuracy  control  problem  considerably  more  difficult 
than  the  white  noise.  Even  predominantly  higher  frequency  noise  can  cause  dif- 
culties  with  a  digital  flight  control  system. 

The  number  of  general  conclusions  that  can  be  made  about  the  normal  accel¬ 
eration  autopilot  design  are  limited  because  they  are  very  much  configuration 
dependent.  However,  some  precautions  should  be  mentioned.  Basic  designs 
should  not  be  predicated  on  simplified  autopilot  representations  in  determining 
autopilot  response  criteria  to  meet  miss  distance  requirements.  In  the  con¬ 
figuration  investigated,  the  actual  autopilot  response  was  well  above  that 
required  to  make  response  effects  negligible.  Details  such  as  body  bending 
and  guidance  gains  must  also  be  considered  in  the  autopilot  design  because  of 
their  interrelationship  with  overall  system  stability  and  time  response.  Some 
type  of  programmed  gain  or  adaptive  configurations  will  be  required  for  the 
normal  acceleration  autopilot  because  of  the  wide  range  of  flight  conditions  to 
which  the  homing  ASM  is  subjected  and  the  rather  stringent  response  require¬ 
ments  during  the  terminal  phase.  As  discussed  in  Section  VI,  a  programmed 
gain  autopilot  would  not  be  easy  to  Implement.  The  bistable  controller  concept 
that  was  investigated  appears  to  have  considerable  promise.  It  can  be  imple¬ 
mented  with  available  components  and  it  provides  good  response  and  steady-state 
error  characteristics. 
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To  provide  additional  information  necessary  to  the  design  of  satisfactory  homing 

ASM's,  the  following  areas  of  investigation  are  recommended: 

(1)  Better  definition,  under  realistic  operating  conditions,  of  the  performance 
characteristics  of  homing  sensors  that  are  potentially  useful  from  a  tactical 
standpoint  is  necessary,  information  is  required  on  acquisition  capability, 
response  characteristics,  and  the  statistical  and  real-time  characteristics 
of  the  output  errors.  The  output  error  data  should  provide  bias  error  dis¬ 
tribution  and  noise  power  spectral  density  characteristics.  Power  spectral 
density  information  should  cover  the  very  low  frequency  range  (less  than 

1  cps)  to  include  bias  shift  effects.  Tape  recordings  of  tracker  noise  would 
be  useful  in  simulation  studies.  Acquisition  range  capability  is  needed  to 
define  midcourse  guidance  accuracy  requirements;  more  detailed  blind 
range  information  is  necessary  to  select  blind  range  guidance  concepts. 

(2)  More  complete  hybrid  computer  simulations  of  optimal  control  concepts 
should  be  conducted  in  which  its  impact  on  the  autopilot  can  be  determined 
(and  vice  versa)  and  effects  of  items  such  as  autopilot  response,  command 
limiting,  and  servo  limiting  can  be  evaluated.  Control  power  requirements 
should  be  compared  with  those  for  proportional  and  pursuit  guidance. 

(3)  Additional  autopilot  studies  in  which  cross-coupling  effects  are  included 
should  be  conducted  to  compare  relative  merits  of  the  bistable  controller 
and  other  adaptive  concepts.  This  effort  should  include  an  investigation  of 
the  use  of  bistable  controllers  with  more  than  one  bistable  signal.  Com¬ 
patibility  of  the  homing  phase  and  midcourse  phase  autopilots  with  the  boost 
phase  control  requirements  should  be  investigated. 
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APPENDIX  I 


STATE-OF-THE-ART  SURVEY 

1.  Introduction.  State-of-the-art  surveys  were  conducted  in  the  areas  of  in¬ 
ertial  guidance  sensors,  optical  homing  guidance  sensors,  and  radar-type  active 
and  passive  homing  sensors.  Original  intentions  to  include  current  ASM  control 
concepts,  terminal  accuracy,  and  effects  of  disturbances  on  terminal  accuracy 
in  the  survey  were  frustrated  by  the  inaccessibility  of  classified  material. 

Three  basic  inertial  navigation  sensor  schemes  were  studied:  strapdown,  floated 
ball  (MIT  "flimbal"  concept),  and  the  gimbaled  platform.  The  survey  resulted 
In  the  selection  of  the  conventional  gimbaled  platform  type  inertial  measurement 
unit  as  a  baseline  sensor  for  the  study.  An  analysis  of  the  errors  generated 
during  midcourse  prior  to  target  acquisition,  intrinsic  to  this  system,  for  three 
representative  air-to-ground  trajectories  and  three  cost  level  instrumentations 
was  accomplished.  The  error  data  generated  from  this  analysis  were  used  as 
initial  position  and  velocity  error  data  for  the  miss  distance  analyses. 

The  results  of  the  survey  of  optical  sensors  indicated  that  all  such  systems  (TV 
or  infrared),  both  current  and  projected  lo'r  the  1970  period,  employ  manual 
target  acquisition.  In  some  systems  the  target  is  acquired  before  launch  from 
the  carrier  aircraft;  in  others  the  sensor  output  is  telemetered  to  the  carrier 
aircraft  and  the  necessary  signals  for  manual  target  acquisition  and  lockon  are 
transmitted  back  to  the  missile.  Once  acquisition  is  achieved  the  system  is  no 
longer  dependent  on  a  man  in  the  loop  and  the  terminal  phase  of  the  flight  is 
completed  automatically.  Because  this  study  was  concerned  primarily  with  the 
terminal  homing  phase  of  flight,  studies  performed  for  optical  homing  sensors 
did  not  include  investigation  of  the  dynamics  and  pecularities  of  manual  target 
acquisition. 

A  gimbaled  TV  tracker  using  centroid  tracking  logic  was  chosen  as  the  baseline 
sensor  because  insufficient  data  was  available  to  evaluate  the  correlation  type 
tracker.  To  get  quantitative  effects  of  target  and  background  characteristics 
and  range  to  target  on  tracker  errors,  an  experiment  was  performed  using  the 
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centroid  tracking  concept  and  photographs  of  various  targets  taken  from  different 
altitudes.  Simulated  homing  flights  yielded  the  effect  of  range-to- target  on 
tracker  error.  The  effect  of  target  contrast  was  not  determined. 

The  study  results  showed  that  radar  cannot  provide  an  operationally  useful  sys¬ 
tem  against  any  wide  variety  of  targets,  due  to  target  scintillation  characteristics 
and  background  reflection  characteristics  of  nearby  objects.  Because  active 
radar  sensors  appear  applicable  only  to  very  restricted  target  situations,  radar 
sensors  were  not  considered  in  the  system  studies. 

2.  Inertial  Measurement  Unit.  The  baseline  inertial  guidance  sensor  chosen  for 
further  study  was  the  conventional  gimbaled  type  of  inertial  platform,  for  which 
information  was  available  on  cost  and  accuracy.  Equipment  variations  for  this 
type  of  inertial  measurement  unit,  relating  cost  to  sensor  accuracy,  were  exer¬ 
cised  with  an  upper  limit  on  cost  of  the  platform  and  associated  electronics  of 
$40,000  (excluding  airborne  computer). 

The  primary  function  of  the  inertial  guidance  system  is  to  furnish  boost  and  mid¬ 
course  navigation  information  for  guidance  purposes,  and  also  in  some  instances 
to  provide  information  on  missile  motion  during  the  homing  phase. 

To  perform  these  functions  the  inertial  navigation  system  and  its  associated 
navigation  and  guidance  computer  muat  have  the  following  characteristics:  (1) 
unlimited  azimuth  freedom  to  handle  missile  launch  In  any  direction;  (2)  Infor¬ 
mation  outputs  of  missile  acceleration,  velocity,  and  attitude;  (3)  guidance  and 
navigation  accuracy  sufficient  to  ensure  target  acquisition  at  the  end  of  the  mid¬ 
course  phase  of  flight;  (4)  compatibility  with  the  weapon  carrier's  master  navi¬ 
gation  system;  and  (5)  be  simple,  reliable,  inexpensive,  light,  small,  and  have 
a  short  starting  time. 

The  three  basic  inertial  navigation  schemes  studied  were  strapdown,  floated 
ball,  and  gimbaled  platform. 

In  the  strapdown  or  analytic  system,  the  inertial  sensors,  consisting  of  rate 
gyros  and  accelerometers,  are  mounted  directly  to  the  missile  structure.  The 
inertial  attitude  of  the  missile  is  determined  by  Integrating  the  rate  gyfro  output; 
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missile  position  and  velocity  in  inertial  coordinates  are  computed  from  the 
accelerometer  outputs,  taking  into  account  the  time  varying  orientation  of  the 
accelerometer  input  axis  with  respect  to  inertial  space.  The  advantages  of  the 
strapdown  system  include:  small  size,  weight,  and  power  requirement,  all  of 
which  result  from  the  absence  of  the  gimbal  structure,  angular  resolvers,  tor¬ 
que  rs,  and  slip  ring  of  gimbaled  systems.  The  system  is  all- attitude  and  certain 
instrument  error  coefficients  can  be  minimized  by  properly  orienting  the  instru¬ 
ments  with  respect  to  the  missile  nominal  acceleration  vector. 

The  disadvantages  of  the  strapdown  system  reside  in  the  large,  high-speed 
computer  to  conduct  the  coordinate  transformation  necessitated  by  the  non- 
commutivity  of  angular  displacement;  the  stringent  demands  on  range  and  reso¬ 
lution  of  the  inertial  instruments,  particularly  the  gyros;  the  susceptibility  of 
the  platform  to  base  motion  coupling,  coning  induced  errors,  misalignment 
introduced  by  thermal  base  motion  coupling,  and  complexity  of  the  airplane 
navigation  interface. 

The  floated  ball  platform  consists  of  displacement  gyroB  and  accelerometers 
mounted  in  a  spherical  ball,  neutrally  buoyant  in  a  fluid  filled  cavity.  External 
communication  with  the  ball  is  provided  by  a  set  of  brush  contacts  to  the  ball, 
multiplexed  communication  signals  being  modulated  on  the  d.c.  power  supply. 

The  platform  is  all-attitude,  compact,  light,  and  with  moderate  power  consump¬ 
tion.  The  stabilized  ball  isolates  the  inertial  instrument  from  the  base  vibration 
and  coning  motion. 

The  disadvantages  of  the  floated  ball  include  the  complex  signal  connection  to  the 
ball,  non-Euler  angle  attitude  sensing,  and  inability  to  achieve  preferential  ori¬ 
entation  of  inertial  sensor  with  respect  to  the  missile  acceleration  vector.  The 
floated  ball  platform  was  regarded  as  high  risk,  because  it  had  not  yet  been  flown. 

The  advantages  of  the  conventional  gimbaled  platform  are  low  risk,  absence  of 
extreme  requirements  on  the  inertial  sensory  compatibility  with  airplane  navi¬ 
gator  Euler  angle -attitude  readout,  minimum  computational  requirements 
(spatial  Integrations  performed  mechanically  by  gimbals),  and  high  degree  of  base 


motion  isolation.  By  use  of  a  four-gimbal  arrangement,  the  platform  can  be 
made  all-attitude. 

Its  disadvantages  included:  size,  weight,  and  power  consumption;  cooling 
requirement;  and  inability  to  achieve  optimal  orientation  of  sensors. 

Three  systems  were  compared  in  a  somewhat  arbitrary  quantitative  manner  in 
the  "Opinion  Table,"  Table  IX.  Each  performance  parameter —  risk,  weight, 
volume,  and  power  consumption —  is  assigned  a  weighting  factor  of  up  to  10  (the 
"maximum  weight"  column)  based  on  its  considered  importance.  Each  of  the 
candidate  systems  is  assigned  a  "relative  weight"  of  up  to  10  that  represents  its 
merit  for  each  of  the  performance  parameters.  The  product  of  the  "maximum 
weight"  and  the  "relative  weight"  represents  the  score.  Total  score  is  the  sum 
of  the  scores  for  the  performance  parameters.  As  shown  in  Table  DC,  the  con¬ 
ventional  gimbaled  platform  received  the  highest  score.  It  was  selected  as  the 
baseline  inertial  navigation  system. 

The  system  selected  for  this  study  consists  of  an  external,  three-gimbal,  roll  - 
limited  stabilized  platform  on  which  are  mounted  an  orthogonal  triad  of  accel¬ 
erometers.  (A  roll-stabilized  missile  is  assumed. )  The  outputs  of  these  accel¬ 
erometers  are  singly  and  doubly  Integrated  to  determine  the  components  of 
missile  position  and  velocity.  There  are  a  number  of  methods  of  mechanizing 
this  system.  Because  the  system  will  probably  be  tracking  the  local  vertical 
before  launch.  It  may  be  well  to  continue  this  type  of  operation  during  missile 
guidance  to  eliminate  starting  transients,  due  to  the  step  change  at  launch,  from 
one  type  of  mechanization  to  another.  However,  oilier  mechanizations  such  as  a 
"space  fixed"  platform  coordinate  Bystem,  with  a  coordinate  transformation  in 
the  navigation  computer,  may  prove  to  be  worth  its  cost,  because  the  use  of 
inaccurate  torquing  components  (intrinsic  to  the  low-cost  grade  platform  con¬ 
sidered  here)  can  be  eliminated.  Final  selection  of  mechanization  should  include 
the  trade  between  errors  and  cost  to  mechanize  and  compute  the  data  for  sensor 
pointing  from  the  platform  coordinates. 


142 


Table  IX:  OPINION  TABLE 


Strapdown  Floated  Ball  Gimbaled 


Maximum 

Relative 

Relative 

Relative 

Parameter 

Weight 

Weight 

Score 

Weight 

Score 

Weight 

Score 

Risk 

10 

5 

50 

3 

30 

10 

100 

Weight,  volume. 

power  consumption 

Extreme  sensor 

3 

10 

30 

8 

24 

4 

12 

requirements 

5 

3 

15 

TO 

50 

10 

50 

Computer  require- 

ments 

Cooling  require¬ 

10 

3 

30 

7 

70 

10 

100 

ments 

2 

10 

20 

7 

14 

4 

8 

Airplane  system 

interface 

5 

4 

20 

8 

40 

10 

50 

Signal  path  equipment 

8 

10 

80 

4 

32 

8 

64 

Base  motion  isolation 

8 

5 

40 

10 

80 

10 

80 

Preferential  instru- 

ment  alignment 

5 

10 

50 

5 

25 

5 

25 

Cost 

10 

10 

100 

_5 

_50 

_7 

70 

TOTAL 

435 

415 

559 

NOTE:  •  The  higher  the  score  the  more  favorable  the  characteristic. 

•  Greatest  possible  weight  =  10. 

•  Product  of  maximum  weight  for  each  parameter  and  relative 
weight  equals  the  score. 


Three  different  trajectories  were  considered  in  the  analysis:  a  long-range 
semiballistic  trajectory,  a  medium-range  low-altitude  skip  trajectory,  and  a 
short- range  semiballistic  trajectory. 

Three  cost  levels  of  IMU  were  considered  in  the  error  analysis  for  each  of  the 
trajectories.  Table  X  presents  the  cost  breakdown  for  each  of  the  IMU's. 


Table  X:  IMU  COST  BREAKDOWN 


Three  Gyros 

$  600 

$  3,000 

$  3,000 

Three  Accelerometers 

1,200 

3,000 

9,000 

Platform  Navigation  Computer 

9,200 

16,000 

18,000 

Electronics 

2,000 

3,000 

4,000 

TOTAL 

$13,000 

$25,000 

$40,000 

The  two  main  sources  of  error  in  this  system  are  initial  condition  errors  and 
instrumentation  errors.  One  representative  set  of  data  for  initial  condition 
errors  at  launch  was  used.  There  are  two  independent  sources  of  initial  head¬ 
ing  error;  the  heading  error  introduced  by  the  weapon  carrier,  and  the  error 
from  the  slaving  misalignment  of  the  inertial  platform.  These  errors  were 
assumed  uncorrelated.  Their  effect  was  to  produce  a  cross-range  error.  Initial 
misalignment  of  the  platform  vertical,  though  considered  in  the  detailed  analysis 
of  errors,  produced  negligible  errors.  The  mapping  error  was  considered  to  be 
1  foot  per  nautical  mile.  To  study  instrumentation  error  propagation  in  the 
inertial  navigation  system  and  assess  the  overall  acouracy  of  the  weapon  system, 
a  previously  developed  digital  computer  inertial  navigation  error  analysis  pro¬ 
gram  was  used.  The  inertial  components  and  their  characteristics  are  tabulated 
in  Table  XI.  Briefly,  the  system  gyro  and  acceleration  analysis  error  models 
included  were:  la  uncertainty  in  hccelerometer  bias,  (b);  the  la  uncertainty  i 

i 

in  accelerometer  scale  factor,  (SF);  the  ia  uncertainty  in  accelerometer  cross¬ 
axis  sensitivity,  (NL);  the  la  uncertainty  in  gyro  random  drift,  (RD);  the  la  ( 

uncertainty  in  mass  unbalance,  (MU);  and  the  la  uncertainty  in  gyro  aniso  elas-  ! 

\ 

tlcity,  (Aniso).  Each  term  of  the  gyro  error  model  was  integrated  separately  \ 
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Table  XI:  VARIOUS  COST  LEVEL  IMU  COMPONENTS 


Aniso~  lo  uncertainty  in  gyro  aniso  elasticity 


to  give  the  components  of  la  uncertainty  in  the  platform  tilt  angle,  and  these 
were  then  introduced  into  the  net  acceleration  model.  Each  term  was  singly 
and  doubly  integrated  to  yield  the  influence  of  each  error  source  on  the  overall 
system  velocity  and  position  accuracy.  Further,  the  program  determined  the 
square  root  of  the  expected  value  of  the  sum  of  squares  of  these  errors.  It  was 
assumed  that  there  were  no  correlated  uncertainties  in  the  error  sources;  hence, 
the  result  was  the  root  square  sum  (RSS)  of  the  errors.  The  hand-calculated 
initial  condition  errors  were  combined  with  the  computer-determined  1MU  errors 
and  the  total  error-at-acquisition  budgets  that  resulted  are  tabulated  in  Tables 
XII,  xni,  and  XIV.  Note  that  in  the  cosi  range  of  $13,000  to  $40,000,  the  accu¬ 
racy  of  the  inertial  navigation  system  is  not  substantially  increased  for  cost 
levels  beyond  $25,  000. 

3.  TV  Honflng  Sensors.  Vendor  data  on  optical  homing  sensors  showing  the 
effect  of  range  to  target,  and  target  versus  background  characteristics  were 
required  as  a  basic  input  to  the  formulation  of  the  optimal  guidance  filter,  but 
the  state-of-the-art  survey  did  not  yield  these  data.  However,  sufficient  de¬ 
scriptive  information  on  TV- type  centroid  trackers  was  obtained  to  enable  such 
accuracy  characteristics  to  be  evaluated  on  a  semi-idealized  basis.  No  infor¬ 
mation  was  found  that  would  allow  a  credible  prediction  of  accuracy  for  corre¬ 
lation-type  trackers;  therefore,  the  TV-type  tracker  using  centroid-tracking 
logic  was  selected  as  the  baseline  sensor  to  be  considered.  Tracker  output 
errors  were  assumed  to  consist  of  bias  error  and  white  noise.  From  all  avail¬ 
able  data,  the  standard  deviation  of  white  noise  was  from  1  to  3  milliradians. 

The  nominal  value  was  chosen  to  be  1  milliradian.  The  nominal  value  of  bias 
error  was  assumed  to  be  1.  74  milliradians. 

The  literature  search  on  TV  sensors  revealed  data  for  only  narrow- field-of- view 
devices  that  were  intended  for  use  with  an  inertial  attitude  or  attitude  rate 
stabilized  gimbaled  mount.  The  wide  field  of  view  devices,  for  which  data  was 
obtained,  provided  the  wide  field  of  view  capability  during  the  acquisition  phase 
only,  and  had  only  a  narrow  field  of  view  capability  in  the  tracking  mode.  Due 
to  lack  of  information  on  the  characteristics  of  wide  field  of  view  trackers,  and 
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Table  XII:  la  ERRORS  FOR  SHORT-RANGE  SE MIB ALLIST1C  TRAJECTORY 


Cost  of  IMU  $13,000  $25,000  $40,000 

Slant  Range 

Acquisition  Total  Error  at  Acquisition 


i^ftl 


15,000 

1065 

1010 

1010 

X 

30,000 

1060 

1010 

1009 

60,000 

1050 

1009 

1005 

15,000 

1091 

1015 

1015 

Position 

y 

30,000 

1091 

1015 

1013 

(~-tt) 

60,000 

1091 

1010 

1009 

15,000 

670 

285 

252 

Z 

30,000 

630 

280 

248 

60,000 

560 

265 

240 

15,000 

5.90 

4.02 

1.27 

X 

30,000 

5.74 

3.98 

1.27 

60,000 

5.42 

3.91 

1.25 

la 

Velocity 

, 

15,000 

19.9 

15.0 

16.8 

Errors 

(ft/sec) 

y 

30,000 

18.  5 

15.0 

16.3 

60,000 

15.  8 

15.0 

15.2 

15,000 

12.  70 

1.69 

2.75 

Z 

30,000 

10.90 

1.  71 

2.21 

60,000 

7.40 

1.74 

1.14 

15,000 

3.030 

0.073 

0.060 

®x 

30,000 

2.  820 

0.070 

0.056 

60,000 

2.400 

0.064 

0.048 

Table 

Tilt 

15,000 

3.  030 

0.074 

0.  063 

30, 000 

2.  820 

0.084 

0.058 

Ry 

60,000 

2.  400 

0.104 

0.049 

(—degree) 

15,000 

3.030 

0.068 

0.047 

ez 

30,000 

2.820  . 

0.080 

0.047 

60,000 

2.400 

0.102 

0.049 

.4“ 


Table  XIII:  la  ERRORS  FOR  LOW- ALTITUDE  SKIP  TRAJECTORY 


Coat  of  IMU  $13,000  $25,000  $40,000 


Slant  Range 

Acquisition  Total  Error  at  Acquisition 

-(-ft) _ 


15,000 

6100 

1390 

1300 

X 

30,000 

4900 

1250 

1210 

60,000 

2500 

1150 

1100 

Position 

15,000 

5000 

1940 

1900 

(~  ftl 

y 

30,000 

4500 

1820 

1720 

60,000 

4000 

1600 

1350 

15,000 

560 

450 

348 

Z 

30,000 

525 

420 

347 

60,000 

430 

380 

347 

15,000 

72.0 

11.3 

10.8 

X 

30,000 

64.0 

10.0 

9.5 

60,000 

55.0 

8.5 

8.0 

lo 

Velocity 

15,000 

72.0 

13.6 

13.2 

Errors 

(ft/sec) 

y 

30,000 

64.0 

12.5 

12.2 

60,000 

57.0 

10.7 

10.5 

15,000 

10.0 

2.6 

1.1 

z 

30,000 

3.6 

1.5 

1.0 

1 

60, 000 

0.3 

0.2 

0.1 

15,000 

7.  460 

0.430 

0.144 

Ox 

30,000 

7.110 

0.410 

0.139 

60,000 

6.620 

0.370 

0. 13% 

Table 

Tilt 

15,000 

7.460 

0.430 

0.144 

8v 

30,000 

7. 110 

0.410 

0.136 

degree) 

60,000 

6.620 

0.370 

0. 132 

15,000 

7.150 

0.430 

0.144 

0z 

30,000 

6.840 

0.410 

0.139 

60,000 

6.450 

0.370 

0.132 
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Table  XIV:  J  a  ERRORS  FOR  LONG-RANGE  SEMIBALLISTIC  TRAJECTORY 


Cost  of  IMU 


$13,000  $25,000  $40,000 


Slant  Range 

Acquisition  Total  Error  at  Acquisition 

_ LzM _ 


15,000 

3790 

1200 

1082 

K 

30,000 

3640 

1190 

1075 

60,000 

3330 

1160 

1060 

Position 

15,000 

2990 

2520 

2480 

(~ft; 

y 

30,000 

2890 

2460 

2410 

60,000 

2700 

2320 

2300 

15,000 

1950 

530 

339 

Z 

30,000 

1930 

525 

334 

60,000 

1890 

500 

240 

15,000 

32.0 

5.3 

4.4 

X 

30,000 

29.0 

4.7 

3.9 

60,000 

23.  3 

3.6 

2.8 

15,000 

32.4 

12.  7 

11.8 

la 

Velocity 

y 

30,000 

29.6 

12.4 

11.7 

Errors 

(ft/sec) 

60,000 

24.0 

12.0 

11.6 

15,000 

15.7 

2.7 

3.4 

z 

30,000 

13.3 

2.6 

3.1 

60,000 

8.  7 

2.5 

2.6 

15,000 

6.39 

0.30 

0.13 

8v 

30,000 

6.28 

0.28 

0.13 

60,000 

6.06 

0.26 

0.12 

Table 
Tilt 
Angle 
(~  degree) 

15,000 

6.39 

0.30 

0. 13 

% 

30,000 

6.28 

0.29 

0.13 

oC,000 

6. 06 

0.27 

0.12 

15,000 

6.24 

0.30 

0.13 

a. 

30,000 

6.14 

0.29 

0.13 

60,000 

5.  93 

0.27 

0.12 

i 


also  the  apparent  equivalence  of  body- fixed  and  gimbaled  sensors  as  conceptually 
applied  to  the  optimal  guidance  task,  the  studies  were  based  on  the  characteristics 
of  narrow-field-of-view  gimbaled  TV  trackers  only.  A  simplified  block  diagram 
of  a  single  axis  implementation  of  this  type  of  tracker  iB  shown  in  Figure  56. 


Figure  56:  TV  TRACKER  BLOCK  DIAGRAM 


An  experimental  program  was  undertaken  to  generate  accuracy  characteristics 
for  a  centroid  tracking  sensor  with  scan  resolution  typical  of  the  Walleye/Condor 
trackers.  To  get  quantitative  effects  of  target  and  background  characteristics 
and  range  to  target  on  tracker  errors,  simulated  homing  flights  were  made  using 
photographs  of  various  representative  targets  taken  from  five  different  altitudes. 
Results  of  tho  range-to-target  effect  on  tracker  error  were  obtained,  but  the 
effect  of  target  contrast  was  not  determined. 


a.  TV  Tracker  Survey  Summary.  The  optical  sensors  considered  included 
imaging  and  nonimaging  systems.  The  general  shortcomings  of  nonimaging 
systems  tend  to  eliminate  them  quickly  from  further  consideration.  Mechanical 
complexity  is  one  major  drawback.  In  order  to  achieve  substantial  background 
rejection  capability  with  a  nonimaging  system,  complex  modulation  and  demodu¬ 
lation  techniques  are  necessary  in  the  electronic  signal  processing,  including 
the  use  of  reticles  and/or  nutating  devices.  Reticle  systems  are  fixed  to  one 
set  of  operating  parameters  —  such  as  field  of  view,  target  size,  modulation 
frequencies,  and  scan  pattern.  Reticles  and  reticle  drive  systems  with  their 
attendant  bulkiness  and  moving  parts  can  be  eliminated  by  use  of  the  vidicon. 

Imaging  tubes  with  electronic  scan  allow  the  use  of  some  rather  straightforward 
signal  processing  to  determine  target  size,  number  of  targets,  and  spatial  posi¬ 
tion  in  a  given  raster  configuration.  Independent  horizontal  and  vertical  deflec¬ 
tion  systems  permit  wide  variation  in  programming  scan  rates  and  raster  patterns 
to  optimize  tracking  capabilities.  A  literature  search  produced  a  great  number 
of  references  to  TV  trackers  and  related  systems.  Examination  of  the  available  j 

material  was  begun  with  a  study  of  the  Navy's  Walleye  missile  guidance  scheme  j 

(an  edge-type  tracker),  the  Norden  proposal  on  a  similar  device  for  the  Condor  i 

missile,  centroid-type  trackers  by  Southern  Research  Institute,  Barnes  Engi-  i 

neering  Company,  and  a  few  others. 

A  typical  TV  tracker  includes  a  camera,  tracker  circuits,  monitor,  and  a  servo 
gimbal  drive  system  to  mechanically  close  the  loop.  For  special  applications, 
accessories  such  as  image  intensifiers,  telescopic  lenses,  optical  filters,  auto¬ 
matic  iris,  and  electronic  exposure  controls  can  be  added  to  the  system.  Target 
position  information  is  derived  by  locating  video  information  relative  to  the 
raster  scanned.  Because  the  beginning  of  each  horizontal  and  vertical  scan  is  • 

known,  the  position  of  the  target  can  be  determined  by  measuring  the  period  of 
time  from  the  start  of  a  frame  and  line  until  the  appearance  of  the  video  pulse 
indicating  target  presence.  Video  pulses  are  located  relative  to  reference  gen-  | 

erator  outputs  that  are  synchronized  with  the  horizontal  and  vertical  sweep  wave-  j 

forms  of  the  camera.  The  reference  generator  outputs  are  used  to  create  x  and  y 
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tracking  gate  crosshairs  that  can  be  positioned  anywhere  within  the  field  of  view 
by  adjusting  a  delay  in  the  sync  information. 

Background  discrimination  is  accomplished  through  the  use  of  gating,  to  blank 
out  a  large  portion  of  the  raster.  This  allows  a  primary  target  to  be  tracked 
while  other  targets,  decoys,  and  background  details  in  the  field  of  view  are  ig¬ 
nored.  Blanking  is  synchronized  with  target  video  so  that  the  video  gate  auto¬ 
matically  tracks  the  target  as  it  moves  within  the  field  of  view.  Video  process¬ 
ing  consists  of  amplification,  clipping,  smoothing,  and  removal  of  undesired 
background  signal.  In  contrast  systems,  the  target  video  is  selected  on  the  basis 
of  amplitude  relative  to  the  background  by  clipping  away  the  lower  amplitude 
information.  Video  peaks  are  then  amplified  and  used  to  trigger  a  position  detector 
flip-flop  circuit.  Target  position  is  determined  from  the  output  pulse  width  of 
the  position  detector  that  has  been  triggered  by  the  target  video  pulse  and  the 
delayed  sync  pulse.  This  information  is  then  integrated  and  amplified  to  provide 
a  d.c.  voltage  output  proportional  to  the  position  of  the  target. 

The  Walleye  tracker,  being  well  along  in  development,  is  probably  most  repre¬ 
sentative  of  the  edge  tracker  state  of  the  art.  One  large  drawback  of  tracker 
systems  that  operate  on  the  edge  or  point  features  of  a  target  is  the  tendency  to 
track  the  corners  of  targets  or  to  shift  to  a  false  target,  such  as  a  bright  point. 
This  may  make  them  highly  susceptible  to  decoying  and  ineffective  against  a 
target- background  complex  containing  many  pointB  of  equal  reflectance.  The 
edge  tracking  philosophy  is  such  that  a  target  will  be  tracked  by  the  edge  that  is 
scanned  earliest,  rather  than  as  a  whole.  Therefore,  at  close  range,  the  edge 
or  corner  of  the  target  becomes  the  aim  point,  and  central  hits  are  of  low 
probability. 

There  are  two  concepts  of  centroid-type  tracking.  One  uses  the  spatial  distri¬ 
bution  of  the  target  center.  When  used  as  the  aimpoint,  tHs  pseudocentroid 
should  produce  central  hits.  This  concept  iB  known  as  "contrast  tracking"  as 
it  is  based  on  an  area-gated  scheme  that  computes  the  location  of  the  centroid 
oi  the  weighted  video  information  lying  inside  the  gate.  The  area-gated  devices 


152 


do  not  track  the  center  of  the  target  proper  aB  the  contrast- contour  trackers 
are  able  to  do. 

The  Ford-Aeronutronies  Multiple  Edge  Adaptive  Tracking  (MEAT)  tracker  and 
the  Southern  Research  Institute  (SRI)  tracker  are  most  representative  of  the 
state  of  the  art  in  contrast- contour  trackers.  The  MEAT  tracker  uses  an  aster¬ 
isk-shaped  scan  pattern  that  can  be  generated  from  the  center  to  outside  edges 
of  the  target,  or,  as  a  result  of  recent  development,  a  reverse  outside  to  edge 
asterisk  pattern  can  be  used.  Thus,  the  tracker  operator  has  the  option  of 
choosing  scan  patterns  to  accommodate  different  target-background  situations. 

The  MEAT  tracker  logic  can  locate  an  apparent  centroid  with  as  little  as  120 
degrees  of  target  periphery  defined  by  a  workable  contrast  edge.  The  SRI 
tracker  has  evolved  through  a  number  of  gating  schemes,  including  a  manually 
set  gate,  a  self- fitting  rectangular  gate,  and  the  most  recently  developed  con¬ 
tour-following  gate  with  shape  memory  that  offers  the  best  immunity  from  mo¬ 
mentary  loss  of  contrast  on  any  edge  of  the  target. 

Correlation- type  trackers  are  the  least  known.  This  method  has  the  potential 
of  being  best  for  use  against  poorly  defined  targets,  because  information  from 
the  entire  scene  is  used.  On  the  other  hand,  a  correlation  tracker  with  a  large 
field  of  view  would  be  ineffective  against  a  target  in  a  uniform  background.  Two 
correlation  tracking  schemes  were  encountered  during  the  survey.  One,  by 
General  Electric,  had  a  field  of  view  smaller  than  the  target  and  a  memory  logic 
capable  of  updating  itself  once  per  second.  The  other  was  a  system  being  devel¬ 
oped  by  General  Precision  Labs  and  Sperry.  This  system  had  a  25-degree  field 
of  view  and  used  a  preinstalled  photo  transparency  of  the  target  scene  as  the 
reference.  The  reference  photo  transparency  is  rotated  back  and  forth  and,  as 
the  image  of  the  scene  on  the  ground  aligns  with  the  photo,  a  bright  correlation 
spot  is  observed.  The  correlation  spot  represents  the  condition  in  which  maxi¬ 
mum  correlation  appears  between  the  bright  points  on  the  ground  and  transparent 
portions  of  the  phototransparency.  At  present,  correlation  devices  have  been  de¬ 
veloped  mainly  as  a  tar  get  acquisition  aid  for  the  pilot.  A  great  deal  more  complexity 
mustbe  added  to  convert  it  into  a  tracker.  At  the  time  of  this  survey,  no  company  had 
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implemented  a  tracker  from  this  concept  as  yet,  although  Emerson  Electric  is 
reported  to  haye  sol .  i  d  the  problem  with  a  more  complex  unit  using  magnetic 
drum  storage  and  a  mechanical  scanning  detector. 

A  summary  of  some  of  the  tracking  techniques  known  is  shown  in  Table  XV. 


Table  XV:  SURVEY  OF  OPTICAL  TRACKING  SYSTEMS 


System 

Military 

Agency 

Tracking 

Discipline 

Tracking 

Characteristics 

Sensor 

Gimbaled 

Sensor 

Field  of 
View 

Remark* 

Walleye 

Condor 

MEAT 

(Phllco) 

Goes  I 

— 

(Phllco) 

Southern 

Research 

Institute 

Ooos  11 
(Q.E.) 

"Hogs" 

(LTV) 

— 

Classified  Data — See  Classified  Supplement,  Figure  8 

This  represents  information  gathered  through  reports  and  visits  to  vendors.  In 
general,  the  majority  of  equipment  being  developed  is  not  beyond  the  breadboard 
stage.  Little  effort  has  been  placed  on  the  mechanics  and  servo  system  involved 
for  gimbaled  use  in  a  missile. 

b.  TV  Tracker  Experiment.  Almost  no  testing  has  been  done  that  would  pro¬ 
duce  usable  performance  or  accuracy  data.  The  material  surveyed  for  the  sub¬ 
ject  contract,  and  visits  with  vendors,  yielded  no  quantitative  information  on  the 
effects  of  target  and  background  characteristics  (i.e. ,  target  detail,  contrast, 
range,  etc.)  on  tracker  errors.  Because  this  information  was  required  as  a  basic 
input  to  the  optimal  guidance  filter  synthesis,  a  program  was  initiated  to  obtain  a 
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reasonable  estimate  of  tracker  errors  due  to  target  and  background  character¬ 
istics  through  the  simulation  of  the  information  processing  characteristics  of  a 
target  centroid  tracker. 

Aerial  transpositive  photographs  were  obtained  for  12  selected  ground  targets 
at  five  different  altitudes.  A  reference  mask  was  photographically  incorporated 
around  each  target,  or  group  of  targets,  to  be  scanned  by  the  simulator.  Targets 
selected  varied  from  sharply  defined  high-contrast  targets  to  poorly  defined  low- 
contrast  targets  approaching  the  marginal  limits  of  manual  target  acquisition 
capability.  These  are  shown  in  Figures  57  and  58. 

The  experimental  equipment  consisted  of  a  flying  spot  scanner  arrangement, 
shown  in  Figures  59  and  60.  Each  aerial  photograph  was  scanned  by  a  250-line, 
1-inch  by  1-inch  raster,  representing  a  25-milliradian  square  field  of  view  that 
gives  a  scan  resolution  typical  of  the  Walleye/Condor  trackers.  A  wide-band 
FM  tape  recorder  was  used  to  record  the  scanner  video  outputs,  the  sweep  gen¬ 
erator  output,  and  the  mask  reference  pulses.  The  analog  data  was  then  digi¬ 
tally  encoded  for  processing  by  a  digital  analysis  program  that  analyzed  the  data 
using  centroid- tracker  logic. 

N. 

The  contour  centroid  of  each  target  was  precisely  measured  with  a  Vanguard 
motion  analyzer  having  an  accuracy  of  0.026%  of  the  25-milliradian  field-of-view 
dimension.  The  deviation  of  the  computed  centroid  from  the  actual  centroid  was 
determined  for  each  target,  based  on  the  difference  between  the  centroid  as  com¬ 
puted  using  the  digital  analysis  of  video  scan  data,  and  the  centroid  as  measured 
by  the  Vanguard  motion  analyzer. 

Analysis  of  the  resulting  data  determined  the  effect  of  altitude  on  the  rms  tracker 
error  as  shown  in  Figure  61.  Each  data  point  represents  tracking  errors  for  a 
minimum  of  10  different  targets.  Two  results  should  be  noted.  First,  the  range 
of  rms  error  is  below  the  1.  74  milliradians  selected  as  a  1CT  error  for  the  study. 
Second,  the  error  increases  significantly  as  altitude  is  decreased.  No  signifi¬ 
cant  correlation  was  obtained  between  target  contrast  and  tracker  accuracy. 


r*i 

e? 

m 

.  SILO  LID 

COVER 

2.  CPANE 

3.  EQUIPMENT 
TRAILER 

4.  MAINTENANCE  5.  MINUTEMAN  6.  CUBICLE 

BUILDING  SILO 


10.  STORAGE  SHED 
VENT 


11.  CONSTRUCTION 
BARRICADE 


12.  MISSILE 

TRANSPORTER/ERECTOR 
TIE-DOWN  PLATE 


Figure  57 :  TWELVE  TEST  TARGETS  FROM  1000  FEET  ALTITUDE 
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Raure  58:  MINUTEMAN  SILO  TARGET  AT  FIVE  DIFFERENT  ALTITUDES 
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Figure 59:  TV  TRACKER  SCANNER  SIMULATION  EQUIPMENT 
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ALTITUDE  -  FT. 

Figure  61:  TV  TRACKER  ACCURACY 


4.  Radar  Homing  Sensors.  Several  possible  means  for  using  active  radar 
trackers  for  terminal  guidance  of  homing  ASM's  were  examined,  and  it  was 
concluded  that  the  types  of  targets  for  which  radar  would  make  a  useful  tracking 
device  were  severely  restricted  as  to  proximity  of  adjacent  targets  and  target 
size.  The  target  acquisition  problem  also  imposes  operationally  undesirable 
requirements  on  the  flight  of  the  carrier  aircraft. 

Airborne  acquisition  radars  are  limited  as  to  the  targets  they  can  acquire  and 
track  by  the  size  of  the  antenna,  the  carrier  wavelength,  range  to  the  target, 
and  clutter  background.  The  mode  of  operation  is  such  that  the  operator  in  an 
airplane  watches  a  radar  image  of  the  ground  until  the  intended  target  comes 
into  view,  lines  up  crosshairs  so  as  to  lock  a  tracking  system  to  the  target,  and 
tuen  allows  radar  information  to  be  available  to  a  missile  that  will  use  it  to  home 
in  on  the  ground  target.  The  slow  reaction  time  of  the  human  operator  coupled 
with  the  great  speed  of  the  missile  and  the  small  space  available  on  the  missile 
for  a  tracking  antenna,  dictate  that  the  only  targets  that  can  be  tracked  must  have 
a  radar  cross  section  of  greater  than  1800  square  feet.  A  metal  sphere,  50  feet 
in  diameter,  has  such  a  cross  section.  Such  a  large  target  has  its  compensations 
however,  because  the  tracking  radar  needs  only  a  3-watt  average  radiated  power 
output. 

The  tracking  and  acquisition  problems  are  quite  different,  the  latter  being  much 
more  difficult  to  solve.  Three  acquisition  systems  were  considered:  missile- 
mounted  search  and  track;  airplane- mounted  synthetic  aperture  sc  .x-ch  and 
missile- mounted  track;  and  airplane- mounted  illumination  and  search  and  missile- 
mounted  track.  The  most  promising  is  the  one  in  which  the  missile  tracks  a  patch 
on  the  ground  illuminated  by  the  airplane  until  the  missile  gets  close  enough  to 
convert  to  active  target  seeker  operation.  For  a  target  with  a  radar  cross  section 
larger  than  other  objects  within  a  600-foot-diameter  circle,  an  airplane- mounted 
forward-looking  radar  that  makes  use  of  antenna  beamwidth  for  azimuth  resolution 
ai  1  time  delay  for  range  resolution  requires  the  aircraft  to  fly  within  10  miles  of 
the  *arget.  (This  assumes  a  5-foot-diameter  scanning  antenna  that  operates  on 
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X-band. )  When  the  target  is  located  by  the  observer,  he  will  switch  the  antenna  to 
an  automatic  tracking  mode  that  always  illuminates  the  target.  The  missile  is  then 
launched,  acquires  the  illuminated  patch,  and  continues  toward  the  target  until  its 
resolution  capability  as  an  active  radar  becomes  better  than  as  a  passive  radar.  This 
occurs  at  about  12,000  feet  slant  range  for  an  antenna  diameter  of  1  foot  at  X-band 
(10  GHz).  This  system  has  not  previously  been  proposed  or  constructed.  How¬ 
ever,  the  techniques  involved  are  straightforward  and  do  not  push  the  state  of 
the  art  far  with  respect  to  individual  components. 

If  it  is  necessary  to  launch  the  missile  while  the  aircraft  is  100  miles  from  the 
target,  a  50- foot  antenna  is  required  on  the  aircraft.  This  may  be  implemented 
by  making  a  fixed  array  out  of  the  wing  and  using  an  electronic  scanning  method, 
such  as  a  butler  matrix  discrete  scan,  to  obtain  azimuth  resolution.  This  sys¬ 
tem  does  not  limit  the  maneuvers  of  the  aircraft. 

In  aqy  of  the  acquisition  modes,  the  tracking  system  would  have  a  capability  of 
homing  in  on  jamming  or  homing  in  on  a  radar  installation.  However,  the  option 
would  have  to  be  selected  before  missile  launch. 

For  target  tracking,  the  resolution  capabilities  of  the  radar  Bystem  are  restricted 
by  the  electrical  antenna  length,  range  to  the  target,  and  operating  frequency. 

The  electrical  length  can  be  made  no  larger  than  the  physical  length  unless  special 
techniques  such  as  synthetic  aperture  antennas  are  employed.  These  special  tech¬ 
niques  are,  in  genera!,  side-looking  and  depend  on  the  vehicle  carrying  the  radar 
to  fly  in  a  straight  line.  A  radar  contained  in  a  terminal  guidance  missile  would 
have  to  use  a  conventional  antenna  due  to  its  anomalous  flight  characteristics. 

The  largest  antenna  likely  to  be  installed  in  an  ASM  is  1  foot  in  d'  eter.  The 
smallest  possible  operating  wavelength  is  0.1  foot  or  X-band.  Therefore,  due 
to  diffraction  limitations,  the  smallest  target  that  can  be  identified  at  a  range  of 
20, 000  feet  would  have  an  average  linear  dimension  of  1000  feet.  This  would  not 
present  a  problem  in  tracking  If  the  target  was  isolated  and  had  a  sky  background, 
because  the  resolution  continues  to  improve  as  the  radar  approaches  the  target. 
However,  with  a  surface  background,  it  is  quite  unlikely  that  a  situation  could  be 
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found  whereby  only  one  target  would  be  found  in  a  1000- foot  circle,  and  that  the 
integrated  surface  reflectivity  would  be  less  than  that  of  the  target.  It  is  expected 
that  the  only  missions  that  would  use  a  radar  for  terminal  guidance  would  be  those 
where  visibility  is  restricted,  because  an  optical  system  would  be  capable  of  op¬ 
erating  in  the  presence  of  electronic  countermeasures  and  dense  target  configu¬ 
ration.  The  types  of  targets  for  which  a  radar  tracker  could  be  used  must  be 
isolated  and  of  high  contrast —  such  as  bridges,  trains,  aircraft,  ground  vehicles, 
or  metal  buildings. 

The  radar  cross  section  of  the  background  is  a  function  of  the  type  of  terrain,  rf 
frequency,  and  the  angle  of  incidence,  which  is  the  angle  between  local  vertical 
and  the  line  of  sight  to  the  target.  Multilobe  radars,  such  as  monopulse  or 
conical  scan,  become  confused  as  to  the  exact  position  of  a  target  that  iB  dis¬ 
tributed  over  several  rf  wavelengths.  Large  amounts  of  position  error  can  be 
obtained  at  almost  any  angle  of  incidence.  However,  if  the  targets  are  small, 
or  the  target  has  a  very  rough  jagged  surface,  large  errors  are  more  likely  at 
small  angles  of  incidence. 

The  noise  and  error  sources  are  numerous.  Ten  of  the  most  significant  error 
sources  are:  receiver  thermal  noise  spectrum;  special  emission  from  hot  micro- 
wave  sources;  multipath  noise;  dynamic  lag  in  the  tracking  servo  loop;  glint; 
scintillation;  friction-backlash-stiction  noise;  propagation  noise  error;  radome 
error;  and  change  in  position  of  antenna  beam  center  as  a  function  of  transmitter 
frequency  change. 

In  conclusion,  target  acquisition  presents  difficulties  that  make  the  use  ot  active 
radar  impractical  for  air- to- surface  ml33ile  systems  of  the  type  being  con¬ 
sidered,  except  for  isolated  high  cross-section  targets.  For  specially  selected 
targets  where  radar  is  practical,  it  is  less  susceptible  to  weather  than  a  TV 
tracking  system.  The  noise  sources  in  h  tracking  radar  are  a  function  of  the 
approach  angle  and  nature  of  the  target,  so  that  specific  missions  must  be  con¬ 
sidered  when  applying  the  Kalman  optimal  filtering  technique  to  a  system  with 
a  radar  tracker. 
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APPENDIX  n 


GUIDANCE  COMPARISON  TECHNIQLES 

The  effectiveness  of  proportional,  pursuit,  and  optimal  guidance  ate  compared  in 
Tables  in  and  IV  by  giving  the  standard  (1  o)  deviation  of  the  position  deviation  nor¬ 
mal  to  the  nominal  trajectory  at  the  time  the  nominal  trajectory  hits  the  target. 

In  these  tables  this  quantity  is  called  the  la  miss  distance.  The  use  of  this  quan¬ 
tity  as  the  index  for  comparison  is  motivated  by  the  discussion  in  Section  V.  3. b. 

To  compute  this  quantity  for  each  of  the  guidance  laws ,  the  ASM  equations  of 
motion  for  each  guidance  law  are  linearized  about  the  nominal  trajectory.  Then 
the  differential  equations  for  the  covariance  matrices  of  the  deviations  in  posi¬ 
tion  and  velocity  from  nominal  values  are  written  in  terms  of  the  coefficients  of 
thesp  linearized  equations.  For  each  of  the  guidance  laws,  the  matrix  differen¬ 
tial  equation  for  the  covariance  matrix  is  solved  to  obtain  the  covariance  matrix 
at  the  time  the  nominal  trajectory  hits  the  target.  A  simple  formula  relates 
this  quantity  and  the  la  miss  distance. 

To  reduce  the  complexity  of  the  guidance  comparison  computation  for  proportional 
and  pursuit  guidance ,  the  ASM  m  otion  was  restricted  to  the  plane  of  the  nominal 
trajectory.  For  the  optimal  guidance  law,  the  computation  has  been  carried  out 
for  the  full  6 -degree -of -freedom  system. 

1.  Computation  of  Proportional  Guidance  Covariance  Matrix.  The  equations  of 
a  proportionally  guided  ASM  moving  in  a  plane  are  given  in  Equations  84  through 
91.  Proportional  guidance  attempts  to  make  the  total  normal  force  on  the  missile 
proportional  to  the  target  tracker  rate.  The  coefficient  of  proportionality  is  a 
negative  constant  times  the  nominal  velocity.  The  normal  force  on  the  vehicle 
consists  of  the  aerodynamic  normal  force ,  Cn  ,  plus  the  gravitational  normal 
force.  In  instrumenting  the  proportional  guidance,  the  aerodynamic  normal  force 
is  directly  measured  by  an  accelerometer.  The  normal  gravitational  force  is 
equal  to  the  weight  times  the  cosine  of  the  angle  the  missile  axis  makes  with  a 
horizontal  plane.  This  angle  is  known  except  for  a  random  bias,  b,  from  an 
inertial  platform  measurement.  Hence,  the  normal  force  which  is  commanded, 
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Fjj,  is  given  by  Equation  90.  The  target  tracker  is  assumed  to  be  a  second-order 
linear  tracker. 
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x  =  -  V  cos  y 

(84) 
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Fn  =  -  3  m  V  T) 

(91) 

When  the  ASM  nears  the  target,  the  target  fills  the  TV  sensor  screen  and  the  tar¬ 
get  tracker  is  unable  to  track  the  target.  For  proportional  and  pursuit  guidance 
the  acceleration  normal  to  the  missile  was  held  equal  to  zero  during  this  blind 
zone. 

When  Equations  84  through  91  are  linearized  about  the  nominal  trajectory,  the 
system 
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The  matrices  A(t)  and  B(t)  are  indicated  in  Figures  62  and  63.  In  Figure  62 
the  expression 


+  m  g  sin  (Q  +  y) 


is  abbreviated  by  "Den." 

In  Equation  92  and  Figure  63  the  change  in  angle  of  attack,  Ao (w,  and  speed,  AVW, 
due  to  winds  have  been  related  to  the  components  of  wind  velocity,  Wx  and  W2,  in 
the  horizontal  and  vertical  directions  by  the  formulas 


wy 

Aaw  =  "  v 


AVW  =  Wv 


sin  y\  f  W 


cos  y 


-  sin  y  cos  y)  \W 


Let  P(t)  denote  the  covariance  matrix  of  the  state  variables  of  Equation  92.  The 
winds,  Wx  and  Wz,  and  the  sensor  noise,  N,  aie  considered  as  independent 
white  noise  processes.  It  follows  from  Reference  4,  Page  250,  that  the  covariance 
matrix  P(t)  satisfies  the  differential  equation 

P  =  A(t)  P  +  P  A(t)T  +  B(t)  Q  B(t)T  (96) 

in  which  Q  is  the  covariance  matrix  of  the  noise  vector,  Wx,  Wz,  N. 

2.  Computation  of  Pursuit  Guidance  Covariance  Matrix.  For  pursuit  guidance 
the  equation 


Fn  =  -  3  m  V  r\  (97] 

of  the  proportionally  guided  system  is  replaced  by  the  equation 

Fn  =  -  c3  {y  +  TJ  +  b).  (98] 

Hence,  the  computation  of  miss  distance  for  pursuit  guidance  is  carried  out  in  a 
similar  manner  to  that  of  proportional  guidance.  Equation  92  holds  for  pursuit 
guidance,  if  the  elements  in  the  4^,  5^,  6^,  and  7th  columns  of  the  3r^  and  4th 
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Figure  62:  MATRIX  A  «)  OF  INFLUENCE  COEFFICIENTS  OF  THE  STATE 
VARIABLES  IN  THE  PROPORTIONALLY  GUIDED  STATE 


GUIDED  ASM 


rows  of  the  matrix  A(t)  are  replaced  by  the  elements  of  the  submatrix  indicated 
in  Figure  64. 


3.  Computation  of  the  Optimal  Guidance  Covariance  Matrix.  Let  Si  denote  the 
state  vector  of  Equation  49,  Section  V,  expressed  in  rectangular  coordinates  and 
S2  the  state  vector  expressed  in  flight  path  coordinates,  that  is: 


Equation  49,  Section  V,  implies  Si  and  S2  are  related  by  the  matrix 


Si  = 


f12<»)  52  =  E  52 


(100) 


in  which  3  has  a  3  x  3  identity  matrix  in  the  upper  left  corner  and  the  matrix, 
Fi2<t),  in  the  lower  right  corner  is  the  3x3  submatrix  in  the  upper  right-hand 
corner  of  F(t)  of  Figure  18. 


F(t)  = 


0 

F2l(t) 


F  i2(t) 
F22(l) 


From  the  theory  of  Kalman  filtering: 


E  j(Si  -  Si)  SiTj  =  0. 


(101) 


and  the  upper  left  6x6  submatrix  of  8  x  8  matrix  £  of  Equation  77  is  equal  to 
E  |(Si  -  Si)  (Si  -  Si)Tj.  in  the  remainder  of  this  section  Z  denotes  this  6x6 
submatrix. 


Define  estimates  of  6V,  tfy ,  6x  by: 


(102) 
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Now  comfxite  E  |  S2  S2Tj 

E  Js2  S2T|  =  E  )e_1  Si  SiT  (S'1)Tj  =  E-1  E  )  Si  SiTj  (E'1)T 

A  A 

Let  Si  =  (Si  -  Si)  +  Si,  substitute  into  the  above  equation  nd  simplify  using 
Equation  101 

E  jS2  S2T|  =  S"1  E  )(Si  -  Si)  (Si  -  Si)T|  (E_1)T  +  E'1  E  j  Si  SiT  j  (E“1)T 
From  the  definition  of  E  given  below  Equation  101 

E)  S2  S2Tj  =  E-1  E  (E_1)T  +  E  js2  S2Tj  (103) 


Equation  79  of  Section  V  may  be  rewritten  as: 
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,  b  is  its  estimate, 


N  is  the  noise  vector 


o 


and  K6(t)  is  the  upper  6x2  submatrix  of  8  x  2  matrix  K(t)  in  Equation  76 
of  Section  V . 


Substitute  for  Si  using  Equation  100  noting  that 
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Solve  this  for  S2 
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s2  =  s  1  S  <s2  -  S2) 


/ E  (S2  -  s2)\ 

+  — K6(t)  H(t)(-  -  £  "  A  +  N 


Substitute  from  Equation  53  into  Equation  49  remembering  that  the  optima!  con- 

A 

trol  results  when  the  optimal  estimate,  S2,  is  substituted  for  the  state,  S2. 

So  =  F(t)  S2  -  G(t)  M"1  G(t)T  U(t)  S2  +  D(t)  W 


Now  use  this  to  help  solve  for  S 
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To  put  this  in  terms  of  S2  and  (S2  -  ^2),  add  and  subtract  F(t)  S2  in  the 

/0  Fi2(t)\ 

second  term.  Remember  that  F(t)  -  [  1  and  partition  G  and  D 
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into  3  x  2  submatrices  G(t)  =  ,  D(t)  = 
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S2=S  S(S2-S2)  +  [  l(S2-S2)  +  |  js2 

^210)  F?.2(t)/  \F2i(t)  F22(t)/ 


G(t)M“1G(t)TU(t)S2J 
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Now  collect  terms  to  give: 


S2  «  [ 


F(t)  -  G(t)  M-1(t)  G(t)T  U(t)]  S2  +  D(t)  W 


F21<t)  *12®  +  F22(t)) 


(S2  -  s2) 


+  s  (t)  K6(t)  H(t) 


S(t)  (S2  -  S2) 


A 

b  -  b 


+  N  . 


Denote  the  covariance  matrix  of  S2  by  P(t).  The  quantity  P(t)  is  given  by 
P(t)  =  E  |s2(t)  S2^(t)|  in  which  E  denotes  the  expected  value  operation.  Sup¬ 
pose  the  winds  are  given  by  a  white  noise  process. 


By  computing 


P  A  A  T  A  A  T 

E  [S2(t)  Sjjtt)1  +  S2(t)  S2(t)A 


and  using  the  relationship  that 


E  [(S2<t)  -  S2(t))  A2(t)T]  =  0  and  E  [(b(t)  -  6<t>)  S2<t)T]  =  0,  (109) 


the  equation 

P(t)  =  [F(t)  -  G(t)  M  _1(t)  GT(t)  U(t)]  P(t) 

+  PW  (F(t)  -  G(t)  M-1(t)  G(t)T 

+  E_1(t)  Kg(t)  B  Kg(t)T  S_1(t)T  +  D(t)  E  [w  WT]  D(t)T  (110 

may  be  obtained.  This  equation  may  bo  solved  with  the  initial  condition  P(t0)  =  0 
to  obtain  E  |s2(tl)  ^2(4l)TJ- 

When  thiB  has  b  ien  done,  the  covariance  matrix  E  |s2(t)  S2(t)*^’ J  ic  expressed  in 
terms  of  previously  computed  quantities  b.  Equation  103.. 

4.  Miss  Distance.  The  formula  for  the  mean  square  deviati tu  normal  to  the 
trajectory  is  given  by 


E  [(6  ;tj)  sin  y(tx)  -  6z(t!)  COB  y(tl))2J- 
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This  formula  may  be  written  in  terms  of  the  variances  and  covariances  of  the 
position  variables  as 

sin2  y(t^)  E  [6x(t/]  -  2  sin  y(tj)  cos  yftj)  E  ^Sxftj)  6z(ti)j 

+  cos2  y(tj)  E  £fiz(tj)2J.  (112) 

Hence,  the  mis3  distance  for  proportional,  pursuit,  and  optimal  guidance  maybe 
computed  from  Equation  112  by  substituting  corresponding  elements  of  the  covari¬ 
ance  matrices  computed  in  Sections  1,2,  and  3  of  this  appendix  for  the  correspond¬ 
ing  covariances  in  Equation  112. 

Let  N(t)  denote  the  column  sector  whose  components  are  rero  except  for  the  first 
and  third  components  that  a  e  sin  y(t)  and -cos  y(t),  respectively.  Equation  112 
may  be  written  in  matrix  notation  by 

N(t])T  E  j  S2(t!)  S^Tj  N(tj).  (113) 

Equation  103  implies,  that  for  optimal  guidance, 

N(t)T  E |  S2(t)  S2(t)T)  N(t)  =  N(t)T  H(t)1  £(t)  H(t)-lT  N(t)  +  N(t)T  P(t)  N(t).  (114) 

Because  N(t)  is  zero  except  for  the  first  and  third  elements,  and  the  upper  right- 
hand  3x3  submatrix  of  H(tj)  is  the  identity  matrix,  the  right  side  of  this  equation 
simplifies  to-. 

N(t)T  £(t)  N(t)  +  N(t)T  P(t)  N(t)  (115) 

The  matrix  L  is  the  covariance  matrix  of  the  difference  between  the  ASM  state  vari¬ 
ables  and  their  estimates.  This  matrix  is  the  solution^  Equation  77  that  was 
solved  numerically  to  determine  the  weighting  coefficients  of  the  optimal  filter 
that  are  given  by  Equation  76.  Notice  that  Equation  77  does  not  depend  on  the 
coefficient  matrix  M  GTU  of  the  optimal  feedback  controller.  The  covariance 
matrix  £  is  an  Indicator  of  how  well  the  optimal  filter  is  estimating  the  >’>i  state 
variables.  Thus  the  first  term  of  Equation  115  can  be  considered  as  the  contri¬ 
bution  of  the  error  in  estimating  the  ASM  state  variables  to  miss  distance.  This 
term  of  Equation  115  will  be  called  filter  error. 

(1)  Modified  as  described  on  page  170. 
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Equation  110  for  the  matrix  P(t)  of  the  second  term  of  Equation  115  does  involve 
the  coefficients  of  the  matrix  G^U  of  the  optimal  feedback  controller.  Note 

that  if  the  optimal  feedback  controller  M-1  G^U  S2  was  replaced  by  a  linear 

.  A 

controller  of  the  form,  A  S2,  the  discussion  of  Section  3  of  this  appendix  would 
be  unchanged  except  that  the  matrix  A  would  replace  M-1  G^U  in  the  formulas. 
Thus  the  second  term  in  Equation  115  may  be  considered  as  the  contribution  of 
the  feedback  controller  to  the  impact  error.  This  term  will  be  called  controller 
error . 

These  two  types  of  error  have  been  plotted  as  functions  of  time  in  Figure  27  as 
discussed  previously  to  show  relative  size  of  these  errors  at  various  times  from 
target  acquisition  to  the  time  the  nominal  trajectory  hits  the  target. 

5 .  Computation  Techniques.  To  compute  the  coefficients  of  the  optimal  filter 
and  optimal  feedback  control  law  and  to  compare  miss  distances  for  optimal, 
proportional,  and  pursuit  guidance,  the  matrix  differential  Equations  54,  77,  96, 
and  110  were  solved  numerically.  Nominal  trajectories  and  the  influence  coeffi¬ 
cient  matrices,  F(t)  and  G(t),  defined  in  Equation  39  were  computed  using  exist¬ 
ing  Boeing  computer  programs.  The  matrices,  F(t)  and  G(t),  which  are  used 
as  coefficients  in  Equations  54,  96,  and  110,  were  fed  as  punched  input  to  the 
programs  for  solving  these  equations.  In  solving  Equation  110,  the  solution 
matrices,  U(t)  and  E(t),  of  Equations  54  and  77  were  also  punched  inputs  to 
the  computer  program. 

Equations  54,  77,  96,  and  110  were  solved  on  a  Univac  1107  computer.  Fortran 
IV  was  used  as  a  computing  language  in  coding  the  program.  Each  of  the  pro¬ 
grams  was  coded  using  less  than  two  hundred  instructions.  A  Romberg  integra¬ 
tion  technique  was  used  in  Integrating  the  equations.  This  integration  technique 
uses  a  global  integrator  and  has  an  adaptive  method  for  varying  integration  step 
size  based  on  accuracy  estimates  of  the  components  of  the  matrix  being  computed. 

These  accuracy  estimates  were  printed  out  to  check  integration  accuracy  and  to 
aid  in  monitoring  the  program.  Several  of  the  computations  were  also  carried 
out  using  a  Runge-Kutta  variable  step  integration  technique.  The  results  of  the 
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computations  with  the  two  dif^rent  methods  were  compared  as  a  check  on  the 
integration  techniques. 

Equations  54,  77,  96,  and  110  are  matrix  differential  equations  of  dimensions 
6x6,  8x8,  7x7,  and  6x6  respectively.  Thus  in  each  of  the  computations 
either  36,  49,  or  64  quantities  are  being  computed  as  functions  of  time.  Integra¬ 
tion  times  for  these  equations  on  the  1107  computer  varied  between  3  and  12 
minutes. 


APPENDIX  in 

OBSERVABILITY  PROBLEMS  IN  THE  OPTIMAL -GUIDANCE  FILTER 


Because  Figures  28  and  29  show  that  the  component  of  position  deviation  normal 
to  the  nominal  trajectory  is  estimated  very  accurately  while  the  vertical  and  hori¬ 
zontal  components  of  position  derivation  are  estimated  much  less  accurately,  it  is 
necessary  to  analyze  the  observation  process  to  understand  the  reasons  for  this 
phenomena  and  to  make  sure  it  does  not  cause  poor  system  performance. 

Insight  about  the  measurement  process  may  be  gained  by  examining  the  effect  of 
linearizing  the  line-of-sight  angle  measurement  equations.  If  the  random  biases, 
bj  and  b2,  and  the  noise,  Nj  and  N2,  were  set  equal  to  their  mean  values  of  zero, 
the  measurement.  Equations  18  and  19,  become 


4>  =  tan  1  ^ 

X 

(116) 

(117) 

Linearizing  these  equations  it  is  seen  that 

6#  =  Ax(t)  fix  +  A2(t)  fiy 

(118) 

fi*  =  B1(t)  fix  +  B2(t)  fiy  +  B3(t)  6z 

(119) 

in  which  Aj(t),  etc. ,  are  the  partial  derivatives  of  the  respective  expressions  on 

the  right  side  of  the  above  equations  evaluated  on  the  nominal  trajectory. 


Now 


Ai(t)  =  - 


(x2  +  y2)  ’ 


A2(t)  = 


X2  +  y2 


Bi(t)  = 


z  x 


(S2  ♦  y2)172  *  y2  .  s2) 


.  B2(t)  = 


z_  y 


(x2  +  y2)1//2  (x2  +  y2  +  z2) 


(120) 


B3(t)  =  - 


Yx2  +  y2‘ 


x2  +  y2  +  z2 
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Now  defining 


cos  ?  = 


V  +  y2 


sin? 


Y  +  9Z 


cos 


j  V  x2  +  y2 

V  X2  +  y2  +  z2 
r  =Vx2  +  y2  +  z2 


sin  $  = 


Y  x2  +  y2  +  z2 


these  formulas  become 


Ai(t) 


sin 1 


r  cos 


$ 


A2(t)  =  r 


cos 


* 


C08  $ 


(121) 


Bl(t)  =  -  4  cos  0  sin  5  t  B2(t)  =  -  4  sin  ?  sin  ®  (122) 

B3(t)  =  4  cos  $  . 

The  vectors 

(-  sin?,  cos  0 ,  0) 

_  -  (123) 

(-  cos  ?  sih  -  sin  ?  sin  $,  cos  $) 

are  mutually  perpendicular  and  are  perpendicular  to  the  vector 

(cos?  cos  $,  sin?  cos  sin  $>)  (124) 

that  points  along  the  line  of  sight  of  the  nominal  trajectory.  The  linearized  meas¬ 
urement  equations  imply  that  the  variations  of  missile  positions  in  these  two  per¬ 
pendicular  directions  times  the  respective  weightings  y^r  cos4>  )and  1/r  are 
being  measured.  The  component  of  position  along  the  nominal  trajectory  is  not 
being  measured.  Hence,  the  only  way  in  which  the  filter  can  correct  errors  in 
the  estimate  of  position  along  the  nominal  trajectory  is  through  the  information 
given  by  these  two  perpendicular  measurements  on  position  along  the  trajectory. 

If  there  is  bending  of  the  nominal  trajectory,  the  three  directions  expressed  by 
the  above  vectors  will  change  along  the  nominal  trajectory.  From  this  it  should 
be  expected  that  the  perpendicular  measurements  will  furnish  a  small  amount  of 


180 


Information  on  position  along  the  nominal.  It  should  be  expected  that  the  correc¬ 
tion  by  the  filter  of  the  estimate  of  this  component  of  position  will  be  slow.  For¬ 
tunately,  the  components  of  position  deviation  perpendicular  to  the  nominal  tra¬ 
jectory  are  the  important  ones  in  defining  miss  distance. 

In  the  discussion  given  below,  it  will  be  shown  that  if  the  nominal  trajectory  is 
a  straight  line,  the  system  is  not  completely  observable  according  to  the  cri¬ 
teria  defined  by  Kalman  (Reference  4,  Page  36).  For  simplicity  of  exposition, 
suppose  the  nominal  trajectory  is  a  straight  line  to  the  target  making  a  45-degree 
angle  with  the  ground,  given  by  the  equations: 


*<t>  =  t  -  tx 

y<t>  =  o 

z(t)  =  tx  -  t 

Evaluating  the  Equations  120  for  Aj(t),  etc. ,  in  this  case  gives 

A1«  -  «  Bi«  “  FTP! j- 

A2«  “  tfh)  B2«  ■  0 

B3W  =  2  (tx  -  t) 


(125) 


(L26) 


The  optimal  filter.  Equation  79,  is  based  on  the  consideration  of  Equations  74  as 
state  equations  and  a  modification  of  Equations  59  as  the  measurement  equations. 
The  system  (Equation  74)  with  measurements  (Equation  59)  would  be  completely 
observable  in  the  terminology  of  Reference  4,  Page  36,  if  the  matrix 

B 

M(0,  s)  =  J  $>T(t,  ti)  HT(t)  R-1(t)  H(t)  $(t,  t!)  dt  (127) 


is  nonsingular  for  some  value  of  s  between  0  and  tj.  The  quantity  <t>  (t,  T)  is  the 
Impulse  response  of  the  system  (Equation  74).  Equation  74  may  be  integrated  to 
show  that; 
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1 


0  0  t  -  T  0 


0 


0 


0 


Consider  the  case  when  the  covariance  matrix,  R(t),  is  a  scalar  times  the  identity 
matrix.  Then  R"1  is  a  scalar  times  the  identity  matrix  and  it  may  be  factored  out 
of  the  expression  for  M(0,  s). 


or 


I 


*T(t,  tj)  HT(t)  = 


For  brevity,  denote  —  (t-tj)  by  T .  Now 


$T(t,  tx)  HT(t)  H(t)  <t>( t,  t!) 


Tc 

0 


?MiT 


fo-tif1 

1 

2 


1 

2 

0 

1 


2T  ° 

0  2 T 

2T  ° 
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From  Equation  131  it  is  easily  seen  that  M(0,  s)  is  singular.  Hence,  the  system 
le  not  completely  observable. 


183 


(131) 


Security  Classification  _ 

DOCUMENT  CONTROL  DATA  •  RAD  ———————— 

fSsOurtty  classification  of  Utlm,  body  of  abstract  and  Indexing  annotation  mu  at  bv  aalarad  wttan  (ht  ortral/  raporl  fa  c  Uamitiad) 


J3  ABSTRACT 

Design  guidelines  were  developed  to  provide  a  basis  for  conducting  design  trades 
for  a  hasdng  type  air- to- surface  slsslle  vlth  high  terelnal  accuracy.  Three 
basic  hoeing  guidance  concepts,  proportional,  pursuit,  and  opt  leal  guidance  ware 
evaluated  on  the  basis  of  1  upset  error.  Two  noalnal  trajectories  were  investigat¬ 
ed.  An  opt  leal  guidance  lor  was  developed  for  an  A8M  with  realistic  aerodynaaie 
and  sensor  characteristics .  This  guidance  las  was  baaed  upon  the  use  of  a  Hainan 
filter  to  obtain  heat  estimates  of  the  ASM  state  variable  errors,  and  a  control 
concept  which  ainlnizes  the  sub  of  the  seen  square  lapeet  error  and  the  integral 
of  a  quadratic  fora  of  the  autopilot  control  variables.  A  linearised  differential 
equation  progrtn  which  oonputed  the  nean  square  lapeet  error  In  the  foxa  of  a 
covariance  aatrix  deviation  perpendicular  to  the  ntartmri  trajectory,  was  used  for 
coaparlson  of  the  guidance  laws.  A  noraal  acceleration  autopilot  was  designed  to 
aeet  the  mission  requirements,  and  advanced  bistable  controller  techniques  were 
applied  to  obtain  a  quasl-sdaptlve  autopilot  which  required  no  gain  changes 
throughout  the  ASM  aldcouree  and  temlnal  phases.  A  Halted  state-of-the-art 
survey  was  endue  ted  of  bearing  and  Inertial  eensore,  and  on-board  digital 
conputers  suitable  for  a  hosting  ASM. 


1  OniGIHATIM  C  activity  (Coipo/»lo  Author; 

The  Boeing  Conpany,  Aero-Space  Group 
P.  0.  Boot  3707 
Seattle,  Wash.  9812  V 


DD  1473 


Security  Classification 


Uncl*e>ifled 


Security  Classification 


Opfciam  control 
Guidance  and  control 
Alr-to-aurface  missiles 


INSTRUCTIONS 


L  ORIGINATING  ACTIVITY:  Enter  the  ntm«  end  addreee 
of  the  contractor,  subcontractor,  grantee.  Department  of  De» 
fenee  activity  or  other  organisation  f corporate  mu  (hot)  issuing 
the  report. 

2e.  REPORT  SECUMTY  CLASSIFICATION:  Enter  'be  over¬ 
all  security  classification  of  the  report.  Indicate  whether 
“Restricted  Data"  it  included  Marie  log  is  to  he  in  accord- 
ance  with  appropriate  security  regulations. 

2b.  GROUP:  Automatic  downgrading  is  specif/ed  in  DoD  Di¬ 
rective  5200.  10  and  Armed  Forces  Industrial  Manual.  Enter 
the  gtoop  comber.  Also,  when  applicable,  ahow  that  optional 
markings  have  been  used  for  Group  3  and  Group  4  at  author¬ 
ized. 

3.  REPORT  TITLE:  Enter  the  complete  report  title  in  all 
capital  letters.  Titles  in  all  cases  should  be  unclassified. 

If  a  meaningful  title  cannot  be  selected  without  classifica¬ 
tion.  show  title  classification  in  all  capitals  in  parenthesis 
immediately  following  tha  title. 

4.  DESCRIPTIVE  NOTES:  If  appropriate,  enter  the  type  of 
report,  e.g,,  interim,  progress,  summary,  annual,  or  final. 

Give  the  inclusive  dates  when  a  specific  reporting  period  is 
covered. 

5.  AUTHOR(S):  Enter  the  name<s)  of  authors)  as  shown  on 
or  in  the  report.  Enter  last  name,  first  name,  middle  initial. 

It  military,  stow  rank  end  branch  of  service.  The  name  of 
the  ftriittipal  author  iu  an  absolute  minimum  requirement. 

6.  REPORT  DATE:  Enter  the  date  of  the  report  es  day, 
month,  yean  or  month,  year.  If  more  than  one  date  appears 
on  the  report.  Use  data  of  publication. 

7a.  TOTAL  NUMBER  OF  PAGES:  The  total  page  count 
should  follow  normal  pagination  procedures,,  he.,  enter  the 
member  of  pages  containing  Information 
76.  NUMBER  OF  REFERENCED  Enter  'be  total  number  of 
references  cited  ir.  the  report. 

8a.  CONTRACT  OR  GRANT  HVUBOt  If  appropriate,  ester 
the  applicable  number  of  the  contract  or  grant  under  which 
the  report  was  written. 

86,  Sc,  fc  Id.  PROJECT  NUMBER:  Enter  the  appropriate 
military  dspretmant^deotlflcatloa,  such  as  project  amber, 
subproject  number,  ays. am  cumbers,  task  number,  etc, 

9a.  ORIGINATOR**  REPORT  NUWER<8):  Enter  the  offi¬ 
cial  report  number  by  which  the  document  will  be  idetSi/ied 
at.d  controlled  by  the  originating  activity.  This  number  must 
be  unique  to  this  report. 

96.  OTHER  REPORT  NUMBER**):  If  the  report  has  bee* 
assigned  any  other  report  numbers  ( either  6y  <6#  originator 
or  by  thm  mponmor),  s!eo  enter  this  number* a). 

10.  /  VAIL ABILrr Y/LIMITATION  NOTICED  Enter  any  lim¬ 
itations  on  further  dissemination  of  the  report,  other  than  those 


imposed  by  security  classification,  using  standard  statements 
such  as:  » 

(1)  “Qualified  requesters  may  obtain  copies  of  this 
report  from  DDG  ’* 

(2)  "Foreign  announcement  and  dissemination  of  this 
report  by  DDC  ia  not  authorized  ” 

(3)  "U.  S.  Government  agencies  may  obtain  copies  of 
this  report  directly  from  DDC.  Other  qualified  DDC 
users  shall  request  through 


(4)  **U.  S.  military  agencies  may  obtain  copies  of  this 

report  directly  from  DDC.  Other  qualified  users 
shall  request  through 


(5)  “All  distribution  of  this  report  is  controlled.  Qual¬ 
ified  DDC  users  shall  request  through 


If  the  rsport  hss  be  on  firniahed  to  the  Office  of  Technical 
Services,  Department  of  Commerce,  for  aale  to  the  public,  indi¬ 
cate  this  feet  and  enter  the  price.  If  known 

IL  SUPPLEMENTARY  NOTES:  Use  for  additional  explana¬ 
tory  notes. 

IX  SPONSORING  MILITARY  ACTIVITY:  Enter  the  name  of 
the  departmental  project  office  or  laboratory  sponsoring  fpey- 
ing  for )  the  research  and  devalopmeat.  Include  adrfcaaa. 

13  ABSTRACT:  Enter  an  abstract  giving  a  brief  and  factual 
summary  of  the  document  indicative  of  the  report,  even  though 
it  may  also  appear  elsewhere  ia  the  body  of  the  technical  re¬ 
port.  If  additional  specs  la  required,  a  continuation  sheet  shell 
be  attached. 

It  is  highly  desirable  that  the  abstract  of  classified  reports 
be  unclassified.  Each  paragraph  of  the  abstract  shall  end  with 
an  indication  of  the  military  security  classification  of  the  in¬ 
formation  in  the  paragraph,  represented  as  fT5>.  (S).  (C),  or  (U). 

There  ia  no  limitation  op  the  length  of  the  abstract.  How¬ 
ever,  the  suggested  length  is  from  150  to  225  words. 

14.  KEY  FORDS:  Kay  words  are  technically  meaningful  terns 
or  short  phrases  that  characterise  s  report  and  may  be  used  as 
index  entries  for  cataloging  the  report.  Kay  words  must  be 
selected  so  that  no  security  classification  is  required-  Identi¬ 
fiers,  such  as  equipment  model  deeination.  trade  name,  military 
project  code  name,  geographic  location,  may  be  used  as  key 
words  but  will  be  followed  by  an  indication  of  technical  con¬ 
text.  The  assignment  of  links,  rulas,  and  weights  is  optional. 


Security  Classification 


